HLA class II peptide-binding-region analysis reveals funneling of polymorphism in action.

BACKGROUND: HLA-class II proteins hold important roles in key physiological processes. The purpose of this study was to compile all class II alleles reported in human population and investigate patterns in pocket variants and their combinations, focusing on the peptide-binding region (PBR). METHODS: For this purpose, all protein sequences of DPA1, DQA1, DPB1, DQB1 and DRB1 were selected and filtered, in order to have full PBR sequences. Proportional representation was used for pocket variants while population data were also used. RESULTS: All pocket variants and PBR sequences were retrieved and analyzed based on the preference of amino acids and their properties in all pocket positions. The observed number of pocket variants combinations was much lower than the possible inferred, suggesting that PBR formation is under strict funneling. Also, although class II proteins are very polymorphic, in the majority of the reported alleles in all populations, a significantly less polymorphic pocket core was found. CONCLUSIONS: Pocket variability of five HLA class II proteins was studied revealing favorable properties of each protein. The actual PBR sequences of HLA class II proteins appear to be governed by restrictions that lead to the establishment of only a fraction of the possible combinations and the polymorphism recorded is the result of intense funneling based on function.

Identifying HLA DRB1-DQB1 alleles associated with Chlamydia trachomatis infection and in silico prediction of potentially-related peptides.

HLA class II (HLA-II) genes' polymorphism influences the immune response to Chlamydia trachomatis (Ct), it is considered a sexually transmitted infection. However, associations between HLA-II alleles and Ct-infection have been little explored in humans; this study was thus aimed at determining HLA-DRB1-DQB1 alleles/haplotypes' effect on Ct-infection outcome in a cohort of Colombian women. Cervical sample DNA was used as template for detecting Ct by PCR and typing HLA-DRB1-DQB1 alleles/haplotypes by Illumina MiSeq sequencing. Survival models were adjusted for identifying the alleles/haplotypes' effect on Ct-outcome; bioinformatics tools were used for predicting secreted bacterial protein T- and B-cell epitopes. Sixteen HLA-DRB1 alleles having a significant effect on Ct-outcome were identified in the 262 women analysed. DRB1*08:02:01G and DRB1*12:01:01G were related to infection-promoting events. Only the DQB1*05:03:01G allele related to clearance/persistence events was found for HLA-DQB1. HLA-DRB1 allele homozygous women were associated with events having a lower probability of clearance and/or early occurrence of persistence. Twenty-seven peptides predicted in silico were associated with protective immunity against Ct; outer membrane and polymorphic membrane protein-derived peptides had regions having dual potential for being T- or B-cell epitopes. This article describes HLA-DRB1-DQB1 alleles/haplotypes related to Ct-infection resolution and the peptides predicted in silico which might probably be involved in host immune response. The data provides base information for developing future studies leading to the development of effective prevention measures against Ct-infection.

NCAM protein and SARS-COV-2 surface proteins: In-silico hypothetical evidence for the immunopathogenesis of Guillain-Barré syndrome.

This study aimed at identifying human neural proteins that can be attacked by cross-reacting SARS-COV-2 antibodies causing Guillain-Barré syndrome. These markers can be used for the diagnosis of Guillain-Barré syndrome (GBS). To achieve this goal, proteins implicated in the development of GBS were retrieved from literature. These human proteins were compared to SARS-COV-2 surface proteins to identify homologous sequences using Blastp. Then, MHC-I and MHC-II epitopes were determined in the homologous sequences and used for further analysis. Similar human and SARS-COV-2 epitopes were docked to the corresponding MHC molecule to compare the binding pattern of human and SARS-COV-2 proteins to the MHC molecule. Neural cell adhesion molecule is the only neural protein that showed homologous sequence to SARS-COV-2 envelope protein. The homologous sequence was part of HLA-A68 and HLA-DQA/HLA-DQB epitopes had a similar binding pattern to SARS-COV-2 envelope protein. Based on these results, the study suggests that NCAM may play a significant role in the immunopathogenesis of GBS. NCAM antibodies can be used as a marker for Guillain-Barré syndrome. However, more experimental studies are needed to prove these results.

Solid-State NMR Investigations of the MHC II Transmembrane Domains: Topological Equilibria and Lipid Interactions.

The major histocompatibility complex class II (MHC II) membrane proteins are key players in the adaptive immune response. An aberrant function of these molecules is associated with a large number of autoimmune diseases such as diabetes type I and chronic inflammatory diseases. The MHC class II is assembled from DQ alpha 1 and DQ beta 1 which come together as a heterodimer through GXXXG-mediated protein-protein interactions and a highly specific protein-sphingomyelin-C18 interaction motif located on DQA1. This association can have important consequences in regulating the function of these membrane proteins. Here, we investigated the structure and topology of the DQA1 and DQB1 transmembrane helical domains by CD-, oriented 2H and 15N solid-state NMR spectroscopies. The spectra at peptide-to-lipid ratios of 0.5 to 2 mol% are indicative of a topological equilibrium involving a helix crossing the membrane with a tilt angle of about 20° and another transmembrane topology with around 30° tilt. The latter is probably representing a dimer. Furthermore, at the lowest peptide-to-lipid ratio, a third polypeptide population becomes obvious. Interestingly, the DQB1 and to a lesser extent the DQA1 transmembrane helical domains exhibit a strong fatty acyl chain disordering effect on the inner segments of the 2H-labelled palmitoyl chain of POPC bilayers. This phosphatidylcholine disordering requires the presence of sphingomyelin-C18 suggesting that the ensemble of transmembrane polypeptide and sphingolipid exerts positive curvature strain.

Critical amino acid variations in HLA-DQB1* molecules confers susceptibility to Autoimmune Thyroid Disease in south India.

The HLA-DQB1* region exhibits complex associations with autoimmune thyroid disease (AITD). AITD patients (Hashimoto's thyroiditis, HT = 180; Graves' disease, GD = 55) and age/sex matched controls (n = 235) were genotyped for DQB1* alleles by PCR-SSP. Alleles DQB1*02:02, *06:03, *06:09, *03:02, and *03:03 showed an increased risk and *02:01, *05:02, and *06:02 showed a protection toward AITD. Multiple sequence alignment was used to find out the amino acid variations within the peptide-binding pockets of susceptible and/or protective DQB1* alleles. We observed susceptible associations for amino acids 'Glu86(P < 0.0007)' and 'Leu87(P < 3.8 × 10-4)' in P1, 'Leu26(P < 4.0 × 10-12)' in P4, 'His9(P < 5.0 × 10-4)' and 'Ala57(P < 3.6 × 10-4)' in P9 toward HT; and 'Gly86(P < 0.0004)' in P1 and 'Asp57(P < 1.9 × 10-4)' in P9 towards GD. Protective associations were observed for amino acids 'Ala86(P < 8.2 × 10-6)' and 'Tyr87(P < 0.0003)' in P1, 'Gly26(P < 4.9 × 10-5)' and 'Ser74(P < 4.9 × 10-5)' in P4, 'Phe9(P < 0.0007)' and 'Ser57(P < 0.0016)' in P9 towards HT. Thus, the present study revealed that DQB1* alleles and putative amino acid residues play an important role in susceptibility toward AITD in south India.

Combined effect of glutamine at position 70 of HLA-DRB1 and alanine at position 57 of HLA-DQB1 in type 1 diabetes: An epitope analysis.

The contribution of specific HLA Class II alleles in type 1 diabetes is determined by polymorphic amino acid epitopes that direct antigen binding therefore, along with conventional allele frequency analysis, epitope analysis can provide important insights into disease susceptibility. We analyzed the highly heterogeneous Cypriot population for the HLA class II loci of T1DM patients and controls and we report for the first time their allele frequencies. Within our patient cohort we identified a subgroup that did not carry the DRB1*03:01-DQA1*05:01-DQB1*02:01 and DRB1*04:xx-DQA1*03:01-DQB1*03:02 risk haplotypes but a novel recombinant one, DRB1*04:XX-DQA1*03:01-DQB1*02:01 designated DR4-DQ2.3. Through epitope analysis we identified established susceptibility (DQB1 A57, DRB1 H13) and resistance (DQB1 D57) residues as well as other novel susceptibility residues DRB1 Q70, DQB1 L26 and resistance residues DRB1 D70, R70 and DQB1 Y47. Prevalence of susceptibility epitopes was higher in patients and was not exclusively a result of linkage disequilibrium. Residues DRB1 Q70, DQB1 L26 and A57 and a 10 amino acid epitope of DQA1 were the most significant in discriminating risk alleles. An extended haplotype containing these epitopes was carried by 92% of our patient cohort. Sharing of susceptibility epitopes could also explain the absence of risk haplotypes in patients. Finally, many significantly associated epitopes were non-pocket residues suggesting that critical immune functions may exist spanning further from the binding pockets.

Unraveling the structural basis for the unusually rich association of human leukocyte antigen DQ2.5 with class-II-associated invariant chain peptides.

Human leukocyte antigen (HLA)-DQ2.5 (DQA1*05/DQB1*02) is a class-II major histocompatibility complex protein associated with both type 1 diabetes and celiac disease. One unusual feature of DQ2.5 is its high class-II-associated invariant chain peptide (CLIP) content. Moreover, HLA-DQ2.5 preferentially binds the non-canonical CLIP2 over the canonical CLIP1. To better understand the structural basis of HLA-DQ2.5's unusual CLIP association characteristics, better insight into the HLA-DQ2.5·CLIP complex structures is required. To this end, we determined the X-ray crystal structure of the HLA-DQ2.5· CLIP1 and HLA-DQ2.5·CLIP2 complexes at 2.73 and 2.20 Å, respectively. We found that HLA-DQ2.5 has an unusually large P4 pocket and a positively charged peptide-binding groove that together promote preferential binding of CLIP2 over CLIP1. An α9-α22-α24-α31-ß86-ß90 hydrogen bond network located at the bottom of the peptide-binding groove, spanning from the P1 to P4 pockets, renders the residues in this region relatively immobile. This hydrogen bond network, along with a deletion mutation at α53, may lead to HLA-DM insensitivity in HLA-DQ2.5. A molecular dynamics simulation experiment reported here and recent biochemical studies by others support this hypothesis. The diminished HLA-DM sensitivity is the likely reason for the CLIP-rich phenotype of HLA-DQ2.5.

Genetic Variation at Exon 2 of the MHC Class II DQB Locus in Blue Whale (Balaenoptera musculus) from the Gulf of California.

The genes of the Major Histocompatibility Complex (MHC) play an important role in the vertebrate immune response and are among the most polymorphic genes known in vertebrates. In some marine mammals, MHC genes have been shown to be characterized by low levels of polymorphism compared to terrestrial taxa; this reduction in variation is often explained as a result of lower pathogen pressures in marine habitats. To determine if this same reduction in variation applies to the migratory population of blue whales (Balaenoptera musculus) that occurs in the Gulf of California, we genotyped a 172 bp fragment of exon 2 of the MHC Class II DQB locus for 80 members of this population. Twenty-two putatively functional DQB allotypes were identified, all of which were homologous with DQB sequences from other cetacean species. Up to 5 putative alleles per individual were identified, suggesting that gene duplication has occurred at this locus. Rates of non-synonymous to synonymous substitutions (ω) and maximum likelihood analyses of models of nucleotide variation provided potential evidence of ongoing positive selection at this exon. Phylogenetic analyses of DQB alleles from B. musculus and 16 other species of cetaceans revealed trans-specific conservation of MHC variants, suggesting that selection has acted on this locus over prolonged periods of time. Collectively our findings reveal that immunogenic variation in blue whales is comparable to that in terrestrial mammals, thereby providing no evidence that marine taxa are subject to reduced pathogen-induced selective pressures.

Questionable expression of unstable DQ heterodimer containing HLA-DQA1*01:07.

Human leukocyte antigens (HLA)-DQA1*01:07 was identified as an HLA-DQ blank specificity that segregated with the serological HLA-A2, -B7, -DR14, -DR52 haplotype, which carried DQB1*05:03. The blank specificity of DQA1*01:07-DQB1*05:03 may be because of lack of reactivity of available typing sera, or disruption of proper assembly of DQ heterodimer. The cDNA sequence of DQA1*01:07 is nearly identical to DQA1*01:04 except for a variant at position 304, which results in the replacement of an arginine with a cysteine at 79α. To determine whether the DQA1*01:07 product can be expressed on cell-surface, we co-expressed DQA1*01:07 with various DQB1*05 or *06 alleles in fibroblast cells. Cell-surface expression of DQ was detectable when DQA1*01:07 was co-expressed with DQB1*06:04 but undetectable with other DQB1*05 and DQB1*06 alleles, including DQB1*05:03, to which DQA1*01:07 was encoded in cis. These data suggest that DQA1*01:07 may act as a phenotypically null allele in the DQA1*01:07-DQB1*05:03 haplotype, while it can be expressed at a low level in the presences of certain DQB1*06 alleles, such as DQB1*06:04, in trans. Based on the null or low expression of DQA1*01:07 as shown in the previous and present studies, DQA1*01:07 has recently been renamed to DQA1*01:07Q, indicating its questionable expression.

Allele and Haplotype Frequencies of Human Leukocyte Antigen-A, -B, -C, -DRB1, and -DQB1 From Sequence-Based DNA Typing Data in Koreans.

BACKGROUND: Data on allele frequencies (AFs) and haplotype frequencies (HFs) of HLA-C and -DQB1 are limited in Koreans. We investigated AFs and HFs of HLA-A, -B, -C, -DRB1, and -DQB1 in Koreans by high-resolution sequence-based typing (SBT). METHODS: Hematopoietic stem cells were obtained from 613 healthy, unrelated donors to analyze HLA-A, -B, -C, -DRB1, and -DQB1 genotypes by using AlleleSEQR HLA-A, -B, -C, -DRB1, and -DQB1 SBT kits (Abbott Molecular, USA), respectively. Alleles belonging to HLA-C*07:01/07:06 group were further discriminated by using PCR-sequence specific primer analysis. AFs and HFs were calculated by direct counting and maximum likelihood method, respectively. RESULTS: In all, 24 HLA-A, 46 HLA-B, 24 HLA-C, 29 HLA-DRB1, and 15 HLA-DQB1 alleles were identified. AFs and HFs of HLA-A, -B, and -DRB1 were similar to those reported previously. For the HLA-C locus, C*01:02 was the most common allele, followed by C*03:03, C*03:04, C*14:02, C*03:02, and C*07:02 (AF ≥7%). AFs of C*07:01 and C*07:06 were 0.16% and 3.18%, respectively. For the HLA-DQB1 locus, DQB1*03:01 was the most common allele, followed by DQB1*03:03, *03:02, *06:01, *05:01, *04:01, and *06:02 (AF ≥7%). AFs of DQB1*02:01 and DQB1*02:02 were 2.12% and 6.69%, respectively. HFs of A*33:03-C*07:06 and C*07:06-B*44:03 were 3.09% and 3.10%, respectively, while those of DRB1*07:01-DQB1*02:02 and DRB1*03:01-DQB1*02:01 were 6.61% and 2.04%, respectively. CONCLUSIONS: This study reported AFs and HFs of HLA, including HLA-C and -DQB1, in Koreans by using high-resolution SBT. These data can be used to resolve ambiguous results of HLA typing for organ and hematopoietic stem cell transplantations.

HLA dosage effect in narcolepsy with cataplexy.

Narcolepsy with cataplexy is a sleep disorder caused by the loss of hypocretin-producing neurons in the hypothalamus. It is tightly associated with a specific human leukocyte antigen (HLA)-allele: HLA-DQB1*06:02. Based on this, an autoimmune process has been hypothesized. A functional HLA-DQ molecule consists of a DQα and a DQß chain. HLA-DQB1*06:02 (DQß) has a strong preference for binding to HLA-DQA1*01:02 (DQα), and together they form the functional DQ0602 dimer. A dosage effect would be expected if the HLA-DQ0602 dimer itself is directly involved in the aetiology. An increased expression of the HLA-DQ0602 dimer is expected in individuals homozygous for HLA-DQB1*06:02-DQA1*01:02, but is also hypothesized in individuals heterozygous for HLA-DQB1*06:02 and homozygous for HLA-DQA1*01:02. To study the impact of the expression of the HLA-DQ0602 dimer on narcolepsy susceptibility, 248 Dutch narcolepsy patients and 1272 Dutch control subjects, all of them positive for DQB1*06:02 (heterozygous and homozygous), were HLA-genotyped with attention not only to DQB1 but also to DQA1*01:02. DQB1*06:02-DQA1*01:02 homozygosity was significantly more often seen in patients compared to controls (O.R. 2.29) confirming previous observations. More importantly, a significantly higher prevalence of homozygosity for DQA1*01:02 was found in HLA-DQB1*06:02 heterozygous patients compared to controls (O.R. 2.37, p < 0.001). The latter finding clearly supports a direct role of the HLA-DQ molecule in the development of disease.

Antigenic peptide molecular recognition by the DRB1-DQB1 haplotype modulates multiple sclerosis susceptibility.

Multiple sclerosis (MS) is an autoimmune disease of the central nervous system that has a notably high incidence in Sardinia. Our study focuses on two HLA class II haplotypes associated with the disease in Sardinia, the rare predisposing DRB1*15:01-DQB1*06:02 and the widespread protective DRB1*16:01-DQB1*05:02. This framework enabled the highlighting of HLA binding pocket specificity and peptide recognition mechanisms by employing molecular dynamics simulations of the whole DRB1-DQB1 haplotype interacting with MBP- and EBV-derived peptides. We analyzed peptide-protein interaction networks and temporal evolution of the original complexes and after key amino acid mutations. The mutation G86V of the protective DRB1 allele exerted its effect mainly in the presence of the EBV viral peptide, with local and long range outcomes. However, the V38A mutation of the protective DQB1 showed a long range effect only in the case of the MBP myelin peptide. Our findings also demonstrate a DRB1/DQB1 complementary molecular recognition of peptides. This mechanism could provide a robust synergistic action and a differential role of DRB1 and DQB1 in tissues and in the time-steps towards autoimmunity. In addition, we demonstrate that negatively charged residues in pockets 4 and 9 play a role in MS susceptibility. Our findings are supported by recent experiments using a closely related MS animal model. Overall, our analysis confirms the role of the DRB1-DQB1 haplotype in conferring disease predisposition and could provide a valuable aid in designing optimal therapeutic peptides for MS therapy.

Interaction between HLA-DRB1-DQB1 haplotypes in Sardinian multiple sclerosis population.

We performed a case-control study in 2,555 multiple sclerosis (MS) Sardinian patients and 1,365 healthy ethnically matched controls, analyzing the interactions between HLA-DRB1-DQB1 haplotypes and defining a rank of genotypes conferring a variable degree of risk to the disease. Four haplotypes were found to confer susceptibility (*13:03-*03:01 OR = 3.3, Pc 5.1 × 10(-5), *04:05-*03:01 OR = 2.1, Pc 9.7 × 10(-8), *15:01-*06:02 OR = 2.0, Pc = 9.1 × 10(-3), *03:01-*02:01 OR = 1.7 Pc = 7.9 × 10(-22)) and protection (*11, OR = 0.8, Pc = 2.7 × 10(-2), *16:01-*05:02 OR = 0.6, Pc = 4.8 × 10(-16), *14:01-4-*05:031 = OR = 0.5, Pc = 9.8 × 10(-4) and *15:02-*06:01 OR = 0.4, Pc = 5.1 × 10(-4)). The relative predispositional effect method confirms all the positively associated haplotypes and showed that also *08 and *04 haplotypes confers susceptibility, while the *11 was excluded as protective haplotype. Genotypic ORs highlighted two typologies of interaction between haplotypes: i) a neutral interaction, in which the global risk is coherent with the sum of the single haplotype risks; ii) a negative interaction, in which the genotypic OR observed is lower than the sum of the OR of the two haplotypes. The phylogenic tree of the MS-associated DRB1 alleles found in Sardinian patients revealed a cluster represented by *14:01, *04:05, *13∶03, *08:01 and *03:01 alleles. Sequence alignment analysis showed that amino acids near pocket P4 and pocket P9 differentiated protective from predisposing alleles under investigation. Furthermore, molecular dynamics simulation performed on alleles revealed that position 70 is crucial in binding of MBP 85-99 peptide. All together, these data suggest that propensity to MS observed in Sardinian population carried by the various HLA-DRB1-DQB1 molecules can be due to functional peculiarity in the antigen presentation mechanisms.

MHC class II association study in eight breeds of dog with hypoadrenocorticism.

Canine hypoadrenocorticism is an endocrine disorder characterised by inadequate secretion of steroid hormones from the adrenal glands. Pathology results from immune-mediated destruction of the adrenal cortex, which is similar to that seen in the human Addison's disease. Both the canine and human diseases have similar clinical presentation, with the diagnosis based on performing a dynamic adrenocorticotropic hormone stimulation test. MHC class II has previously been associated with the human and canine diseases. In the current study, we conducted an MHC class II association study in eight breeds of dog with diagnoses of hypoadrenocorticism. We demonstrated significant differences in dog leukocyte antigen (DLA) haplotype frequencies in six of these breeds: Cocker spaniel, Springer spaniel, Labrador, West Highland white terrier (WHWT), Bearded collie, and Standard poodle. In the Springer spaniel, the DLA-DRB1*015:01--DQA1*006:01--DQB1*023:01 haplotype was significantly associated with disease risk (p = 0.014, odds ratio (OR) = 5.14) and showed a similar trend in the Cocker spaniel. This haplotype is related to one associated with hypoadrenocorticism in the Nova Scotia duck tolling retriever. Similar haplotypes shared between breeds were demonstrated, with DLA-DRB1*001:01--DQA1*001:01--DQB1*002:01 more prevalent in both affected Labrador (p = 0.0002, OR = 3.06) and WHWT (p = 0.01, OR = 2.11). Other haplotypes that have not previously been associated with the disease were identified. The inter-breed differences in DLA haplotypes associated with susceptibility to canine hypoadrenocorticism could represent divergent aetiologies. This could have implications for clinical diagnosis and future comparative studies. Alternatively, it may suggest that the gene of interest is closely linked to the MHC.

Structural and functional studies of trans-encoded HLA-DQ2.3 (DQA1*03:01/DQB1*02:01) protein molecule.

MHC class II molecules are composed of one α-chain and one ß-chain whose membrane distal interface forms the peptide binding groove. Most of the existing knowledge on MHC class II molecules comes from the cis-encoded variants where the α- and ß-chain are encoded on the same chromosome. However, trans-encoded class II MHC molecules, where the α- and ß-chain are encoded on opposite chromosomes, can also be expressed. We have studied the trans-encoded class II HLA molecule DQ2.3 (DQA1*03:01/DQB1*02:01) that has received particular attention as it may explain the increased risk of certain individuals to type 1 diabetes. We report the x-ray crystal structure of this HLA molecule complexed with a gluten epitope at 3.05 Å resolution. The gluten epitope, which is the only known HLA-DQ2.3-restricted epitope, is preferentially recognized in the context of the DQ2.3 molecule by T-cell clones of a DQ8/DQ2.5 heterozygous celiac disease patient. This preferential recognition can be explained by improved HLA binding as the epitope combines the peptide-binding motif of DQ2.5 (negative charge at P4) and DQ8 (negative charge at P1). The analysis of the structure of DQ2.3 together with all other available DQ crystal structures and sequences led us to categorize DQA1 and DQB1 genes into two groups where any α-chain and ß-chain belonging to the same group are expected to form a stable heterodimer.

Gluten-specific T cells cross-react between HLA-DQ8 and the HLA-DQ2α/DQ8ß transdimer.

Because susceptibility to celiac disease is associated strongly with HLA-DQ2 (DQA1*05/DQB1*02) and weakly with HLA-DQ8 (DQA1*03/DQB1*03), a subset of patients carries both HLA-DQ2 and HLA-DQ8. As a result, these patients may express two types of mixed HLA-DQ2/8 transdimers (encoded by DQA1*05/DQB1*03 and DQA1*03/DQB1*02) in addition to HLA-DQ2 and HLA-DQ8. Using T cells from a celiac disease patient expressing HLA-DQ8trans (encoded by DQA*0501/DQB*0302), but neither HLA-DQ2 nor HLA-DQ8, we demonstrate that this transdimer is expressed on the cell surface and can present multiple gluten peptides to T cell clones isolated from the duodenum of this patient. Furthermore, T cell clones derived from this patient and HLA-DQ2/8 heterozygous celiac disease patients respond to gluten peptides presented by HLA-DQ8trans, as well as HLA-DQ8, in a similar fashion. Finally, one gluten peptide is recognized better when presented by HLA-DQ8trans, which correlates with preferential binding of this peptide to HLA-DQ8trans. These results implicate HLA-DQ8trans in celiac disease pathogenesis and demonstrate extensive T cell cross-reactivity between HLA-DQ8 and HLA-DQ8trans. Because type 1 diabetes is strongly associated with the presence of HLA-DQ8trans, our findings may bear relevance to this disease as well.