Genetic factors are indicated in the development of type 1 diabetes (DM1). Recently, nucleotide variants of BACH2 and SOD2 have been associated with this chronic condition. Therefore, the purpose of the present study was to investigate the contribution of BACH2 rs3757247 and SOD2 rs4880 (Ala16Val) polymorphisms to the risk of DM1 and diabetes long-term complications. Selected polymorphic variants of BACH2 and SOD2 were investigated in a group of 141 patients with DM1 and in a group of age, gender-matched healthy subjects (n=369) using a high-resolution melting curve method. There was no evidence for either allelic or genotypic association with the risk of DM1 and diabetes chronic complications for analysed polymorphisms. In addition, no interaction between BACH2 and SOD2 variants in the development of this condition was observed. However, the frequency of BACH2 rs3757247 AG and AA genotypes was statistically different between DM1 patients with retinopathy and healthy individuals (odds ratio, 2.455; 95% confidence interval, 0.999–6.035; P=0.044), but this result did not survive multiple testing corrections. The present study did not confirm the involvement of BACH2 rs3757247 and SOD2 rs4880 polymorphisms in the development of DM1 and diabetes long-term complications. Further studies in a larger population sample are required.
type 1 diabetesBACH2 rs3757247 polymorphismSOD2 rs4880 polymorphismIntroduction
Type 1 diabetes (DM1) is a multifactorial disease that results from autoimmune destruction of β-cells in the pancreas (1). DM1 characterises the largest range of concordance rates in monozygotic twins among all autoimmune diseases, and this provides convincing evidence that genetic factors are strongly involved in the pathogenesis of DM1 (2). Initially, based on family-based linkage analyses, the strong association between human leukocyte antigen (HLA) class II genes encoding cell-surface antigen-presenting proteins and DM1 was found (3). With the advent of high-throughput single-nucleotide polymorphism (SNP) genotyping array technology studies, searching for the novel DM1 loci became possible (4). Based on the results of these studies, >60 new loci were identified across the human genome, which were associated with the risk of DM1 (5). In addition, advanced genetic technologies have become an important strategy to identify genetic factors contributing to the development of chronic diabetic complications (6,7).
Recently, follow-up analysis of a genome-wide association study was performed on DM1 patients and revealed that an SNP in the intron of the BTB and CNC homology 1, basic leucine zipper transcription factor 1 (BACH2) gene (rs3757247) was associated with the risk of DM1 (8).
The protein product of BACH2 plays a role as a key regulator of nucleic acid-triggered antiviral responses in human cells (9). It has been shown that Bach2 is expressed in the early stages of B-cell differentiation and is suppressed in terminally differentiated plasma cells by repressing the expression of B lymphocyte-induced maturation protein-1, the key regulator of plasma cell differentiation. The observed expression profile may serve to sensitize cells to oxidative stress and/or allow for the effective elimination of cells at early stages of differentiation (10). BACH2 is also critical for class switch recombination and somatic hyper mutations (SHM) of immunoglobulin genes. Analysis of transgenic Bach2-/- mice revealed that BACH2 is a key regulator of antibody response, which is required for T-cell independent and dependent immunoglobulin G (IgG) responses and SHM (11). Additionally, BACH2 may contribute to the development of DM1 by regulation of proinflammatory cytokine-induced apoptotic pathways in pancreatic β-cells by crosstalk with the PTPN2 gene and activation of JNK1 and BIM (12).
The missense T to C nucleotide gene variant of SOD2 (rs4880) is another SNP associated with the risk of DM1 and its chronic complications (13). The gene product of SOD2 is the free radical scavenging mitochondrial SOD2 that catalyses the dismutation of superoxide into oxygen and hydrogen peroxide (14). The presence of the C allele at the −9 position of the SOD2 gene results in a substitution of alanine for valine (Ala16Val) in the mitochondrial targeting sequence of the enzyme. This substitution is associated with less efficient transport of SOD2 into the mitochondrial matrix and decreases the efficiency of SOD2 to neutralise superoxide radical (15). The production of superoxide in DM1 is enhanced by frequent episodes of hyperglycaemia (16). This results in the dysfunction and apoptosis of endothelial cells, which contribute to the development of chronic complications (17). Therefore, the SOD2 CC variant may increase the risk of DM1 and its chronic complications, which was shown in certain clinical studies (18,19).
Based on the results of the studies showing the significance of the association between the BACH2 gene rs3757247 and SOD2 gene rs4880 polymorphisms and risk of DM1 and diabetes chronic complication pathogenesis (8,15), the analysis of the distribution of selected polymorphic variants of BACH2 and SOD2 in a Polish DM1 population was performed.
Materials and methodsStudy group
The study group was composed of 141 DM1 patients from the Department of Internal Medicine and Diabetology, Poznan University of Medical Sciences (Poznan, Poland). DM1 in all the patients was diagnosed according to the criteria from the American Diabetes Association on the basis of classical symptoms, blood glucose concentration >11.1 mmol/l and C-peptide concentration <0.5 µg/l (20). The mean age of DM1 patients was 33±10 years and the mean diabetes duration was 15±8 years. All the patients underwent a screening for diagnosis of microvascular complications, such as diabetic retinopathy, diabetic kidney disease and diabetic neuropathy.
Diabetic retinopathy was diagnosed by direct ophthalmoscopy through dilated pupils followed by fundus photography in all the patients. Diabetic retinopathy was graded according to the classification of the American Academy of Ophthalmology into classifications of: No retinopathy; mild, nonproliferative retinopathy (mild, moderate and severe); and proliferative retinopathy (21).
Diabetic kidney disease was detected at the stage of albuminuria (urinary albumin excretion rate between 30 and 300 mg/24 h) in two samples collected over a 3-month period after exclusion of secondary causes of microproteinuria. Diabetic kidney disease was defined as the presence of albuminuria in connection with diabetes of >10 years in duration or with diagnosed diabetic retinopathy (22).
Diabetic neuropathy was diagnosed in patients with two or more of the following components: The presence of neuropathy symptoms, the absence of ankle tendons on reflexes and abnormal scores for pressure and/or vibration perception (20).
The control group consisted of 369 healthy, age- and gender-matched individuals randomly selected from blood donors and healthy volunteers. All the study participants were Caucasians of Polish origin. The study received approval from the Ethics Committee of Poznan University of Medical Sciences (no. 607/12) and written informed consent was obtained from all the participants.
Biochemical parameters
Analysis of biochemical parameters was performed in an accredited medical laboratory. The serum concentration of total cholesterol (TCH) was determined using the commercially available assay kits (Roche Diagnostics GmbH, Mannheim, Germany). The glycated hemoglobin A1c (HbA1c) value was measured using high-performance liquid chromatography with the Variant Hemoglobin A1c Program (Bio-Rad, Hercules, CA, USA) (23).
SNP selection and genotyping
BACH2 rs3757247 and SOD2 rs4880 polymorphic variants were selected based on their association with DM1 and diabetes complications in previous studies. Characteristics of those SNPs are presented in Table I. DNA was isolated from peripheral blood lymphocytes by salt extraction. Genotyping was carried out by high-resolution melting (HRM) curve analysis on the LightCycler® 480 system (Roche Diagnostics GmbH) with the use of the 5x HOT FIREPol® EvaGreen® HRM mix (Solis BioDyne, Tartu, Estonia). For quality control, ∼10% of the randomly chosen samples were re-genotyped. Samples that failed genotyping were not repeated and were excluded from statistical analyses. Primer sequences and conditions for HRM analysis are presented in Table II.
Statistical analysis
For each SNP, the Hardy-Weinberg equilibrium was assessed by the χ2 test in patients and controls and P<0.05 was considered to indicate a statistically significant difference. The differences in allele and genotype frequencies between cases and controls were determined using standard χ2 or Fisher tests. SNPs were tested for the association with DM1 or chronic diabetic complications using the Cochran-Armitage trend test. The strength of the association was measured by odds ratios (ORs) with 95% confidence intervals (CIs). The data were analysed under recessive and dominant inheritance models. A statistical adjustment for multiple comparisons was accomplished using the Bonferroni correction.
Statistical significance was interpreted as P<0.05. The gene-gene interaction among tested SNPs was analysed using the logistic regression and epistasis option in the PLINK software (http://pngu.mgh.harvard.edu/purcell/plink/). PLINK generates a model based on allele dosage for each SNP and considers allelic by allelic epistasis. Values of P<0.05 were considered to indicate a statistically significant difference. To assess which parameter independently influences the presence of chronic complications in DM1 patients, multivariable regression analysis was performed as follows: Age of DM1 verses HbA1c level verses DM1 duration verses variants of BACH2 rs3757247 and SOD2 rs4880 polymorphisms.
ResultsPatient characteristics
The basic characteristics of the patients are presented in Table III. There were 81 females and 60 males among the Polish patients with DM1. The mean of TCH concentration was within the reference level [4.0 (3.9–4.4) mmol/l], whereas the mean level of HbA1c [8.0 (7.1–8.9) %] showed that the majority of the patients did not maintain good glycaemic control.
Retinopathy was diagnosed in 38 patients (27%), diabetic kidney disease in 8 patients (6%) and neuropathy in 14 patients (10%). Of these patients, 25 had retinopathy only (18%), 4 patients had only neuropathy (3%) and 16 patients (11%) had more than one microvascular complication.
Associations of the tested polymorphisms
None of the tested polymorphisms showed evidence of deviation from Hardy-Weinberg equilibrium in either DM1 patients or controls (P>0.05). Genotype counts, genotypic OR and corresponding 95% CI calculations for BACH2 rs3757247 and SOD2 rs4880 are reported in Table IV. There was no evidence for either allelic or genotypic association with the risk of DM1 for the analysed nucleotide variants. The minor allele frequency was 0.48 in cases and 0.44 in controls for BACH2 rs3757247 (pallelic=0.327), and 0.48 (allele T) in cases and 0.48 (allele C) in controls for SOD2 rs4880 (pallelic=0.332). The calculated OR for individuals with the combined AA and AG genotypes of the BACH2 rs3757247 polymorphism compared to the GG homozygotes was 1.554 (95% CI, 0.989–2.441). However, this result was not statistically significant (P=0.055, pcorr=0.110).
The gene-gene interaction analysis showed no significant epistatic interaction between the BACH2 and SOD2 variants in the development of DM1 (Table V). The OR for interaction was 1.066 (P=0.763).
No significant association between the BACH2 rs3757247 or SOD2 rs4880 genotypes with different chronic diabetic complications was identified (Table VI). However, there was an exception of the borderline association with diabetic retinopathy for the BACH2 rs3757247 variant. Compared to the GG homozygotes, individuals with the AG or AA genotypes (dominant model of inheritance) had an OR of 2.455 (95% CI, 0.999–6.035; P=0.044). However, the significance of this result did not survive multiple testing corrections (pcorr=0.088). The result of multivariate regression analysis has shown that only the duration of DM1 will independently influence the presence of chronic complications in DM1 patients (β=0.473, R2=0.300, P<0.05).
Discussion
Previous results of large meta-analysis studies show that the BACH2 rs3757247 and SOD2 rs4800 polymorphisms played a role in the development of DM1 (8,13). The A allele of rs3757247 polymorphism in BACH2 increased the risk of DM1 (8), whereas the C allele of the rs4880 polymorphism in the SOD2 gene had protective effects on the risk of DMI (13). In the present study, this was inconclusive. Most likely, one of the reasons of a negative association between the rs37577247 polymorphism in BACH2 and DM1 risk is the extremely modest contribution of non-HLA loci to the development of DM1 associated with the HLA region (24). In addition, it was shown that the rs3757247 polymorphism of BACH2 facilitated autoimmunity rather than increased the risk of DM1 itself (25,26). In the genome-wide association analysis of autoantibody positivity in DM1 patients, the association between the minor allele of the BACH2 rs3757247 variant and higher positivity of thyroid peroxidase autoantibodies was identified (25). This led to the conclusion that this polymorphic variant increased the risk of Graves' disease (26), not DM1 itself. Additionally, it was shown that Bach2-/- mice present deficient T-cell-independent and dependent IgG responses, which shows that Bach2 is important for the immune reaction and played an important role in autoimmune response (11).
Notably, the rs3757247 polymorphism of BACH2 may be associated with a risk of retinopathy. This finding may arise from the role of BACH2 in the induction of apoptosis under oxidative stress. BACH2 is a proapoptotic factor as it inhibits anti-oxidative and anti-apoptotic genes in the response of enhanced production of reactive oxygen species (ROS) (27). Oxidative stress is associated with the diabetic condition and is one of the most common mechanisms leading to the development of retinopathy (28), therefore we believe that the rs3757247A allele of BACH2 could be involved in enhanced apoptosis of retinal cells. However, the significance of the present result did not survive multiple testing corrections. Consequently, an additional large-scale study should be performed to conclude whether the rs3757247 polymorphism of BACH2 represents novel loci in the pathogenesis of DM1 and/or diabetes retinopathy.
Thus far, there is only one study showing the association of the rs4880 polymorphism of SOD2 with the risk of DM1 (12). The majority of the studies indicate that this polymorphism impacted the development of microvascular complications in patients with diabetes (29–31). The Val/Val genotype of rs4880 was shown to increase the risk of retinopathy in DM1 patients; however, no association was found for nephropathy (29). Notably, studies performed on Finish and Swedish populations of DM1 patients indicated that the homozygosity for the SOD2 Val allele is associated with an increased risk of nephropathy (30). Similar results were shown for Japanese type 2 diabetic patients (31). In the present study, no associations of the SOD2 rs4880 polymorphism with any of the diagnosed microvascular complications in the studied group were identified. These observations are consistent with the findings of el-Masry et al (19), who demonstrated that the Val/Val genotype of rs4880 is not a significant factor in diabetic patients with nephropathy. In addition, no association was identified between rs4880 polymorphisms and the risk of microvascular complications in DM1 or DM2 patients. Notably, the association between macroangiopathy and the lower frequency of the C (Ala) allele of SOD2 rs4880 was identified in this study (32). The Val isoform of SOD2 may lead to a decreased resistance against ROS produced in the mitochondria and to oxidative damage of proteins caused by less efficient SOD2 transport into the mitochondria (33). Therefore, we expected to find the association between the TT genotype of the SOD2 gene and the presence of long-term complications, caused partially by enhanced oxidative stress (28). One of the explanations for the negative result could be due to a relatively small sample size of the DM1 group. It has to be considered that due to this limitation, the statistical analysis may have failed to demonstrate any significant differences. Another reason of failure to replicate an association of SNPs for candidate genes could be too low a number of selected SNPs to achieve full gene coverage.
By contrast, the result of linear regression showed that the DM1 duration independently influenced the development of microvascular complications in DM1 patients, indicating the conventional risk factors as: Diabetes duration and high glycosylated level play a primarily role in the development of microangiopathy in DM1 (34).
In conclusion, neither BACH2 rs3757247 nor SOD2 rs4880 increase the risk of DM1 or chronic diabetic complications in Polish DM1 patients. To confirm this result, further investigations on larger groups should be performed.
Acknowledgements
The study was supported by a grant for young scientists from the University of Medical Sciences in Poznan (grant no. 502-14-02 228 370-09 003).
ReferencesWatkinsRAEvans-MolinaCBlumJSDiMeglioLAEstablished and emerging biomarkers for the prediction of type 1 diabetes: A systematic reviewTransl Res164110121201410.1016/j.trsl.2014.02.00424662515StankovKBencDDraskovicDGenetic and epigenetic factors in etiology of diabetes mellitus type 1Pediatrics13211121122201310.1542/peds.2013-165224190679SteckAKRewersMJGenetics of type 1 diabetesClin Chem57176185201110.1373/clinchem.2010.14822121205883Wellcome Trust Case Control ConsortiumGenome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controlsNature447661678200710.1038/nature0591117554300BakayMPandeyRHakonarsonHGenes involved in type 1 diabetes: An updateGenes Basel4499521201310.3390/genes403049924705215SandholmNSalemRMMcKnightAJNew susceptibility loci associated with kidney disease in type 1 diabetesPLoS Genet8e1002921201210.1371/journal.pgen.100292123028342HuangYCLinJMLinHJChenCCChenSYTsaiCHTsaiFJGenome-wide association study of diabetic retinopathy in a Taiwanese populationOphthalmology118642648201110.1016/j.ophtha.2010.07.02021310492GrantSFQuHQBradfieldJPFollow-up analysis of genome-wide association data identifies novel loci for type 1 diabetesDiabetes58290295200910.2337/db08-102218840781OyakeTItohKMotohashiHHayashiNHoshinoHNishizawaMYamamotoMIgarashiKBach proteins belong to a novel family of BTB-basic leucine zipper transcription factors that interact with MafK and regulate transcription through the NF-E2 siteMol Cell Biol166083609519968887638MutoATashiroSTsuchiyaHKumeAKannoMItoEYamamotoMIgarashiKActivation of Maf/AP-1 repressor Bach2 by oxidative stress promotes apoptosis and its interaction with promyelocytic leukemia nuclear bodiesJ Biol Chem2772072420733200210.1074/jbc.M11200320011923289MutoATashiroSNakajimaOHoshinoHTakahashiSSakodaEIkebeDYamamotoMIgarashiKThe transcriptional programme of antibody class switching involves the repressor Bach2Nature429566571200410.1038/nature0259615152264MarroquíLSantinIDos SantosRSMarselliLMarchettiPEizirikDLBACH2, a candidate risk gene for type 1 diabetes, regulates apoptosis in pancreatic β-cells via JNK1 modulation and crosstalk with the candidate gene PTPN2Diabetes6325162527201410.2337/db13-144324608439TianCFangSDuXJiaCAssociation of the C47T polymorphism in SOD2 with diabetes mellitus and diabetic microvascular complications: A meta-analysisDiabetologia54803811201110.1007/s00125-010-2004-521181397ChristiansonDWStructural chemistry and biology of manganese metalloenzymesProg Biophys Mol Biol67217252199710.1016/S0079-6107(97)88477-59446936SuttonAKhouryHPrip-BuusCCepanecCPessayreDDegoulFThe Ala16Val genetic dimorphism modulates the import of human manganese superoxide dismutase into rat liver mitochondriaPharmacogenetics13145157200310.1097/00008571-200303000-0000412618592FiorentinoTVPriolettaAZuoPFolliFHyperglycemia-induced oxidative stress and its role in diabetes mellitus related cardiovascular diseasesCurr Pharm Des1956955703201310.2174/138161281131932000523448484de M BandeiraSda FonsecaLJda S GuedesGRabeloLAGoulartMOVasconcelosSMOxidative stress as an underlying contributor in the development of chronic complications in diabetes mellitusInt J Mol Sci1432653284201310.3390/ijms1402326523385234ChistyakovDASavostanovKVZotovaEVNosikovVVPolymorphisms in the Mn-SOD and EC-SOD genes and their relationship to diabetic neuropathy in type 1 diabetes mellitusBMC Med Genet24200110.1186/1471-2350-2-411299047el-MasryTMZahraMAel-TawilMMKhalifaRAManganese superoxide dismutase alanine to valine polymorphism and risk of neuropathy and nephropathy in Egyptian type 1 diabetic patientsRev Diabet Stud27074200510.1900/RDS.2005.2.7017491681American Diabetes AssociationStandards of medical care in diabetes - 2007Diabetes Care30(Suppl 1)S4S41200710.2337/dc07-S00417192377WilkinsonCPFerrisFLIIIKleinRELeePPAgardhCDDavisMDillsDKampikAPararajasegaramRVerdaguerJTProposed international clinical diabetic retinopathy and diabetic macular edema disease severity scalesOphthalmology11016771682200310.1016/S0161-6420(03)00475-513129861KDOQIKDOQI clinical practice guidelines and clinical practice recommendations for diabetes and chronic kidney diseaseAm J Kidney Dis49(Suppl 2)S12S154200710.1053/j.ajkd.2006.12.00517276798JeppssonJOKoboldUBarrJApproved IFCC reference method for the measurement of HbA1c in human bloodClin Chem Lab Med407889200210.1515/CCLM.2002.01611916276SteckAKWongRWagnerBJohnsonKLiuERomanosJWijmengaCNorrisJMEisenbarthGSRewersMJEffects of non-HLA gene polymorphisms on development of islet autoimmunity and type 1 diabetes in a population with high-risk HLA-DR, DQ genotypesDiabetes61753758201210.2337/db11-122822315323de GraaffLCMartín-MartorellPBaanJBallieuxBSmitJWRadderJKLong-term follow-up of organ-specific antibodies and related organ dysfunction in type 1 diabetes mellitusNeth J Med696671201121411842PlagnolVHowsonJMSmythDJGenome-wide association analysis of autoantibody positivity in type 1 diabetes casesPLoS Genet7e1002216201110.1371/journal.pgen.100221621829393TashiroSMutoATanimotoKTsuchiyaHSuzukiHHoshinoHYoshidaMWalterJIgarashiKRepression of PML nuclear body-associated transcription by oxidative stress-activated Bach2Mol Cell Biol2434733484200410.1128/MCB.24.8.3473-3484.200415060166TarrJMKaulKChopraMKohnerEMChibberRPathophysiology of diabetic retinopathyISRN Ophthalmol2013343560201310.1155/2013/34356024563789HovnikTDolzanVBratinaNUPodkrajsekKTBattelinoTGenetic polymorphisms in genes encoding antioxidant enzymes are associated with diabetic retinopathy in type 1 diabetesDiabetes Care3222582262200910.2337/dc09-085219752172MöllstenAMarklundSLWessmanMSvenssonMForsblomCParkkonenMBrismarKGroopPHDahlquistGA functional polymorphism in the manganese superoxide dismutase gene and diabetic nephropathyDiabetes56265269200710.2337/db06-069817192491NomiyamaTTanakaYPiaoLNagasakaKSakaiKOgiharaTNakajimaKWatadaHKawamoriRThe polymorphism of manganese superoxide dismutase is associated with diabetic nephropathy in Japanese type 2 diabetic patientsJ Hum Genet48138141200312624725FlekacM1SkrhaJHilgertovaJLacinovaZJarolimkovaMGene polymorphisms of superoxide dismutases and catalase in diabetes mellitusBMC Med Genet930200810.1186/1471-2350-9-3018423055WangLIMillerDPSaiYLiuGSuLWainJCLynchTJChristianiDCManganese superoxide dismutase alanine-to-valine polymorphism at codon 16 and lung cancer riskJ Natl Cancer Inst9318181821200110.1093/jnci/93.23.181811734599PedroRARamonSAMarcBBJuanFBIsabelMMPrevalence and relationship between diabetic retinopathy and nephropathy, and its risk factors in the North-East of Spain, a population-based studyOphthalmic Epidemiol17251265201010.3109/09286586.2010.49866120642348
Characteristics of the polymorphisms genotyped in the SOD2 and BACH2 genes.
Gene symbol
Chromosome location
rs number
SNP function
Allelesa
MAFb
BACH2
6q15
rs3757247
intron
A/G (REV)
0.47
SOD2
6q25.3
rs4880
missense (p.Val16Ala)
C/T (REV)
0.46
According to the Single Nucleotide Polymorphism database (dbSNP); underlining denotes the minor allele.