Evidence of new risk genetic factor to systemic lupus erythematosus: the UBASH3A gene

The ubiquitin associated and Src-homology 3 (SH3) domain containing A (UBASH3a) is a suppressor of T-cell receptor signaling, underscoring antigen presentation to T-cells as a critical shared mechanism of diseases pathogenesis. The aim of the present study was to determine whether the UBASH3a gene influence the susceptibility to systemic lupus erythematosus (SLE) in Caucasian populations. We evaluated five UBASH3a polymorphisms (rs2277798, rs2277800, rs9976767, rs13048049 and rs17114930), using TaqMan® allelic discrimination assays, in a discovery cohort that included 906 SLE patients and 1165 healthy controls from Spain. The SNPs that exhibit statistical significance difference were evaluated in a German replication cohort of 360 SLE patients and 379 healthy controls. The case-control analysis in the Spanish population showed a significant association between the rs9976767 and SLE (Pc = 9.9E-03 OR = 1.21 95%CI = 1.07-1.37) and a trend of association for the rs2277798 analysis (P = 0.09 OR = 0.9 95%CI = 0.79-1.02). The replication in a German cohort and the meta-analysis confirmed that the rs9976767 (Pc = 0.02; Pc = 2.4E-04, for German cohort and meta-analysis, respectively) and rs2277798 (Pc = 0.013; Pc = 4.7E-03, for German cohort and meta-analysis, respectively) UBASH3a variants are susceptibility factors for SLE. Finally, a conditional regression analysis suggested that the most likely genetic variation responsible for the association was the rs9976767 polymorphism. Our results suggest that UBASH3a gene plays a role in the susceptibility to SLE. Moreover, our study indicates that UBASH3a can be considered as a common genetic factor in autoimmune diseases.

Islet-brain1/C-Jun N-terminal kinase interacting protein-1 (IB1/JIP-1) promoter variant is associated with Alzheimer's disease.

Islet-brain1 (IB1) or c-Jun NH2 terminal kinase interacting protein-1 (JIP-1), the product of the MAPK8IP1 gene, functions as a neuronal scaffold protein to allow signalling specificity. IB1/JIP-1 interacts with many cellular components including the reelin receptor ApoER2, the low-density lipoprotein receptor-related protein (LRP), kinesin and the Alzheimer's amyloid precursor protein. Coexpression of IB1/JIP-1 with other components of the c-Jun NH2 terminal-kinase (JNK) pathway activates the JNK activity; conversely, selective disruption of IB1/JIP-1 in mice reduces the stress-induced apoptosis of neuronal cells. We therefore hypothesized that IB1/JIP-1 is a risk factor for Alzheimer's disease (AD). By immunocytochemistry, we first colocalized the presence of IB1/JIP-1 with JNK and phosphorylated tau in neurofibrillary tangles. We next identified a -499A>G polymorphism in the 5' regulatory region of the MAPK8IP1 gene. In two separate French populations the -499A>G polymorphism of MAPK8IP1 was not associated with an increased risk to AD. However, when stratified on the +766C>T polymorphism of exon 3 of the LRP gene, the IB1/JIP-1 polymorphism was strongly associated with AD in subjects bearing the CC genotype in the LRP gene. The functional consequences of the -499A>G polymorphism of MAPK8IP1 was investigated in vitro. In neuronal cells, the G allele increased transcriptional activity and was associated with an enhanced binding activity. Taken together, these data indicate that the increased transcriptional activity in the presence of the G allele of MAPK8IP1 is a risk factor to the onset of in patients bearing the CC genotype of the LRP gene.

Examining ERBB2 as a candidate gene for susceptibility to leprosy (Hansen’s disease) in Brazil

Leprosy remains prevalent in Brazil. ErbB2 is a receptor for leprosy bacilli entering Schwann cells, which mediates Mycobacterium leprae-induced demyelination and the ERBB2 gene lies within a leprosy susceptibility locus on chromosome 17q11-q21. To determine whether polymorphisms at the ERBB2 locus contribute to this linkage peak, three haplotype tagging single nucleotide polymorphisms (tag-SNPs) (rs2517956, rs2952156, rs1058808) were genotyped in 72 families (208 cases; 372 individuals) from the state of Pará (PA). All three tag-SNPs were associated with leprosy per se [best SNP rs2517959 odds ratio (OR) = 2.22; 95% confidence interval (CI) 1.37-3.59; p = 0.001]. Lepromatous (LL) (OR = 3.25; 95% CI 1.37-7.70; p = 0.007) and tuberculoid (TT) (OR = 1.79; 95% CI 1.04-3.05; p = 0.034) leprosy both contributed to the association, which is consistent with the previous linkage to chromosome 17q11-q21 in the population from PA and supports the functional role of ErbB2 in disease pathogenesis. To attempt to replicate these findings, six SNPs (rs2517955, rs2517956, rs1810132, rs2952156, rs1801200, rs1058808) were genotyped in a population-based sample of 570 leprosy cases and 370 controls from the state of Rio Grande do Norte (RN) and the results were analysed using logistic regression analysis. However, none of the associations were replicated in the RN sample, whether analysed for leprosy per se, LL leprosy, TT leprosy, erythema nodosum leprosum or reversal reaction conditions. The role of polymorphisms at ERBB2 in controlling susceptibility to leprosy in Brazil therefore remains unclear.

LINGO1 and LINGO2 variants are associated with essential tremor and Parkinson disease.

Genetic variation in the leucine-rich repeat and Ig domain containing 1 gene (LINGO1) was recently associated with an increased risk of developing essential tremor (ET) and Parkinson disease (PD). Herein, we performed a comprehensive study of LINGO1 and its paralog LINGO2 in ET and PD by sequencing both genes in patients (ET, n=95; PD, n=96) and by examining haplotype-tagging single-nucleotide polymorphisms (tSNPs) in a multicenter North American series of patients (ET, n=1,247; PD, n= 633) and controls (n=642). The sequencing study identified six novel coding variants in LINGO1 (p.S4C, p.V107M, p.A277T, p.R423R, p.G537A, p.D610D) and three in LINGO2 (p.D135D, p.P217P, p.V565V), however segregation analysis did not support pathogenicity. The association study employed 16 tSNPs at the LINGO1 locus and 21 at the LINGO2 locus. One variant in LINGO1 (rs9652490) displayed evidence of an association with ET (odds ratio (OR) =0.63; P=0.026) and PD (OR=0.54; P=0.016). Additionally, four other tSNPs in LINGO1 and one in LINGO2 were associated with ET and one tSNP in LINGO2 associated with PD (P<0.05). Further analysis identified one tSNP in LINGO1 and two in LINGO2 which influenced age at onset of ET and two tSNPs in LINGO1 which altered age at onset of PD (P<0.05). Our results support a role for LINGO1 and LINGO2 in determining risk for and perhaps age at onset of ET and PD. Further studies are warranted to confirm these findings and to determine the pathogenic mechanisms involved.

A microRNA Signature Associated with Early Recurrence in Breast Cancer.

Recurrent breast cancer occurring after the initial treatment is associated with poor outcome. A bimodal relapse pattern after surgery for primary tumor has been described with peaks of early and late recurrence occurring at about 2 and 5 years, respectively. Although several clinical and pathological features have been used to discriminate between low- and high-risk patients, the identification of molecular biomarkers with prognostic value remains an unmet need in the current management of breast cancer. Using microarray-based technology, we have performed a microRNA expression analysis in 71 primary breast tumors from patients that either remained disease-free at 5 years post-surgery (group A) or developed early (group B) or late (group C) recurrence. Unsupervised hierarchical clustering of microRNA expression data segregated tumors in two groups, mainly corresponding to patients with early recurrence and those with no recurrence. Microarray data analysis and RT-qPCR validation led to the identification of a set of 5 microRNAs (the 5-miRNA signature) differentially expressed between these two groups: miR-149, miR-10a, miR-20b, miR-30a-3p and miR-342-5p. All five microRNAs were down-regulated in tumors from patients with early recurrence. We show here that the 5-miRNA signature defines a high-risk group of patients with shorter relapse-free survival and has predictive value to discriminate non-relapsing versus early-relapsing patients (AUC = 0.993, p-value<0.05). Network analysis based on miRNA-target interactions curated by public databases suggests that down-regulation of the 5-miRNA signature in the subset of early-relapsing tumors would result in an overall increased proliferative and angiogenic capacity. In summary, we have identified a set of recurrence-related microRNAs with potential prognostic value to identify patients who will likely develop metastasis early after primary breast surgery.

ESR1 gene promoter region methylation in free circulating DNA and its correlation with estrogen receptor protein expression in tumor tissue in breast cancer patients.

BACKGROUND Tumor expression of estrogen receptor (ER) is an important marker of prognosis, and is predictive of response to endocrine therapy in breast cancer. Several studies have observed that epigenetic events, such methylation of cytosines and deacetylation of histones, are involved in the complex mechanisms that regulate promoter transcription. However, the exact interplay of these factors in transcription activity is not well understood. In this study, we explored the relationship between ER expression status in tumor tissue samples and the methylation of the 5' CpG promoter region of the estrogen receptor gene (ESR1) isolated from free circulating DNA (fcDNA) in plasma samples from breast cancer patients. METHODS Patients (n = 110) with non-metastatic breast cancer had analyses performed of ER expression (luminal phenotype in tumor tissue, by immunohistochemistry method), and the ESR1-DNA methylation status (fcDNA in plasma, by quantitative methylation specific PCR technique). RESULTS Our results showed a significant association between presence of methylated ESR1 in patients with breast cancer and ER negative status in the tumor tissue (p = 0.0179). There was a trend towards a higher probability of ESR1-methylation in those phenotypes with poor prognosis i.e. 80% of triple negative patients, 60% of HER2 patients, compared to 28% and 5.9% of patients with better prognosis such as luminal A and luminal B, respectively. CONCLUSION Silencing, by methylation, of the promoter region of the ESR1 affects the expression of the estrogen receptor protein in tumors of breast cancer patients; high methylation of ESR1-DNA is associated with estrogen receptor negative status which, in turn, may be implicated in the patient's resistance to hormonal treatment in breast cancer. As such, epigenetic markers in plasma may be of interest as new targets for anticancer therapy, especially with respect to endocrine treatment.

Identification of a novel splice-site mutation in the CYP1A2 gene.

AIMS: To identify the molecular basis for a low CYP1A2 metabolic status, as determined by a caffeine phenotyping test, in a 71-year-old, nonsmoking, Caucasian woman who presented with very high clozapine concentrations despite being administered a standard dose of the drug. METHODS: The nucleotide sequence of the 7 exons, exon-intron boundaries and 5'-flanking region of the CYP1A2 gene was analysed by direct sequencing. RESULTS: Only one heterozygous point mutation was identified in the donor splice site of intron 6 (3534G > A) of CYP1A2. This mutation could cause abnormal RNA splicing and therefore lead to a truncated nonfunctional enzyme. No other carrier of this mutation was identified in a population of 100 unrelated healthy Caucasians. CONCLUSIONS: This is the first report of a splice-site mutation affecting the CYP1A2 gene. This polymorphism is a likely explanation for the low CYP1A2 activity associated with high clozapine concentrations in this patient.

Gene expression profiling in breast cancer: understanding the molecular basis of histologic grade to improve prognosis.

BACKGROUND: Histologic grade in breast cancer provides clinically important prognostic information. However, 30%-60% of tumors are classified as histologic grade 2. This grade is associated with an intermediate risk of recurrence and is thus not informative for clinical decision making. We examined whether histologic grade was associated with gene expression profiles of breast cancers and whether such profiles could be used to improve histologic grading. METHODS: We analyzed microarray data from 189 invasive breast carcinomas and from three published gene expression datasets from breast carcinomas. We identified differentially expressed genes in a training set of 64 estrogen receptor (ER)-positive tumor samples by comparing expression profiles between histologic grade 3 tumors and histologic grade 1 tumors and used the expression of these genes to define the gene expression grade index. Data from 597 independent tumors were used to evaluate the association between relapse-free survival and the gene expression grade index in a Kaplan-Meier analysis. All statistical tests were two-sided. RESULTS: We identified 97 genes in our training set that were associated with histologic grade; most of these genes were involved in cell cycle regulation and proliferation. In validation datasets, the gene expression grade index was strongly associated with histologic grade 1 and 3 status; however, among histologic grade 2 tumors, the index spanned the values for histologic grade 1-3 tumors. Among patients with histologic grade 2 tumors, a high gene expression grade index was associated with a higher risk of recurrence than a low gene expression grade index (hazard ratio = 3.61, 95% confidence interval = 2.25 to 5.78; P < .001, log-rank test). CONCLUSIONS: Gene expression grade index appeared to reclassify patients with histologic grade 2 tumors into two groups with high versus low risks of recurrence. This approach may improve the accuracy of tumor grading and thus its prognostic value.

The systemic lupus erythematosus IRF5 risk haplotype is associated with systemic sclerosis

Systemic sclerosis (SSc) is a fibrotic autoimmune disease in which the genetic component plays an important role. One of the strongest SSc association signals outside the human leukocyte antigen (HLA) region corresponds to interferon (IFN) regulatory factor 5 (IRF5), a major regulator of the type I IFN pathway. In this study we aimed to evaluate whether three different haplotypic blocks within this locus, which have been shown to alter the protein function influencing systemic lupus erythematosus (SLE) susceptibility, are involved in SSc susceptibility and clinical phenotypes. For that purpose, we genotyped one representative single-nucleotide polymorphism (SNP) of each block (rs10488631, rs2004640, and rs4728142) in a total of 3,361 SSc patients and 4,012 unaffected controls of Caucasian origin from Spain, Germany, The Netherlands, Italy and United Kingdom. A meta-analysis of the allele frequencies was performed to analyse the overall effect of these IRF5 genetic variants on SSc. Allelic combination and dependency tests were also carried out. The three SNPs showed strong associations with the global disease (rs4728142: P  = 1.34×10(-8), OR  = 1.22, CI 95%  = 1.14-1.30; rs2004640: P  = 4.60×10(-7), OR  = 0.84, CI 95%  = 0.78-0.90; rs10488631: P  = 7.53×10(-20), OR  = 1.63, CI 95%  = 1.47-1.81). However, the association of rs2004640 with SSc was not independent of rs4728142 (conditioned P  = 0.598). The haplotype containing the risk alleles (rs4728142*A-rs2004640*T-rs10488631*C: P  = 9.04×10(-22), OR  = 1.75, CI 95%  = 1.56-1.97) better explained the observed association (likelihood P-value  = 1.48×10(-4)), suggesting an additive effect of the three haplotypic blocks. No statistical significance was observed in the comparisons amongst SSc patients with and without the main clinical characteristics. Our data clearly indicate that the SLE risk haplotype also influences SSc predisposition, and that this association is not sub-phenotype-specific.

Infrequent promoter methylation of the MGMT gene in liver metastases from uveal melanoma.

Uveal melanoma is associated with a high mortality rate once metastases occur, with over >90% of metastatic patients dying within less than 1 year from metastases to the liver. The intraarterial hepatic (iah) administration of the alkylating agent fotemustine holds some promise with response rates of 36% and median survival of 15 months. Here, we investigated whether the DNA-repair-protein MGMT may be involved in the variability of response to fotemustine and temozolomide in uveal melanoma. Epigenetic inactivation of MGMT has been demonstrated to be a predictive marker for benefit from alkylating agent therapy in glioblastoma. We found a methylated MGMT promoter in 6% of liver metastases from 34 uveal melanoma patients. The mean MGMT activity measured in liver metastases with negligible liver tissue content was significantly lower than in liver tissue (146 versus 523 fmol/mg protein, p = 0.002). Expression of the MGMT protein was detectable in 50% of 88 metastases by immunohistochemistry on a tissue microarray. Expression was heterogeneous, and in accordance with MGMT activity data, usually lower than in the surrounding liver. Differential MGMT activity/expression between metastasis and liver tissue and more efficient depletion of MGMT with higher doses of alkylating agent therapy using iah delivery may provide the pharmacologic window for the higher response rate. However, these results do not support MGMT methylation status or protein expression as predictive markers for treatment outcome to iah chemotherapy with alkylating agents.

Gene-lifestyle interaction and type 2 diabetes: the EPIC interact case-cohort study.

BACKGROUND Understanding of the genetic basis of type 2 diabetes (T2D) has progressed rapidly, but the interactions between common genetic variants and lifestyle risk factors have not been systematically investigated in studies with adequate statistical power. Therefore, we aimed to quantify the combined effects of genetic and lifestyle factors on risk of T2D in order to inform strategies for prevention. METHODS AND FINDINGS The InterAct study includes 12,403 incident T2D cases and a representative sub-cohort of 16,154 individuals from a cohort of 340,234 European participants with 3.99 million person-years of follow-up. We studied the combined effects of an additive genetic T2D risk score and modifiable and non-modifiable risk factors using Prentice-weighted Cox regression and random effects meta-analysis methods. The effect of the genetic score was significantly greater in younger individuals (p for interaction  = 1.20×10-4). Relative genetic risk (per standard deviation [4.4 risk alleles]) was also larger in participants who were leaner, both in terms of body mass index (p for interaction  = 1.50×10-3) and waist circumference (p for interaction  = 7.49×10-9). Examination of absolute risks by strata showed the importance of obesity for T2D risk. The 10-y cumulative incidence of T2D rose from 0.25% to 0.89% across extreme quartiles of the genetic score in normal weight individuals, compared to 4.22% to 7.99% in obese individuals. We detected no significant interactions between the genetic score and sex, diabetes family history, physical activity, or dietary habits assessed by a Mediterranean diet score. CONCLUSIONS The relative effect of a T2D genetic risk score is greater in younger and leaner participants. However, this sub-group is at low absolute risk and would not be a logical target for preventive interventions. The high absolute risk associated with obesity at any level of genetic risk highlights the importance of universal rather than targeted approaches to lifestyle intervention.

Functional and molecular evolutionary analysis of the "Brevix Radix" gene family in mono-and dicotyledonous plants

ABSTRACTIn contrast to animals, plants cannot move from their place of birth and, therefore, need to adapt to their particular habitat in order to survive. Thus, plant development is remarkably plastic, making plants an ideal system for the isolation of genes that account for intraspecific natural variation and possibly environmental adaptation. However, to date, this approach mostly identified null alleles and missed mutations with subtle effects. For instance, BREVIS RADIX (BRX) has been isolated as a key regulator of root growth through a naturally occurring loss-of-function allele in the Arabidopsis thaliana accession Uk-1 and is the founding member of a highly-conserved plant-specific gene family.In this work, we show that a strong selective pressure is acting on the BRX gene family and dates back before the monocot-dicot divergence. However, functional diversification is observed mainly in dicotyledon BRX family genes and is correlated with acceleration in the evolutionary rates in the N-terminal regions. Population genetic data revealed that BRX is highly conserved across Arabidopsis accessions and presents signatures of adaptation. Interestingly, a seven amino acid deletion polymorphism in BRX sequence was found in a few accessions, which seems to be responsible for their enhanced primary root growth. Nevertheless, BRX might not only be active in the root, as suggested by its expression in the shoot. Indeed, leaves and cotyledons of brx mutants are significantly smaller than wild- type. This phenotype is a direct consequence of the absence of BRX function in the shoot rather than an indirect effect of an altered root system growth. Interestingly, cotyledons of brx plants reflect the same physiological defects as the root. Moreover, phenotypes in BRX gain-of-function plants, such as epinastic leaves and increased epidermal cell size, could be associated with an increase in leaf brassinosteroid content.Collectively, these results indicate that BRX contributes to local adaptation by ubiquitously regulating plant growth, probably through the modulation of brassinosteroid biosynthesis.RÉSUMÉContrairement à la plupart des animaux, les plantes ne peuvent se mouvoir et doivent ainsi s'adapter à leur environnement pour survivre. Pour cette raison, elles représentent un système idéal pour l'identification de gènes contribuant à la variation naturelle intra- spécifique, ainsi qu'à l'adaptation. Cependant, cette approche a, jusqu'à présent, surtout permis d'isoler des allèles nuls et non des mutations conférant des effets plus subtiles. C'est le cas du gène Β REVIS RADIX (BRX), un régulateur clé de la croissance racinaire, qui a été identifié grâce à un allèle non-fonctionnel présent dans l'accession naturelle d'Arabidopsis thaliana Uk-1. BRX et ses homologues des plantes mono- et dicotylédones forment une famille très conservée et spécifique aux plantes.Dans ce travail, nous démontrons que la famille de gènes BRX est soumise à une forte pression de sélection qui remonte avant la divergence entre mono- et dicotylédones. Cependant, une diversification fonctionnelle a été observée chez les gènes des dicotylédones et corrèle avec une accélération de la vitesse d'évolution dans leur région N- terminale. Une analyse génétique de différentes accessions naturelles d'Arabidopsis a révélé que BRX est hautement conservé et présente des signatures d'adaptation. Remarquablement, un polymorphisme de délétion de sept acides aminés a été détecté dans quelques accessions et a pour conséquence une plus forte croissance de la racine primaire. Néanmoins, il semble que le rôle de BRX ne se limite pas qu'à la racine, comme indiqué par son expression dans les parties aériennes de la plante. En effet, les mutants brx présentent des cotylédons et des feuilles significativement plus petits que le type sauvage, une conséquence directe de l'absence d'activité de BRX dans ces organes. Nous avons aussi noté que les cotylédons des mutants brx, à l'instar des racines, ont une perception altérée de l'auxine et peuvent être complémentés par l'application exogène de brassinostéroïdes. De plus, dans des plantes présentant un gain de fonction BRX, les feuilles sont épinastiques et les cellules de leur épiderme plus grandes. Ces phénotypes sont accompagnés d'une augmentation de la concentration de brassinostéroïdes dans les feuilles. Conjointement, ces résultats démontrent que BRX contribue à une adaptation locale de la plante par la régulation générale de sa croissance, probablement en modulant la biosynthèse des brassinostéroïdes.

Genomic acute myeloid leukemia-associated inv(16)(p13q22) breakpoints are tightly clustered.

The inv(16) and related t(16;16) are found in 10% of all cases with de novo acute myeloid leukemia. In these rearrangements the core binding factor beta (CBFB) gene on 16q22 is fused to the smooth muscle myosin heavy chain gene (MYH11) on 16p13. To gain insight into the mechanisms causing the inv(16) we have analysed 24 genomic CBFB-MYH11 breakpoints. All breakpoints in CBFB are located in a 15-Kb intron. More than 50% of the sequenced 6.2 Kb of this intron consists of human repetitive elements. Twenty-one of the 24 breakpoints in MYH11 are located in a 370-bp intron. The remaining three breakpoints in MYH11 are located more upstream. The localization of three breakpoints adjacent to a V(D)J recombinase signal sequence in MYH11 suggests a V(D)J recombinase-mediated rearrangement in these cases. V(D)J recombinase-associated characteristics (small nucleotide deletions and insertions of random nucleotides) were detected in six other cases. CBFB and MYH11 duplications were detected in four of six cases tested.

Blood pressure loci identified with a gene-centric array.

Raised blood pressure (BP) is a major risk factor for cardiovascular disease. Previous studies have identified 47 distinct genetic variants robustly associated with BP, but collectively these explain only a few percent of the heritability for BP phenotypes. To find additional BP loci, we used a bespoke gene-centric array to genotype an independent discovery sample of 25,118 individuals that combined hypertensive case-control and general population samples. We followed up four SNPs associated with BP at our p < 8.56 × 10(-7) study-specific significance threshold and six suggestively associated SNPs in a further 59,349 individuals. We identified and replicated a SNP at LSP1/TNNT3, a SNP at MTHFR-NPPB independent (r(2) = 0.33) of previous reports, and replicated SNPs at AGT and ATP2B1 reported previously. An analysis of combined discovery and follow-up data identified SNPs significantly associated with BP at p < 8.56 × 10(-7) at four further loci (NPR3, HFE, NOS3, and SOX6). The high number of discoveries made with modest genotyping effort can be attributed to using a large-scale yet targeted genotyping array and to the development of a weighting scheme that maximized power when meta-analyzing results from samples ascertained with extreme phenotypes, in combination with results from nonascertained or population samples. Chromatin immunoprecipitation and transcript expression data highlight potential gene regulatory mechanisms at the MTHFR and NOS3 loci. These results provide candidates for further study to help dissect mechanisms affecting BP and highlight the utility of studying SNPs and samples that are independent of those studied previously even when the sample size is smaller than that in previous studies.

Intronic variants in the NFKB1 gene may influence hearing forecast in patients with unilateral sensorineural hearing loss in Meniere's disease.

Meniere's disease is an episodic vestibular syndrome associated with sensorineural hearing loss (SNHL) and tinnitus. Patients with MD have an elevated prevalence of several autoimmune diseases (rheumatoid arthritis, systemic lupus erythematosus, ankylosing spondylitis and psoriasis), which suggests a shared autoimmune background. Functional variants of several genes involved in the NF-κB pathway, such as REL, TNFAIP3, NFKB1 and TNIP1, have been associated with two or more immune-mediated diseases and allelic variations in the TLR10 gene may influence bilateral affectation and clinical course in MD. We have genotyped 716 cases of MD and 1628 controls by using the ImmunoChip, a high-density genotyping array containing 186 autoimmune loci, to explore the association of immune system related-loci with sporadic MD. Although no single nucleotide polymorphism (SNP) reached a genome-wide significant association (p<10(-8)), we selected allelic variants in the NF-kB pathway for further analyses to evaluate the impact of these SNPs in the clinical outcome of MD in our cohort. None of the selected SNPs increased susceptibility for MD in patients with uni or bilateral SNHL. However, two potential regulatory variants in the NFKB1 gene (rs3774937 and rs4648011) were associated with a faster hearing loss progression in patients with unilateral SNHL. So, individuals with unilateral MD carrying the C allele in rs3774937 or G allele in rs4648011 had a shorter mean time to reach hearing stage 3 (>40 dB HL) (log-rank test, corrected p values were p = 0.009 for rs3774937 and p = 0.003 for rs4648011, respectively). No variants influenced hearing in bilateral MD. Our data support that the allelic variants rs3774937 and rs4648011 can modify hearing outcome in patients with MD and unilateral SNHL.