Allele frequency spectrum of known ankylosing spondylitis associated variants in a Swedish population.

Objective: The genetic predisposition to ankylosing spondylitis (AS) has been most widely studied in cohorts with European ancestry. However, within Europe, disease prevalence is higher in Sweden. Given this, we aimed to characterize known AS susceptibility variants in a homogeneous Swedish data set, assessing reproducibility and direction of effect.Method: The power to detect association within an existing Swedish targeted sequencing study (381 controls; 310 AS cases) was examined, and a set of published associations (n = 151) was intersected with available genotypes. Association to disease was calculated using logistic regression accounting for population structure, and HLA-B27 status was determined with direct polymerase chain reaction genotyping.Results: The cases were found to be 92.3% HLA-B27 positive, with the data set showing ≥ 80% predictive power to replicate associations, with odds ratios ≥ 1.6 over a range of allele frequencies (0.1-0.7). Thirty-four markers, representing 23 gene loci, were available for investigation. The replicated variants tagged MICA and IL23R loci (p < 1.47 × 10-3), with variable direction of effect noted for gene loci IL1R1 and MST1.Conclusion: The Swedish data set successfully replicated both major histocompatibility complex (MHC) and non-MHC loci, and revealed a different replication pattern compared to discovery data sets. This was possibly due to population demographics, including HLA-B27 frequency and measured comorbidities.

Roan, ticked and clear coat patterns in the canine are associated with three haplotypes near usherin on CFA38.

White coat patterning is a feature of many dog breeds and is known to be coded primarily by the gene micropthalmia-associated transcription factor (MITF). This patterning in the coat can be modified by other factors to produce the attractive phenotypes termed 'ticked' and 'roan' that describe the presence of flecks of color that vary in distribution and intensity within otherwise 'clear' white markings. The appearance of the pigment in the white patterning caused by ticking and roaning intensifies in the weeks after birth. We applied genome-wide association to compare English Cocker Spaniels of roan phenotype (N = 34) with parti-color (non-roan) English Cocker Spaniels (N = 9) and identified an associated locus on CFA 38, CFA38:11 057 040 (Praw  = 8.9 × 10-10 , Pgenome  = 2.7 × 10-5 ). A local case-control association in English Springer Spaniels comparing 11 ticked and six clear dogs identified indicative association with a different haplotype, CFA38:11 122 467G>T (Praw  = 1.7 × 10-5 ) and CFA38:11 124 294A>C (Praw  = 1.7 × 10-5 ). We characterize three haplotypes in Spaniels according to their putative functional variant profiles at CFA38:11 111 286C>T (missense), CFA38:11 131 841-11 143 239DUP.insTTAA (using strongly linked marker CFA38:11 143 243C>T) and CFA38:11 156 425T>C (splice site). In Spaniels, the haplotypes work as an allelic series including alleles (t, recessive clear; T, dominant ticked/parti-color; and TR , incomplete dominant roan) to control the appearance of pigmented spots or flecks in otherwise white areas of the canine coat. In Spaniels the associated haplotypes are t (CCT), T (TCC) and TR (TTT) for SNP markers on CFA38 at 11 111 286C>T, 11 143 243C>T and 11 156 425T>C respectively. It is likely that other alleles exist in this series and together the haplotypes result in a complex range of patterning that is only visible when dogs have white patterning resulting from the epistatic gene Micropthalmia-associated transcription factor (the S-locus).

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.

Presence of thoracic and lumbar vertebral malformations in pugs with and without chronic neurological deficits.

Congenital vertebral malformations (CVMs) are common in brachycephalic dogs such as the pug, and are often considered incidental findings. However, specific CVMs have been suggested to be associated with neurological deficits in pugs. The objective of this study was to investigate the clinical importance of CVMs in the pug by comparing computed tomography studies of the thoracolumbar spine from pugs without neurological deficits with those from pugs with a confirmed T3-L3 spinal cord lesion and neurological deficits consistent with a chronic T3-L3 myelopathy. A total of 57 pugs were recruited into the study from Sweden (n=33), United Kingdom (n=21) and Norway (n=3); 30 with neurological deficits and 27 without. Focal T3-L3 pathology was confirmed in all pugs with neurological deficits by magnetic resonance imaging (n=29) and/or pathology (n=15). Computed tomography studies of the thoracolumbar spine from pugs with and without neurological deficits were compared to investigate possible associations between presentation of neurological deficits consistent with chronic T3-L3 pathology and signalment variables, presence of CVMs and type of CVMs. Congenital vertebral malformations were as common in pugs with, as in pugs without, neurological deficits. Regardless of neurological status, the majority of pugs (96%) presented with one or more CVM. An association between presence, or type of CVM in the T1-L3 vertebral column, and neurological deficits consistent with T3-L3 pathology could not be confirmed.

Interactions between CXCR4 and CXCL12 promote cell migration and invasion of canine hemangiosarcoma.

The CXCR4/CXCL12 axis plays an important role in cell locomotion and metastasis in many cancers. In this study, we hypothesized that the CXCR4/CXCL12 axis promotes migration and invasion of canine hemangiosarcoma (HSA) cells. Transcriptomic analysis across 12 HSA cell lines and 58 HSA whole tumour tissues identified heterogeneous expression of CXCR4 and CXCL12, which was associated with cell movement. In vitro, CXCL12 promoted calcium mobilization, cell migration and invasion that were directly proportional to surface expression of CXCR4; furthermore, these responses proved sensitive to the CXCR4 antagonist, AMD3100, in HSA cell lines. These results indicate that CXCL12 potentiates migration and invasion of canine HSA cells through CXCR4 signalling. The direct relationship between these responses in HSA cells suggests that the CXCR4/CXCL12 axis contributes to HSA progression.

Extended exome sequencing identifies BACH2 as a novel major risk locus for Addison's disease.

BACKGROUND: Autoimmune disease is one of the leading causes of morbidity and mortality worldwide. In Addison's disease, the adrenal glands are targeted by destructive autoimmunity. Despite being the most common cause of primary adrenal failure, little is known about its aetiology. METHODS: To understand the genetic background of Addison's disease, we utilized the extensively characterized patients of the Swedish Addison Registry. We developed an extended exome capture array comprising a selected set of 1853 genes and their potential regulatory elements, for the purpose of sequencing 479 patients with Addison's disease and 1394 controls. RESULTS: We identified BACH2 (rs62408233-A, OR = 2.01 (1.71-2.37), P = 1.66 × 10-15 , MAF 0.46/0.29 in cases/controls) as a novel gene associated with Addison's disease development. We also confirmed the previously known associations with the HLA complex. CONCLUSION: Whilst BACH2 has been previously reported to associate with organ-specific autoimmune diseases co-inherited with Addison's disease, we have identified BACH2 as a major risk locus in Addison's disease, independent of concomitant autoimmune diseases. Our results may enable future research towards preventive disease treatment.

Absolute quantification reveals the stable transmission of a high copy number variant linked to autoinflammatory disease.

BACKGROUND: Dissecting the role copy number variants (CNVs) play in disease pathogenesis is directly reliant on accurate methods for quantification. The Shar-Pei dog breed is predisposed to a complex autoinflammatory disease with numerous clinical manifestations. One such sign, recurrent fever, was previously shown to be significantly associated with a novel, but unstable CNV (CNV_16.1). Droplet digital PCR (ddPCR) offers a new mechanism for CNV detection via absolute quantification with the promise of added precision and reliability. The aim of this study was to evaluate ddPCR in relation to quantitative PCR (qPCR) and to assess the suitability of the favoured method as a genetic test for Shar-Pei Autoinflammatory Disease (SPAID). RESULTS: One hundred and ninety-six individuals were assayed using both PCR methods at two CNV positions (CNV_14.3 and CNV_16.1). The digital method revealed a striking result. The CNVs did not follow a continuum of alleles as previously reported, rather the alleles were stable and pedigree analysis showed they adhered to Mendelian segregation. Subsequent analysis of ddPCR case/control data confirmed that both CNVs remained significantly associated with the subphenotype of fever, but also to the encompassing SPAID complex (p < 0.001). In addition, harbouring CNV_16.1 allele five (CNV_16.1|5) resulted in a four-fold increase in the odds for SPAID (p < 0.001). The inclusion of a genetic marker for CNV_16.1 in a genome-wide association test revealed that this variant explained 9.7 % of genetic variance and 25.8 % of the additive genetic heritability of this autoinflammatory disease. CONCLUSIONS: This data shows the utility of the ddPCR method to resolve cryptic copy number inheritance patterns and so open avenues of genetic testing. In its current form, the ddPCR test presented here could be used in canine breeding to reduce the number of homozygote CNV_16.1|5 individuals and thereby to reduce the prevalence of disease in this breed.

Effect of Breed on Plasma Endothelin-1 Concentration, Plasma Renin Activity, and Serum Cortisol Concentration in Healthy Dogs.

BACKGROUND: There are breed differences in several blood variables in healthy dogs. OBJECTIVE: Investigate breed variation in plasma endothelin-1 (ET-1) concentration, plasma renin activity, and serum cortisol concentration. ANIMALS: Five-hundred and thirty-one healthy dogs of 9 breeds examined at 5 centers (2-4 breeds/center). METHODS: Prospective observational study. Circulating concentrations of ET-1 and cortisol, and renin activity, were measured using commercially available assays. Absence of organ-related or systemic disease was ensured by thorough clinical investigations, including blood pressure measurement, echocardiography, ECG, blood and urine analysis. RESULTS: Median ET-1 concentration was 1.29 (interquartile range [IQR], 0.97-1.82) pg/mL, median cortisol concentration 46.0 (IQR, 29.0-80.8) nmol/L, and median renin activity 0.73 (IQR, 0.48-1.10) ng/mL/h in all dogs. Overall, breed differences were found in ET-1 and cortisol concentrations, and renin activity (P < .0001 for all). Pair-wise comparisons between breeds differed in 67% of comparisons for ET-1, 22% for cortisol, and 19% for renin activity, respectively. Within centers, breed differences were found at 5/5 centers for ET-1, 4/5 centers for cortisol, and 2/5 centers for renin activity. Newfoundlands had highest median ET-1 concentration, 3 times higher than Cavalier King Charles Spaniels, Doberman Pinschers, and Dachshunds. Median renin activity was highest in Dachshunds, twice the median value in Newfoundlands and Boxers. Median cortisol concentration was highest in Finnish Lapphunds, almost 3 times higher than in Boxers. CONCLUSIONS AND CLINICAL IMPORTANCE: Breed variation might be important to take into consideration when interpreting test results in clinical studies.

Breed distribution of SOD1 alleles previously associated with canine degenerative myelopathy.

BACKGROUND: Previous reports associated 2 mutant SOD1 alleles (SOD1:c.118A and SOD1:c.52T) with degenerative myelopathy in 6 canine breeds. The distribution of these alleles in other breeds has not been reported. OBJECTIVE: To describe the distribution of SOD1:c.118A and SOD1:c.52T in 222 breeds. ANIMALS: DNA from 33,747 dogs was genotyped at SOD1:c.118, SOD1:c.52, or both. Spinal cord sections from 249 of these dogs were examined. METHODS: Retrospective analysis of 35,359 previously determined genotypes at SOD1:c.118G>A or SOD1:c.52A>T and prospective survey to update the clinical status of a subset of dogs from which samples were obtained with a relatively low ascertainment bias. RESULTS: The SOD1:c.118A allele was found in cross-bred dogs and in 124 different canine breeds whereas the SOD1:c.52T allele was only found in Bernese Mountain Dogs. Most of the dogs with histopathologically confirmed degenerative myelopathy were SOD1:c.118A homozygotes, but 8 dogs with histopathologically confirmed degenerative myelopathy were SOD1:c.118A/G heterozygotes and had no other sequence variants in their SOD1 amino acid coding regions. The updated clinical conditions of dogs from which samples were obtained with a relatively low ascertainment bias suggest that SOD1:c.118A homozygotes are at a much higher risk of developing degenerative myelopathy than are SOD1:c.118A/G heterozygotes. CONCLUSIONS AND CLINICAL IMPORTANCE: We conclude that the SOD1:c.118A allele is widespread and common among privately owned dogs whereas the SOD1:c.52T allele is rare and appears to be limited to Bernese Mountain Dogs. We also conclude that breeding to avoid the production of SOD1:c.118A homozygotes is a rational strategy.

Breed differences in natriuretic peptides in healthy dogs.

BACKGROUND: Measurement of plasma concentration of natriuretic peptides (NPs) is suggested to be of value in diagnosis of cardiac disease in dogs, but many factors other than cardiac status may influence their concentrations. Dog breed potentially is 1 such factor. OBJECTIVE: To investigate breed variation in plasma concentrations of pro-atrial natriuretic peptide 31-67 (proANP 31-67) and N-terminal B-type natriuretic peptide (NT-proBNP) in healthy dogs. ANIMALS: 535 healthy, privately owned dogs of 9 breeds were examined at 5 centers as part of the European Union (EU) LUPA project. METHODS: Absence of cardiovascular disease or other clinically relevant organ-related or systemic disease was ensured by thorough clinical investigation. Plasma concentrations of proANP 31-67 and NT-proBNP were measured by commercially available ELISA assays. RESULTS: Overall significant breed differences were found in proANP 31-67 (P < .0001) and NT-proBNP (P < .0001) concentrations. Pair-wise comparisons between breeds differed in approximately 50% of comparisons for proANP 31-67 as well as NT-proBNP concentrations, both when including all centers and within each center. Interquartile range was large for many breeds, especially for NT-proBNP. Among included breeds, Labrador Retrievers and Newfoundlands had highest median NT-proBNP concentrations with concentrations 3 times as high as those of Dachshunds. German Shepherds and Cavalier King Charles Spaniels had the highest median proANP 31-67 concentrations, twice the median concentration in Doberman Pinschers. CONCLUSIONS AND CLINICAL IMPORTANCE: Considerable interbreed variation in plasma NP concentrations was found in healthy dogs. Intrabreed variation was large in several breeds, especially for NT-proBNP. Additional studies are needed to establish breed-specific reference ranges.

Molecular profiling reveals prognostically significant subtypes of canine lymphoma.

We performed genomewide gene expression analysis of 35 samples representing 6 common histologic subtypes of canine lymphoma and bioinformatics analyses to define their molecular characteristics. Three major groups were defined on the basis of gene expression profiles: (1) low-grade T-cell lymphoma, composed entirely by T-zone lymphoma; (2) high-grade T-cell lymphoma, consisting of lymphoblastic T-cell lymphoma and peripheral T-cell lymphoma not otherwise specified; and (3) B-cell lymphoma, consisting of marginal B-cell lymphoma, diffuse large B-cell lymphoma, and Burkitt lymphoma. Interspecies comparative analyses of gene expression profiles also showed that marginal B-cell lymphoma and diffuse large B-cell lymphoma in dogs and humans might represent a continuum of disease with similar drivers. The classification of these diverse tumors into 3 subgroups was prognostically significant, as the groups were directly correlated with event-free survival. Finally, we developed a benchtop diagnostic test based on expression of 4 genes that can robustly classify canine lymphomas into one of these 3 subgroups, enabling a direct clinical application for our results.

Diabetes mellitus in elkhounds is associated with diestrus and pregnancy.

BACKGROUND: Female Elkhounds are shown to be at increased risk for diabetes mellitus, and occurrence of diabetes during pregnancy has been described in several cases. HYPOTHESIS: Onset of diabetes mellitus in Elkhounds is associated with diestrus. ANIMALS: Sixty-three Elkhounds with diabetes mellitus and 26 healthy controls. METHODS: Medical records from 63 Elkhounds with diabetes were reviewed and owners were contacted for follow-up information. Blood samples from the day of diagnosis were available for 26 dogs. Glucose, fructosamine, C-peptide, growth hormone (GH), insulin-like growth factor-1, progesterone, and glutamate decarboxylase isoform 65-autoantibodies were analyzed and compared with 26 healthy dogs. Logistic models were used to evaluate the association of clinical variables with the probability of diabetes and with permanent diabetes mellitus after ovariohysterectomy (OHE). RESULTS: All dogs in the study were intact females and 7 dogs (11%) were pregnant at diagnosis. The 1st clinical signs of diabetes mellitus occurred at a median of 30 days (interquartile range [IQR], 3-45) after estrus, and diagnosis was made at a median of 46 days (IQR, 27-62) after estrus. Diabetes was associated with higher concentrations of GH and lower concentrations of progesterone compared with controls matched for time after estrus. Forty-six percent of dogs that underwent OHE recovered from diabetes with a lower probability of remission in dogs with higher glucose concentrations (odds ratio [OR], 1.2; P=.03) at diagnosis and longer time (weeks) from diagnosis to surgery (OR, 1.5; P=.05). CONCLUSIONS: Diabetes mellitus in Elkhounds develops mainly during diestrus and pregnancy. Immediate OHE improves the prognosis for remission of diabetes.

Association of dopamine- and serotonin-related genes with canine aggression.

Human-directed canine aggression was studied using 50 aggressive and 81 non-aggressive dogs. We examined 62 single nucleotide polymorphisms (SNPs) occurring in or in the close vicinity of 16 neurotransmitter-related genes. Allelic associations with aggression were identified for DRD1, HTR1D, HTR2C and SLC6A1. Risk or protective haplotypes for aggressive behaviour based on 2-5 SNPs were identified. The frequency of aggressive dogs varied significantly between the haplotypes within loci and the odds ratios of aggression in dogs with risk haplotypes compared with protective haplotypes varied from 4.4 (HTR2C) to 9.0 (SLC6A1). A risk haplotype across the neurotransmitter receptor gene HTR1D harboured a non-synonymous SNP with a potential effect on protein function. We identified no haplotypes in complete association with the recorded phenotypes, supporting a complex inheritance of aggression.

A population study of a mutation allele associated with cone-rod dystrophy in the standard wire-haired dachshund.

Cone-rod dystrophy in the standard wire-haired dachshund (SWHD) is inherited as a simple autosomal recessive trait and the recently discovered mutation is widespread within the SWHD population in Norway and other Scandinavian countries. The gene frequency was estimated to be 4.8%. On the basis of the assumption that the size of the ancestral haplotype around a mutation is inversely correlated with the number of generations since the mutation arose, we have found that the mutation is of a relatively recent origin. The conserved haplotype was found to be 8 Mb in size and therefore we estimate that the mutation arose roughly eight generations (approximately 37 years) ago. This indicates that the mutation arose after breed separation.

A genome assembly-integrated dog 1 Mb BAC microarray: a cytogenetic resource for canine cancer studies and comparative genomic analysis.

Molecular cytogenetic studies have been instrumental in defining the nature of numerical and structural chromosome changes in human cancers, but their significance remains to be fully understood. The emergence of high quality genome assemblies for several model organisms provides exciting opportunities to develop novel genome-integrated molecular cytogenetic resources that now permit a comparative approach to evaluating the relevance of tumor-associated chromosome aberrations, both within and between species. We have used the dog genome sequence assembly to identify a framework panel of 2,097 bacterial artificial chromosome (BAC) clones, selected at intervals of approximately one megabase. Each clone has been evaluated by multicolor fluorescence in situ hybridization (FISH) to confirm its unique cytogenetic location in concordance with its reported position in the genome assembly, providing new information on the organization of the dog genome. This panel of BAC clones also represents a powerful cytogenetic resource with numerous potential applications. We have used the clone set to develop a genome-wide microarray for comparative genomic hybridization (aCGH) analysis, and demonstrate its application in detection of tumor-associated DNA copy number aberrations (CNAs) including single copy deletions and amplifications, regional aneuploidy and whole chromosome aneuploidy. We also show how individual clones selected from the BAC panel can be used as FISH probes in direct evaluation of tumor karyotypes, to verify and explore CNAs detected using aCGH analysis. This cytogenetically validated, genome integrated BAC clone panel has enormous potential for aiding gene discovery through a comparative approach to molecular oncology.

Evaluation of the serotonergic genes htr1A, htr1B, htr2A, and slc6A4 in aggressive behavior of golden retriever dogs.

Aggressive behavior displays a high heritability in our study group of Golden Retriever dogs. Alterations in brain serotonin metabolism have been described in aggressive dogs before. Here, we evaluate whether four genes of the canine serotonergic system, coding for the serotonin receptors 1A, 1B, and 2A, and the serotonin transporter, could play a major role in aggression in Golden Retrievers. We performed mutation screens, linkage analysis, an association study, and a quantitative genetic analysis. There was no systematic difference between the coding DNA sequence of the candidate genes in aggressive and non-aggressive Golden Retrievers. An affecteds-only parametric linkage analysis revealed no strong major locus effect on human-directed aggression related to the candidate genes. An analysis of 41 single nucleotide polymorphisms (SNPs) in the 1 Mb regions flanking the genes in 49 unrelated human-directed aggressive and 49 unrelated non-aggressive dogs did not show association of SNP alleles, genotypes, or haplotypes with aggression at the candidate loci. We completed our analyses with a study of the effect of variation in the candidate genes on a collection of aggression-related phenotypic measures. The effects of the candidate gene haplotypes were estimated using the Restricted Maximum Likelihood method, with the haplotypes included as fixed effects in a linear animal model. We observed no effect of the candidate gene haplotypes on a range of aggression-related phenotypes, thus extending our conclusions to several types of aggressive behavior. We conclude that it is unlikely that these genes play a major role in the variation in aggression in the Golden Retrievers that we studied. Smaller phenotypic effects of these loci could not be ruled out with our sample size.

Integrated cytogenetic BAC map of the genome of the gray, short-tailed opossum, Monodelphis domestica.

The generation of high-quality genome assemblies for numerous species is advancing at a rapid pace. As the number of genome assemblies increases, so does our ability to investigate genome relationships and their contributions to unraveling complex biological, evolutionary, and biomedical processes. A key process in the generation of a genome assembly is to determine and verify the precise physical location and order of the large sequence blocks (scaffolds) that result from the assembly. For organisms of relatively recent common ancestry this process may be achieved largely through comparative sequence alignment. However, as the evolutionary distance between species lengthens, the use of comparative sequence alignment becomes increasingly less reliable. Simultaneous cytogenetic mapping, using multicolor fluorescence in-situ hybridization (FISH) analysis, offers an alternative means to define the cytogenetic location and relative order of DNA sequences, thereby anchoring the genome sequence to the karyotype. In this article we report the molecular cytogenetic locations of 415 bacterial artificial chromosome (BAC) clones that served to anchor sequence scaffolds of the gray, short-tailed opossum (Monodelphis domestica) to its karyotype, which enabled accurate integration of these regions into the genome assembly.

A second-generation association study of the 5q31 cytokine gene cluster and the interleukin-4 receptor in asthma.

We have analyzed a dense set of single-nucleotide polymorphisms (SNPs) and microsatellites spanning the T-helper cytokine gene cluster (interleukins 3, 4, 5, 9, and 13, interferon regulatory factor-1, colony-stimulating factor-2, and T-cell transcription factor-7) on 5q31 and the gene encoding the interleukin-4 receptor (IL4R) on 16p12 among Finnish families with asthma. As shown by haplotype pattern mining analysis, the number of disease-associated haplotype patterns differed from that expected for the 129Q allele polymorphism in IL13 for high serum total immunoglobulin (Ig) E levels, but not for asthma. The same SNP also yielded the best haplotype associations. For IL4R, asthma-associated haplotype patterns, most spanning the S411L polymorphism, showed suggestive association. However, these haplotypes consisted of the major alleles for the intracellular part of the receptor and were very common among both patients and controls. The minor alleles 503P and 576R have been reported to be associated with decreased serum IgE levels and changes in the biological activity of the protein, especially when inherited together. In the Finnish population, these two polymorphisms segregated in strong linkage disequilibrium. Our data support previous findings regarding L4R, indicating that 503P and 576R may act as minor protecting alleles for IgE-mediated disorders.