Potential contribution of gut microbiota in the development of autoantibodies in T1D children carrying HLA-DRB1/DQB1 risk alleles: an experimental and in silico analysis.

This study aimed to investigate the prevalence of insulin autoantibody (IAA), glutamic acid decarboxylase antibody (GADA), and insulinoma-associated antigen-2 antibody (IA-2A) in type 1 diabetes (T1D) children based on the presence of predisposing HLA-II alleles. Additionally, to assess the sequence homology between autoantigens of islet cells and selected proteins derived from gut bacteria in terms of their binding capacities to HLA risk alleles, HLA-DRB1/DQB1 alleles were determined by PCR-SSOP in 111 T1D children (probands) along with 222 parents and 133 siblings. Autoantibodies were measured by ELISA, and in silico analysis was run as follows: protein extraction, homology and epitope prediction, peptide alignment, and HLA-peptide docking. Higher significant frequencies of DRB1*03:01, DQB1*02:01, and DQB1*03:02 alleles and DRB1*03:01 ~ DQB1*02:01 haplotype and lower frequencies of DRB1*11:01, DRB1*14:01, and DQB1*03:01 alleles were found in probands compared to parents and siblings. DRB1*11:01 ~ DQB1*03:01, DRB1*14:01 ~ DQB1*05:03, and DRB1*15:01-DQB1*06:02 haplotypes were significantly less frequent in the probands compared to parents. Out of 111 probands, 21 were seronegative, 90 tested positive for one autoantibody, and 15 showed the concurrent presence of three autoantibodies. Logistic regression analysis revealed that DRB1*04 ~ DQB1*03:02 haplotype was associated with the induction of GADA and IA-2A, while DRB1*11:01 ~ DQB1*03:01 was associated with seronegativity. Epitopes derived from GAD and gut bacteria showed strong binding capacities to HLA risk alleles. Due to the sequence similarities between gut bacteria-derived proteins and islet cell autoantigens and their potential for binding to HLA risk alleles, dysbiosis of gut microbiota can be considered another risk factor for the development of T1D, especially in genetically susceptible individuals.

Family-based association of HLA-DRB1 and DQB1 alleles and haplotypes in a group of Iranian Type 1 diabetes children.

This family-based study was conducted in a group of Iranians with Type 1 diabetes (T1D) to investigate the transmission from parents of risk and non-risk HLA alleles and haplotypes, and to estimate the genetic risk score for this disease within this population. A total of 240 T1D subjects including 111 parent-child trios (111 children with T1D, 133 siblings, and 222 parents) and 330 ethnically matched healthy individuals were recruited. High-resolution HLA typing for DRB1/DQB1 loci was performed for all study subjects (n = 925) using polymerase chain reaction-sequence-specific oligonucleotide probe method. The highest predisposing effect on developing T1D was conferred by the following haplotypes both in all subjects and in probands compared to controls: DRB1*04:05-DQB1*03:02 (Pc = 2.97e-06 and Pc = 6.04e-10, respectively), DRB1*04:02-DQB1*03:02 (Pc = 5.94e-17 and Pc = 3.86e-09, respectively), and DRB1*03:01-DQB1*02:01 (Pc = 8.26e-29 and Pc = 6.56e-16, respectively). Conversely, the major protective haplotypes included DRB1*13:01-DQB1*06:03 (Pc = 6.99e-08), DRB1*15:01-DQB1*06:02 (Pc = 2.97e-06) in the cases versus controls. Also, DRB1*03:01-DQB1*02:01/DRB1*04:02|05-DQB1*03:02 and DRB1*03:01-DQB1*02:01/DRB1*03:01-DQB1*02:01 diplotypes conferred the highest predisposing effect in the cases (Pc = 8.65e-17 and Pc = 6.26e-08, respectively) and in probands (Pc = 5.4e-15 and Pc = 0.001, respectively) compared to controls. Transmission disequilibrium test showed that the highest risk was conferred by DRB1*04:02-DQB1*03:02 (Pc = 3.26e-05) and DRB1*03:01-DQB1*02:01 (Pc = 1.78e-12) haplotypes and the highest protection by DRB1*14:01-DQB1*05:03 (Pc = 8.66e-05), DRB1*15:01-DQB1*06:02 (Pc = 0.002), and DRB1*11:01-DQB1*03:01 (Pc = 0.0003) haplotypes. Based on logistic regression analysis, carriage of risk haplotypes increased the risk of T1D development 24.5 times in the Iranian population (p = 5.61e-13). Also, receiver operating characteristic curve analysis revealed a high predictive power of those risk haplotypes in discrimination of susceptible from healthy individuals (area under curve: 0.88, p = 5.5e-32). Our study highlights the potential utility of genetic risk assessment based on HLA diplotypes for predicting T1D risk in individuals, particularly among family members of affected children in our population.

Determination of HLA class II risk alleles and prediction of self/non-self-epitopes contributing Hashimoto's thyroiditis in a group of Iranian patients.

One of the probable hypotheses for the onset of autoimmunity is molecular mimicry. This study aimed to determine the HLA-II risk alleles for developing Hashimoto's thyroiditis (HT) in order to analyze the molecular homology between candidate pathogen-derived epitopes and potentially self-antigens (thyroid peroxidase, TPO) based on the presence of HLA risk alleles. HLA-DRB1/-DQB1 genotyping was performed in 100 HT patients and 330 ethnically matched healthy controls to determine the predisposing/protective alleles for HT disease. Then, in silico analysis was conducted to examine the sequence homology between epitopes derived from autoantigens and four potentially relevant pathogens and their binding capacities to HLA risk alleles based on peptide docking analysis. We identified HLA-DRB1*03:01, *04:02, *04:05, and *11:04 as predisposing alleles and DRB1*13:01 as a potentially predictive allele for HT disease. Also, DRB1*11:04 ~ DQB1*03:01 (Pc = 0.002; OR, 3.97) and DRB1*03:01 ~ DQB1*02:01 (Pc = 0.004; OR, 2.24) haplotypes conferred a predisposing role for HT. Based on logistic regression analysis, carrying risk alleles increased the risk of HT development 4.5 times in our population (P = 7.09E-10). Also, ROC curve analysis revealed a high predictive power of those risk alleles for discrimination of the susceptible from healthy individuals (AUC, 0.70; P = 6.6E-10). Analysis of peptide sequence homology between epitopes of TPO and epitopes derived from four candidate microorganisms revealed a homology between envelop glycoprotein D of herpes virus and sequence 151-199 of TPO with remarkable binding capacity to HLA-DRB1*03:01 allele. Our findings indicate the increased risk of developing HT in those individuals carrying HLA risk alleles which can also be related to herpes virus infection.