Analysis of the frequency of single nucleotide polymorphisms in cytokine genes in patients with New Onset Diabetes After Transplant

New Onset Diabetes After Transplantation (NODAT) is a serious metabolic complication. While β-cell dysfunction is considered the main contributing factor in the development of NODAT, the precise pathogenesis is not well understood. Cytokines are thought to be involved in the inflammation of islet β-cells in diabetes; however, few studies have investigated this hypothesis in NODAT. A total of 309 kidney transplant recipients (KTRs) were included in this study. An association between kidney transplants, and the development of diabetes after transplant (NODAT) was investigated. Comparison was made between KTRs who develop diabetes (NODAT cases) or did not develop diabetes (control), using key cytokines, IL-6 G (− 174)C, macrophage mediator; IL-4 C (− 490)T, T helper (Th)-2 cytokine profile initiator; Th-1 cytokine profile initiator interferon-γ T (+ 874) A gene and TGF β1 C (+ 869) T gene polymorphisms were investigated. The genes were amplified using well-established polymerase chain reaction (PCR) techniques in our laboratory. Compared to the AA and AT genotypes of interferon gamma (IFNG), there was a strong association between the TT genotype of IFNG and NODAT kidney transplant recipients (KTRs) versus non-NODAT KTRs (p = 0.005). The AA genotype of IFNG was found to be predominant in the control group (p = 0.004). Also, significant variations of IL6 G (− 174) C, IL-4 C (− 590) T, interferon-γ T (+ 874) A gene and transforming growth factor β1 C (+ 869) T may contribute to NODAT. Our data is consistent with theTh-1/T-reg pathway of immunity. Further larger pan Arab studies are required to confirm our findings.

www.nature.com/scientificreports/ scarce, with the exception of few epidemiologic studies [9][10][11] suggesting that 25-30% of KTRs eventually develop NODAT [9][10][11][12] . Specifically, the incidence rate of NODAT, following kidney transplantation, has been found in to be as high as 27% in Saudi Arabia 9 , 30% in Bahrain 10 , 22.2% in Egypt 11 , and 25.6% in Kuwait 12,13 . Genetic predisposition to NODAT likely involves an inherited defect in the peripheral tolerance to T-cell infiltration along with inadequate insulin secretion [14][15][16] . There is clear evidence showing that the different stages of NODAT are discrete entities, and that progression from one stage to the next is dependent on immune regulation 16 . T-cells are important for orchestrating the immune response and can be categorized into subsets according to their phenotypic characteristics resulting from polarization of naïve T-cells (i.e. Th-1, Th-2 and regulatory T-cell; T-reg) [17][18][19] . Each T-cell subset has a unique functional role highlighted by their capacity to produce pro-inflammatory and anti-inflammatory cytokines in response to various immune challenges. Cytokines are a group of pharmacologically active polypeptides that possess autocrine, paracrine, and juxtacrine effects with characteristic features 19 . Given the critical role of cytokines in regulating immune responses, subtle differences in cytokine expression may have a major effect on the outcomes of those responses 20 . Single nucleotide polymorphisms (SNP) in the regulatory region of the cytokine genes have been shown to correlate with individual variations in cytokine production [9][10][11][12][13][14][15][16][17][18][19][20][21][22][23] . Thus, variations in susceptibility to NODAT may be influenced by inter-individual genetic variations of cytokine genes [19][20][21][22][23] . Research on the association of SNPs have now become a potential means not only for better understanding the etiopathogenesis of the disease, but also as a probable marker for disease susceptibility and severity. It has been shown that interleukin-6(IL-6) derived from macrophages play a key role in inflammation 19,24 . This process is regulated via two major arms of the immune system: intra-cellular suppressor of cytokine signaling (SOCS) and T-reg cells 25 .
IL-6 gene (Online Mendelian Inheritance in Man #147620) maps to chromosome 7p21 26 . IL-6 a glycoprotein composed of 212 amino acids with a signal peptide of 27 amino acids and two potential NH2-linked glycosylation sites 26 . The molecular weight ranges from 21 to 28 kDa. IL-6 is a key cytokine with pleiotropic interactions in different human diseases pathogenesis 19 . IL-6 activates the immune system and enhances inflammatory response but also carries anti-inflammatory properties 27 . Understanding such pleiotropic effects of IL-6 in in the pathogenesis of the disease may help determine the progression, severity and duration of the disease. IL-6 is currently considered an important target for clinical interventions 28,29 . Hence, better understanding of IL-6 role in the pathogenesis of NODAT can provide rational and therapeutic intervention 29 . The guanine (G)/cytosine (C) polymorphism in the IL-6 proximal promoter region, at position IL-6 G (− 174) C, regulates transcription of the IL-6 gene and subsequently in IL-6 levels of plasma and serum 19,30 ; the same polymorphism was found to be associated with fasting insulin levels,insulin sensitivity, and the integrated area under the curve of serum glucose concentrations 30,31 .
Interferon-gamma (IFN-γ), also known as type II interferon or macrophage-activating factor 32 , is a multipotent cytokine with an approximate molecular weight of 17 kDa. IFN-γ is secreted by activated T-cells and natural killer cells (NK), and it modulates many facets of the hosts immune response 33 . It enhances the hosts defense and innate immune response through promoting inflammation via the JAK-STAT1 signaling pathway and inhibiting both interleukin (IL)-4 (the initiator of the Th-2 cytokine cascade) and the promoter of the Th-1 cytokine profile 22 . TheDNA sequence of the human IFN-γ gene (GenBank AF330164) shows the presence of a variable-length CA repeat in the first intron of the gene where allele 2 (12 CA repeats) produces a high level of IFN-γ protein 18 . This sequence is immediately adjacent and correlated with allele T for thiamine of the polymorphism at position + 874 34 . It is well established that the IFN-γ gene polymorphism T (+ 874) A (IFNG) of the first intron is correlated with serum level of IFN-γ production and mRNA expression in vitro and in vivo; the TT genotype correlates with high levels of IFN-γ production, and TA and AA are correlated with intermediate and low production level 18,22,[34][35][36][37][38] . This polymorphism coincides with a putative NF-κB binding site that may mediate high production of IFN-γ 39,40 . Importantly, IFN-γ is believed to play an important role in the autoimmune pathogenesis of type 1 diabetes 22,41 . Inhibition of the IFN-γ function in non-obese diabetic (NOD) mice using either IFN-γ-specific antibodies 42 or soluble IFN-γ receptors (IFN-γ-R) 43 reduced the incidence of spontaneous diabetes and also prevented the transfer of diabetes via splenocytes from NOD donor mice 44 . Furthermore, transgenic expression of IFN-γ by β-cells induced autoimmunity, resulting in overt diabetes in otherwise diabetes-resistant mice 44 . The resistance of IFN-γ-R deficient animals to CD4+ T-cell mediated diabetes was shown to be primarily due to a lack of IFN-γ-R expression by β-cells. Kuriya et al. demonstrated that IFN-γ-R deficiency has distinct effects on CD4+ compared with CD8+ T-cell mediated diabetes, and that IFN-γ may play a critical role in CD4+ T-cell mediated destruction of β-cells 45 .
IL-4 is a 20-KDa glycoprotein encoded by the IL-4 gene on chromosome 5q23.31. It is secreted by helper T-cells (CD4) type 2 (Th-2), by NK cells, and by cells of the innate immune system: mast cells, basophils, and eosinophils 46 . IL-4 regulates proliferation, apoptosis, gene expression, and differentiation in many hematopoietic cells; it directs the immunoglobulin (Ig0) class switch to IgG1 and IgE, downregulates the production of Th-1 cells and is a critical mediator of the Th 1/ Th-2 cytokine milieu counterbalance 22,46 . IL-4 appears to protect human islets from cytotoxic damage induced by proinflammatory and Th-1 cytokines. Another study showed that long-term exposure of rat pancreatic islets to IL-4 resulted in an inhibitory action to certain islet functions 47 . These phenomena occur by the modulation of the homing of autoreactive cells to inflammatory sites and the stabilization of a protective Th-2-mediated environment in the thymus, spleen, and pancreatic islets. Thus, IL-4 treatment favors the expansion of regulatory CD4+ Th-2 cells in vivo and prevents the onset of insulitis and type 1 diabetes mediated by autoreactive Th-1 cells 22,48 . It has been suggested that IL-4 protects human islets from cytotoxic damage induced by proinflammatory and Th-1 cytokines. The local expression of IL-4 in the pancreatic islets of NOD mice (ins-IL-4 mice) restricted the activation of autoreactive T-cells and promoted complete protection against spontaneous diabetes 49 . IL-4 T (− 590) C, rs2243250, is a functional promoter of gene polymorphism, where C is substituted by T 22  www.nature.com/scientificreports/ corresponds to high level of IL-4 serum and mRNA production in vivo as well as in vitro. The same is true for IL-4 C (− 590) C for low levels, and TC for putative intermediate protein levels 22,50 . Transforming growth factor-β1 (TGF-β1) belongs to a family of multifunctional growth factors which have profound regulatory effects on many developmental and physiological processes 51 .The human TGF-β1 gene is located on chromosome 19q13.1-13.3 52 , and more than ten polymorphic loci are presently known that are distributed across exons, introns, and the 5′flanking region 53 . A single nucleotide polymorphism (SNP) of codon 10 in the TGF-β1 T (+ 869) C gene(TFGB), T (proline) to C (leucine) is associated with different diseases including type 1 and type 2 diabetes 54 . There is well-established evidence regarding the association of codon 10 with varying levels of TGF-β1 synthesis in vitro and in vivo 23 . Increases or decreases in the production of TGF-β1 have been linked to numerous diseases including atherosclerosis, and fibrotic diseases of the kidney, liver, and lung 23 .
There is increasing evidence to show that SOCS proteins may be involved in the development of diabetes and its associated complications 55 . In addition, SOCS-modulating properties have been attributed to pharmacological agents that are currently used for the treatment of diabetes 56 .
Our results indicate that susceptibility to NODAT might be monitored by genotyping KTRs who developed diabetes compared with KTRs who did not develop diabetes for selected major Th 1(IFNG)/ Th-2 (IL-4)/TGFB, T-reg in addition to IL-6, macrophage derived cytokines. We have investigated well-established SNPs in major cytokines, according to their site of origin, in order to identify our future line of research.

Materials and methods
This piece of research was conducted following written approval from Dasman Diabetes Institute Research and Ethical Committee reference: RA 2015-013. Oral and written informed consent was obtained from all participants prior to data collection. The methods were performed in accordance with the relevant rules and regulations. This work was conducted in collaboration between Dasman Diabetes Institute and Hamed Al-Essa Organ Transplant Center of Kuwait. KTRs with NODAT and non-NODAT were recruited from May 2015 until March 2017. A total of 309 KTRs were recruited: 155 patients in the non-diabetic, non-NODAT cohort and 154 patients who developed NODAT following kidney transplantation (Table 1).
Genomic DNA was extracted from a 5 ml sample of fresh peripheral blood using the QIAamp DNA Mini Kit (QIAGEN) according to the manufacturer's instructions. Extracted DNA sample had a final concentration of 55-365 ng/ml, and samples were stored at −20°C prior to use. Allele-specific sequence primers were used to amplify IL-6G, IFNG, IL4G and TGFG using well-established techniques available in our laboratories ( Table 2) 22 , and PCR products were separated on a 2% agarose gel. To ensure the quality of the experiments, 10% of all samples were genotyped twice, producing 100% reproducibility.

Diagnosis of NODAT.
All kidney transplant recipients were screened for fasting plasma glucose and glycated hemoglobin, according to Kidney Disease: Improving Global Outcomes (KDIGO) guidelines 57 . Abnormal results were confirmed using an oral glucose tolerance test. In these cases, diabetes management was introduced immediately that included diet, exercise, oral agents and/or insulin, in addition to regular blood glucose monitoring at home. Immunosuppression protocol. The immunosuppression protocol consisted of five doses of anti-thymocyte globulin (Sanofi US, Bridgewater, NJ, USA) for high-risk patients (re-transplants, prior pregnancy, blood transfusion, HLA-antibody positive, and/or more than four HLA mismatches), or two doses of IL-2 receptor blocker (Basiliximab, Novartis, Inc., Switzerland) for low-risk patients. Maintenance therapy consisted of prednisolone, mycophenolic acid, and a calcineurin inhibitor (CNI). The dose of CNI was gradually decreased over 12 months guided by a 12-h trough level.
Acute cellular rejection was treated using intravenous methylprednisolone sodium succinate (solumedrol, 1 g daily for 3 days) or thymoglobulin (1 mg/kg for 7-10 days) for steroid-resistant rejection. Antibody-mediated rejection was treated using plasma exchange, intravenous immunoglobulin (2 g/kg), and rituximab. All rejection episodes were biopsy-proven according to the Banff criteria (2015) 58 . Patients who received thymoglobulin as anti-rejection therapy were managed by universal chemoprophylaxis for both cytomegalovirus (CMV) and Pneumocystis Jirovecii Pneumonia (PJP). Valganciclovir was used as CMV secondary prophylaxis for one month, while those who developed CMV viremia during this period were managed with a therapeutic dose for three weeks, followed by three months prophylaxis. Trimethoprim was used for one month as a prophylaxis for PJP. Associated infections were recorded if patients required hospital admission. Details of the patients who developed CMV infection or rejection episodes during the study period were recorded. Clinical data were collected with special emphasis on patient age, sex, donor type, immunosuppressive therapy, dialysis type and duration, primary kidney disease, pretransplant comorbidities, details of rejection episodes, post-transplant infections and graft and patient outcome.
Statistical analysis. The Hardy-Weinberg equilibrium was assessed using a chi-square test based on a comparison of the observed and expected genotypes. Statistical analyses were performed using SPSS software (SPSS, version 20.0, IBM Corporation, Armonk, NY, USA). The sample size was calculated to accept a marginal error of 6.6% (95% confidence interval) in a normally-distributed population. Allelic and genotypic distribution of IL6G, IFNG, IL-4G and TGFBG between both cohorts were compared using a paired-sample t-test, independent sample t-test, chi-square test, Fisher's exact test, and ANOVA, as appropriate. Results are expressed as mean ± standard deviation, and differences were considered significant at p < 0.05.  (Table 1). We noticed that the two groups were homogenous in variability of original kidney disease, dialysis type, donor type, and type of immunosuppression (both induction and maintenance). As well as that, all subjects had no history of diabetes as demonstrated by normal fasting and postprandial blood glucose levels prior to transplantation. There were no differences in ethnicity of our cohorts, 55% Kuwaiti, 45% non-Kuwaiti (Table 3). Moreover, no statistically significant differences were noticed in pre-transplant comorbidities in both cohorts, especially   www.nature.com/scientificreports/ hypertension, history tuberculosis treatment, ischemic heart disease, bone disease, anemia, and hyperlipidemia, HLA class 1 -A and -B and class II DR, p > 0.05, Table 3. There was a significantly higher prevalence of patients with chronic hepatitis C virus (HCV) infection in NODAT (12 cases in NODAT compared to two cases in control, p = 0.006, Table 3). In addition, there were significantly more CMV IgG positive patients for in NODAT than control (153 vs. 148, respectively; p = 0.02). However, the two groups were comparable regarding pre transplant CMV IgM (p > 0.05). Evaluation of immediate post-transplant graft function showed that the number of patients with slow graft function was significantly higher in control (p = 0.03). Overweight, as measured by mean basal body mass index (BMI), was significantly higher in the NODAT group compared to the control group (28.07 ± 5.5 vs. 26.11 ± 7, respectively; p = 0.01), while follow-up mean BMI was comparable in both cohorts irrespective of age, p = 0.21, see (Table 3). The number of cases with post-transplant CMV viremia was significantly higher in control compared to NODAT (p = 0.01, Table 4). Moreover, there was no significant difference between NODAT and control groups regarding both graft and patient outcome (p > 0.05).
Frequencies of cytokine genotypes according to protein production. Patients with of IL-6(CC), IFN-G (AA), and TGF-B (CC), corresponding to low protein production, were significantly higher in control (p = 0.05, p = 0.01, and p = 0.001, respectively, Fig. 1). IL-6(GG), IFN-G (TT), and TGF-B (TT) genotypes, correlating to high protein production, were significantly higher in NODAT (p = 0.05, p = 0.003 and p = 0.002, respectively). On the other hand, patients with IL-4 genotype (TT) that correspond to high IL-4 production were significantly higher in the control group while the low producers were significantly higher in NODAT, p = 0.001.
Age specificity and genotypic frequency. The allelic distribution of the IL-6, IL-4, IFNG and TGFB among both studied cohorts were significant, p = 0.006, p < 0.0001, p = 0.0002 and p < 0.0001 (Fig. 2) however, these effects were only observed when the alleles were paired in genotypes.
Genotype frequencies of IL-6, IL-4, IFNG and TGFB did not show any statistical gender specificity among NODAT versus control. IL-6 GG gene frequency was predominant in control participants < 40 years of age, p < 0.0001, while it was most frequent in patients 40-60 years of age in the NODAT group, p = 0.02. Although its distribution was greater in patients > 60 years of age in the NODAT group, this data did not rank to significance, p = 0.06 (Table 5). IL-6 GC was predominantly more frequent in control patients > 40 years of age versus NODAT, p=0.003 (  www.nature.com/scientificreports/ 0.002. The same was true for IL-6 CC genotype distribution. Its frequency was more in the controls among > 40 years of age patients, p = 0.0001, and significantly increased in > 60 years of age NODATs, p = 0.005, Table 5. IL-4 TT was more abundant in the controls among patients < 40 years of age when compared to NODAT, p = 0.01. However, TT distribution was comparable in 40-60 years of age groups. No statistically significant results were found between control and NODAT in the > 60 years of age group and this is likely due to the small number of cases, Table 4. Likewise, IL-4 TC was more frequent among older participants than 40 years in control, p < 0.0001, while its distribution was shifted in older age groups towards NODAT, p = 0.04, (among 40-60 years of age) and p < 0.0001 among > 60 years of age, Table 5. FNG AA, AT, TT was more dominant in controls less than 40 years, p = 0.0002, p < 0.0001, p < 0.0001 respectively. Its genotype frequencies were shifted towards NODAT in 40-60 years of age groups, p = 0.04, 0.01, 0.04 respectively (Table 5). In the same order, the genotype distributions of IFNG were predominantly more in NODAT than the control with > 60 years of age, p = 0.02, p = 0.02, p = 0.009 respectively ( Table 5).
The distribution of TGFB followed almost the same pattern as IFNG. The distribution of all genotypes were predominantly more in control among those older than 40 years, p < 0.0001, p = 0.001, p < 0.0001. The distributions were shifted towards dominance in NODAT among older age groups, p = 0.03, p = 0.4, p = 0.0007 in 40-60 Table 4. Virology compared with frequencies of IL-6, IL-4, IFNG and TGFB in NODAT versus control. IL-6 CG was more frequent in NODAT than control, p = 0.03; while IL-6 GG showed to be more frequent among control than NODAT, p = 0.003 in post-transplant CMV. IL-4 TC also was more frequent in NODAT than control, p = 0.006. On the other hand, IL-4 TC was more frequent on control than NODAT, p = 0.009. IFNG AA and AT were more frequent among NODAT compared to control in HCV positive individuals, p = 0.05, 0.007 respectively. While post-transplant CMV was not associated to the disease, TGFB did not show any statistically significant association with NODAT nor control in HCV positive individuals. While TGFB TC was significantly more in control subjects, p = 0.02. www.nature.com/scientificreports/ years of age and p = 0.002, p = 0.001and p = 0.03 among > 60 years of age groups ( Table 5). The frequency of the above cytokines was investigated according to their HLA class I -A and -B and class II DR but no statistical association was noted.

Discussion
The presence of NODAT has increasingly emerged as an important factor in the long-term outcome of KTRs, including mortality 59 . In the current study, the patient demographics of NODAT and control cohorts were indistinguishable. There were more males than females in both cohorts. This was similar to several other studies which reported that organ-recipients were mainly males. It is possible this reflects a gender bias in the incidence of transplant-related pathologies 60,61 . However, some differences were noted, patients with chronic HCV and CMV infections were significantly more prevalent amongst the NODAT group. Furthermore, in our study, age appears to be an important contributing factor in the development of NODAT. More KTRs < 40 years of age were found in the control group, while those > 40 years of age were more susceptible to NODAT. These findings are in line with several reports from Europe and the USA 1-3,62,63 . There was a 90% increase of relative risk in kidney transplant patients aged 45-59 years and a 160% increase in patients older than 60 years compared to patients between 18-44 years old 63 . However, our data did not support findings concluded from studies from Egypt, Bahrain, Saudi Arabia, and Iran [9][10][11]64 . This could be explained by transracial differences or a lack of large systematic studies. Associations of IL6 GG, IL-4 CC, IFNG TT and TGFB TT genotypes were superseded by age group and bio-physiology of patients. This might confirm identified differences in innate immune system in older KTRs compared to younger ones 65 .
Interestingly, patients with IL-6(GG), IFN-G (TT), and TGFB (TT) genotypes, which were associated with high protein production, were significantly higher in the NODAT cohort. While, IL-4(CC), related to low protein production, was significantly associated with the disease. On the contrary, HCV infection has been shown to be associated with IFNG and the development of type 2 diabetes mellitus in general population 66 .    67 , while other studies reported that CMV was not a risk factor for NODAT 62,63,68 . In pretransplant evaluation we found a significantly higher number of patients with positive CMV IgG in NODAT cohort, while CMV IgM was comparable between both cohorts. Interestingly, there were significantly more patients with post-transplantation CMV viremia in control cohort, especially among IL-6 (GG), IL-4 (TC) and TGFB (TC), but no statistical association was found between IFNG and CMV in either cohorts. This could be explained by the use of routine antiviral chemoprophylaxis among our patients after transplant which nullified CMV as a risk for NODAT in our cohort. The significant association of genetic polymorphisms in IL-6, IFNG, TGFB in NODAT could be due to Th1 mediated immunity which occurs as a result of deviation of SOCS/T-reg Fox P3 balance. Our results support a previously reported deviation in SOCS/T-reg Fox P3 balance causing immune dysregulation by IFN-γ. The association of IL-6 (GG), IFNG (TT) and TGFB (TT) high protein producers and IL4 (CC), low IL-4 protein producer might go some way to explain these findings.
KTRs gender did not appear to effect genetic susceptibility of either IL-6 G (− 147) C, IL-4 C (− 590) T, IFNG, and TGFB genotypes. However, there was a clear association with age and the distribution of IL-6 G (− 147) C, IL-4 C (− 590) T, IFNG, and TGFB genotypes among NODAT and control groups. The associations of IL6 G (− 147) G, IL-4 C (− 590) C, IFNG TT and TGFB TT were superseded by age group and bio-physiology of patients. In short, IL6 G (− 147) G, IL-4C (− 590) C, IFNG TT and TGFB TT were not independently associated to NODAT. The current data might indicate that the latter immune signature complex can provide a clue to the pathogenesis of NODAT. Although throughout our induction procedure, we aimed to suppress T-cells, thereby reducing the chance of kidney rejection.

Conclusion
The ability to perform noninvasive testing to transplant recipients will provide tools to identify patient risk of NODAT and the individualization of immune suppression regimens to improve outcomes after transplantation. The pathogenesis of NODAT is Th-1 cell-mediated variations, while IFNG, IL-4, TGF-β1 and IL-6may play a crucial role in that mechanism. Interestingly no association was found between the above cytokine genes and well-established HLA loci for NODAT. The latter rules out the autoimmune nature of NODAT.
Tailoring of immunosuppressive agents may be used to target patients according to their genetic makeup. The results of our research might provide a suitable platform for a larger multicenter study, focusing on the Arab population, to evaluate the role of cytokine genes in NODAT to confirm our findings and better understand the prevalence of NODAT in KTR patients. Table 5. Age at onset compared with frequencies of IL-6, IL-4, IFNG and TGFB in NODAT versus control. IL-6 GG was more frequent among controls than NODAT in < 40 years old individuals, p ≤ 0.0001; while, the frequency was more in NODAT than controls among 40-60 years old individuals, p = 0.03. IL-6 CG was more frequent among control than NODAT, p = 0.0001 in < 40 years individuals; while, it was significantly associated with the disease in > 60 years old, p = 0.002. In the same manner, IL-6 CC showed to more frequent in control in < 40 years old individuals, p = 0.0001; while the same genotype showed to be significantly associated to the disease in > 60 years old individuals, p = 0.005. IL-4 TT was more frequent in control than NODAT among < 40 years old individuals, p = 0.01; while it 699 was more in NODAT than control among > 60 years individuals, p = 0.03. In the same manner, IL-4 TC was significantly increased in control subjects in < 40 years old individuals, p ≤ 0.0001; while, its association was deviated towards NODAT in 40-60 and > 60 years old individuals, p = 0.04, < 0.0001 respectively. The same was true in IL-4 CC which was significantly increased in controls, p ≤ 0.0001; while its association was shifted towards NODAT in 40-60 years individuals, p = 0.006. In the same manner, IFNG AA also showed the same association with the disease in which it has been highly associated in controls in < 40 and its association was deviated towards NODAT in 40-60 and > 60 years old, p = 0.0002, 0.04 and 0.02 respectively. Interestingly, the same association was repeated among IFNG AT and TT. Where genotypes were highly associated with control subjects in < 40 and deviated then towards NODAT in older ages, p ≤ 0.0001, 0.01, 0.002, < 0.0001 and 0.03 respectively. Interestingly, the same manner has been followed up in TGFB CC, TC and TT genotypes. All were highly associated with controls in < 40 years old, p ≤ 0.0001, 0.0001 and < 0.0001 respectively. The genotypes association was deviated towards NODAT occurrence in 40-60 and > 60 years individuals, p = 0.03, 0.0002, 0.001, 0.0007 and 0.03 respectively. The above scenarios indicate lack of independent association of IFNG, IL-4, IL-6association of our key cytokines were not independent as they were following the old age association manner. www.nature.com/scientificreports/