- Open Access
The combined role of allelic variants of IRS-1 and IRS-2 genes in susceptibility to type2 diabetes in the Punjabi Pakistani subjects
Diabetology & Metabolic Syndrome volume 11, Article number: 64 (2019)
Diabetes mellitus is a multifactorial disorder characterized by a high level of glucose in the blood. Both genetic and environmental factors interact to cause diabetes. Insulin receptor substrate (IRS) proteins have a significant part in insulin signaling pathways. We aimed to investigate the relationship of type 2 diabetes with a Gly972Arg (G972R) variant of the IRS-1 gene and Gly1057Asp (G1057D) polymorphism of IRS-2 gene in the population of Punjab, Pakistan.
We collected 926 samples, 500 healthy controls (fasting blood sugar < 99 mg/dL, random blood sugar < 126 mg/dL) and 426 cases with diabetes (fasting blood sugar > 99 mg/dL, random blood sugar > 126 mg/dL). Several anthropometric measurements were measured. Statistical analysis was performed by using SPSS to determine the allele group/genotype frequency of the selected variants in the study population.
The genotyping results of G972R by RLFP-PCR showed the allelic frequency of G = 0.68 and R = 0.32 in controls while G = 0.71 and R = 0.29 in the cases. The minor R allele had a slightly higher frequency in the cases than the controls (OR = 0.86, CI 0.706–1.052, p = 0.17). The genotyping results of G1057D showed allelic frequency G = 0.74 and D = 0.26 in the controls while G = 0.961 and D = 0.29 in the cases. The minor D allele appeared to be a risk allele for this SNP although the difference in the allele frequencies was not statistically significant (OR = 1.55, CI 0.961–1.41, p = 0.108). The combined genotype analysis showed that the difference in the allele and genotype frequencies reached statistical difference between the cases and the controls as well as the odds ratio substantially increased when the R allele (G972R) was present together with D allele (G1057D) in any combination. When the association of single variants with the lipid traits was observed, only D allele (G1057D) showed significant association with TG, HDL and LDL, however when the analysis was repeated for combined genotypes using general linear model, many more significant associations between the genotype where D allele and R allele are together, were seen with many lipid traits.
In conclusion, the single nucleotide polymorphisms with low-modest effect size may not affect the phenotype individually but when in combination, the effect becomes stronger and more visible, therefore, for the SNP association studies, the more the number of SNPs included in the analysis, the more meaningful the results.
Diabetes is characterized by chronic hyperglycemia leading to the imbalance of fat, protein and carbohydrate metabolism due to defective insulin action, secretion or both . The prevalence of diabetes is increasing worldwide due to aging, population growth, urbanization, physical inactivity and obesity [2, 3]. The sequelae of diabetes include long term damage and dysfunctioning of various organs leading to their failure . The characteristic symptoms of diabetes mellitus involve polyuria, blurred vision, thirst and weight loss [5, 6]. The persons with either type 1 or type 2 diabetes require insulin treatment at some stage of the disease. Type 2 is the most prevalent form of diabetes that is characterized by insulin resistance [7, 8].
The insulin receptor substrates (IRSs) are the cytoplasmic receptors, important in insulin signaling . They also play a pivotal role in maintaining basic cellular functions like survival, growth and metabolism. After the activation of these proteins, several signaling cascades are activated. This leads to the activation of many downstream effectors. Finally, they gave the insulin signal to branching series of intracellular pathways. There are four members of the IRS family and are similar in their overall architecture .
When insulin binds to its receptor on cell surface, IRS-1, a docking protein is phosphorylated and activates phosphatidyl inositol 3-kinase (PI3-kinase). This initiates a cascade of intracellular signaling resulting in a variety of responses by the cell. One of these responses is the activation of GLUT (glucose transport) receptor resulting in increased uptake of glucose by the cell . The individuals with non-insulin dependent diabetes mellitus (NIDDM) have decreased expression and function of IRS-1 in the fat cells. The Gly972Arg polymorphism is the most commonly studied variant of IRS-1 . Downstream IRS-1, the main docking protein found in the cells is IRS-2 . It acts as the backup protein for the intracellular spread of the insulin signal including the activation of PI3-kinase. However, the hallmark of insulin resistance is that, for the activation of IRS-2, higher concentration of insulin is required. Gly1057Asp (G1057D) is a widely reported variant in IRS-2 gene whereby glycine (G) is replaced by aspartate (D) at 1057 position in the gene [14, 15]. These SNPs were selected for analysis in Pakistani subjects because (1) Pakistani population represents a unique ethnic group which allows the study of concentrated risk genetic markers (2) the selected SNPs have been investigated in the Caucasians so their analysis can provide the information on the relationship of these genetic variants with daibetes and (3) these SNPs have not been previously investigated in the Pakistani population and the current study is the first report of their study in our population. The biochemical risk factors included were lipid profile parameters namely total cholesterol (TC), very low density lipoprotein (VLDL), LDL-C, high density lipoprotein cholesterol (HDL-C) and triglyceride (TG) levels.
926 subjects (426 diabetic cases and 500 normal controls) were recruited from the local population using a case control design between July 2014 and April 2015. The subject recruitment, inclusion and exclusion criteria have been described previously [3, 16]. Inclusion criteria for diabetic subjects were (i) diabetes diagnosed according to etiologic classification of diabetes by the International Diabetes Federation (IDF) and (ii) confirmation that all the grandparents of the subjects are of Pakistani origin. The exclusion criteria were the presence of any infectious disease, conditions where phlebotomy is contraindicated, age below 10 years, body mass index (BMI) ≤ 18.5 kg/m2, pregnancy, handicapped/mentally disturbed individuals, obese (obesity defined as: BMI > 23 kg/m2 as overweight and BMI > 26 kg/m2 as obese, criteria used for Asian populations as described previously , cancer and ethnicity other than Pakistani. The controls were healthy subjects from the general population with normal blood sugar levels ( < 126 mg/dL random or < 99 mg/dL fasting) . All the subjects were genetically unrelated and gave written informed consent. The employed procedures were according to the Helsinki declaration and an ethical approval was obtained from the institutional ethical board.
Venous blood was collected from the subjects in the fasting state. 5 mL blood was collected from the median cubital vein using aseptic measures. 2.5 mL of blood was dispensed in an EDTA vial and used for DNA isolation and 2.5 mL was poured in a gel clot activator containing vial and used for biochemical analysis.
The DNA was isolated by manual method. The isolated DNA was stored at 4 °C or at − 20 °C for further use. The presence of DNA was determined by Agarose gel electrophoresis. After isolation, all DNA samples were quantified before amplification using Epoch Biotek microplate reader. The quantification was performed to equalize the final concentration of DNA samples up to approximately 10 ng/μL. The PCR conditions consisted of initial denaturation at 95 °C for 2 min, followed by 35 cycles of denaturation at 95 °C for 30 s, annealing at 57 °C for 45 s for G1057D variant and at 61 °C for 45 s for G972R variant, and extension at 72 °C for 1 min followed by a final extension at 72 °C for 5 min. PCR was done in advanced primus 96 thermal cycler. For Gly972Arg, the sequence of primers was: F: CTTCTGTCAGGTGTCCATCC, R: CTCTGCAGCAATGCCTGTTC. The 311 bp PCR product was digested with the help of BstN1 restriction enzyme and produced 3 fragments of 207 bp, 81 bp and 23 bp. For G1057D variant, the sequence of primers was: F: AGCTCCCCCAAGTCTCCTAA, R: GGCCACACCAAAAGCCATCT. The 291 PCR product was digested using the HhaI enzyme and produced 2 fragments one of 268 bp and other of 23 bp. The products were seen by running 2% agarose gel electrophoresis and observing the gel under the U.V light (Additional file 1: Figures S1, S2).
For conducting statistical analysis, MS Excel and SPSS (Statistical Package of Social Sciences) Package were used. Different variables such as Anthropometric, Biochemical and Genotype were analyzed by calculating mean and standard Deviation. Further the study population was tested for Hardy Weinberg Equilibrium (HWE), Chi-square test was also used to test the significance of difference of Allele and Genotype frequencies among diabetic and non-diabetic. To check the relationship of the SNPs with diabetes, t-test and logistic regression was conducted. Further the association of Polymorphism with Anthropometric and Biochemical parameters was tested through one way ANOVA. An adjusted p-value less than 0.025 was used as a significance cutoff for all such analysis due to inclusion of two SNPs at a time.
The general characteristics of the subjects have been described elsewhere  (Table 1). The study population was tested for Hardy–Weinberg equilibrium by a chi-squared test and was found to be in equilibrium (HWE p-value was 0.698 for G972R and was 0.377 for G1057D).
Allele/genotype frequencies of G972R and G1057D polymorphisms
After SNP genotyping, the genotypic and allelic frequencies were calculated. For Gly972Arg, in controls 230 (46%) individuals had GG genotype, 221 (44.2%) had GR genotype and 49 (9.8%) had RR genotype. Whereas, in cases 226 (53.05%) individuals had GG genotype, 153 (35.91%) had GD genotype and 47 (11.03%) had DD genotype. The odd ratio was 1.83 with 95% confidence interval (CI) 0.98–1.74. The p-value was 0.067 which shows an insignificant difference in allele/genotype frequencies between the cases and the controls. Allelic and genotypic frequencies of G1057D polymorphism showed that in controls 270 (54%) individuals had GG genotype, 200 (40.0%) individuals had GD genotype and 30 (6%) individuals had DD genotype. In diabetic subjects 231 (54.2%) individuals had GG genotype, 140 (32.8%) individuals had GD genotype and 55 (12.9%) individuals had DD genotype. The OR ratio of genotypic frequencies was 1.62 and 95% confidence interval was 1.22–2.16. The p-value of genotypic frequencies was 0.108 (Table 2). The combined genotype analysis revealed that the coexistence of any combination of the two variants with one or both risk alleles (GG + GD, GG + DD, GR + GD, GR + DD, RR + GD, RR + DD) had significant effect on the outcome as indicated by the odds ratio which increased proportionately when the heterozygous and homozygous genotypes of G972R was analyzes with heterozygous and homozygous genotypes of G1057D (Table 3).
Effect of the G972R and G1057D polymorphisms on biochemical parameters
Biochemical parameters were checked across the three genotypes for both variants individually by one way ANOVA as well as in combination by multivariate general linear model (GLM). For single variant analysis, G972R did not show statistically significant association with any of the lipid traits. For G1057D variant, the presence of risk allele significantly increased TG and LDL and lowered HDL although there was no effect on TC (Table 4). When the genotypes were analyzed in combination, the number of significant associations increased considerably. It was observed that when the common homozygous heterozygous and risk homozygous genotypes of G972R were examined against common homozygous, heterozygous and risk homozygous genotypes of G1057D. The results of this analysis are shown in Table 5. The table shows that the presence of homozygous risk genotype of G972R in any combination with G1057D (GG + GD, GG + DD, GR + GD, GR + DD, RR + GD, RR + DD) had the greatest effect on the change of lipid parameters (significant rise in TC, TG, LDL-C and a decrease in HDL-C).
The health care burden of diabetes is increasing at an alarming rate globally. The proper management for diabetes is therefore required to lessen its outcomes and to prevent the further complications [18, 19]. The Insulin Receptor Substrates (IRSs) of cytoplasmic adaptor proteins functioning in signal transduction in response to insulin binding at the cell surface consists of IRS-1, IRS-2, IRS-3, and IRS-4, named in order of identification [20, 21]. The identification of IRS-1 and IRS-2 gene polymorphisms associated with type2 diabetes was a landmark achievement in the field of diabetes genetics. The Gly972Arg is the most extensively studied polymorphism in IRS-1 gene, present in a highly conserved region of the protein that is involved in the phosphorylation of the substrate . In an in vitro cell culture study of human mesangial cells, the expression of G972R caused a significant reduction in the insulin-stimulated phosphorylation of IRS-1 and Akt kinase . The association studies of G972R SNP report positive as well as negative findings depending upon sample size, genotyping techniques and ethnicity. A previous study reported that carriers of Gly972Arg are at 25% higher risk of havig diabetes as compared to non carriers . Another study reported the association of Gly972Arg with type 2 diabetes in Mexican population . We found contradictory results compared to the Caucasian Dutch population where Gly972Arg did not appear to play a role in increased risk of type 2 diabetes . However, our results are in accordance with the previous published study to report the association of Gly972Arg with type 2 diabetes [27, 28].
The G1057D variant of IRS-2 gene has been reported in different groups but results are controversial as well. Certain studies have indicated that this variant is not linked to type 2 diabetes in Chinese, Finish, Turkish and Tunisian populations [29,30,31]. While some studies have shown that G1057D variant of IRS-2 gene predisposes mice to type 2 diabetes . Another study reported IRS-2 polymorphisms as one of the genetic causes of type 2 diabetes . Studies that included individuals from Kurdish Ethnic Group (in West Iran) and Asian Indians reported association of G1057D polymorphism of IRS-2 gene with type 2 diabetes [34, 35]. Similarly, genetic causes that lead to type 2 diabetes in Iranian and Indian population included G1057D polymorphism [36, 37].
In the current study, we determined the association of the variants of IRS-1 Gly972Arg (G972R) and IRS-2 Gly972Arg (G72R) with diabetes. The genotyping results of G972R showed that the minor R allele had a slightly higher frequency in the cases than the controls. The genotyping results of G1057D showed the minor D allele appeared to be a risk allele for this SNP although the difference in the allele frequencies was not statistically significant. The combined genotype analysis showed that the difference in the allele and genotype frequencies reached statistical difference between the cases and the controls as well as the odds ratio substantially increased when the R allele (G972R) was present together with D allele (G1057D) in any combination. When the association of single variants with the lipid traits was observed, only D allele (G1057D) showed significant association with TG, HDL and LDL, however when the analysis was repeated for combined genotypes using general linear model, many more significant associations between the genotype where D allele and R allele are together, were seen with many lipid traits.
There are certain limitations of our study as the sample size was relatively small and the samples were restricted to the Punjab province only. The mean age of controls and cases is significantly different that may cause bias in the study. There are other genes and polymorphisms that affect this association, but could not be included in the current investigation. However, despite these limitations, the results of the current study demonstrate the role of common variants in the progression of complex diseases like diabetes in the Pakistani population and show that if more variants are included in an investigation, more conclusive results can be obtained about the impact of the genetic markers on the phenotype. To our knowledge, the present study is the first attempt to report the association of IRS variants with diabetes in the local population.
In conclusion, there are various environmental factors and genes that work together in predisposition of an individual to complex diseases like diabetes. Moreover the single nucleotide polymorphisms can play an important role in the etiology of the complex diseases, however, is important to consider their combined effect on the phenotype.
Availability of data and materials
All the necessary information has been provided along with the manuscript, however, the corresponding author can be contacted for any information related to this paper.
Mooradian AD. Dyslipidemia in type 2 diabetes mellitus. Nat Clin Pract Endocrinol Metab. 2009;5(3):150–9.
Pinhas-Hamiel O, Zeitler P. The global spread of type 2 diabetes mellitus in children and adolescents. J Pediatr. 2005;146(5):693–700.
Shahzadi S, Shabana, Sarwar S, Shahid SU. Genome-wide implicated risk variants of TCF7L2 gene contribute to type 2 diabetes susceptibility by modulating serum lipids in Pakistani population. Int J Diabetes Dev Ctries. 2019;39(2):302–7.
Group CM-CTDS. Closing the gap: effect of diabetes case management on glycemic control among low-income ethnic minority populations. Diabetes Care. 2004;27(1):95–103.
Lowell BB, Shulman GI. Mitochondrial dysfunction and type 2 diabetes. Science. 2005;307(5708):384–7.
Saltiel AR, Kahn CR. Insulin signalling and the regulation of glucose and lipid metabolism. Nature. 2001;414(6865):799–806.
Devendra D, Liu E, Eisenbarth GS. Type 1 diabetes: recent developments. BMJ. 2004;328(7442):750–4.
Mohan V, Sandeep S, Deepa R, Shah B, Varghese C. Epidemiology of type 2 diabetes: Indian scenario. Indian J Med Res. 2007;125(3):217.
Yang Q, Graham TE, Mody N, Preitner F. Serum retinol binding protein 4 contributes to insulin resistance in obesity and type 2 diabetes. Nature. 2005;436(7049):356.
Muoio DM, Newgard CB. Molecular and metabolic mechanisms of insulin resistance and β-cell failure in type 2 diabetes. Nat Rev Mol Cell Biol. 2008;9(3):193–205.
Jakobsen SN, Hardie DG, Morrice N, Tornqvist HE. 5′-AMP-activated protein kinase phosphorylates IRS-1 on Ser-789 in mouse C2C12 myotubes in response to 5-aminoimidazole-4-carboxamide riboside. J Biol Chem. 2001;276(50):46912–6.
Mahmutovic L, Bego T, Sterner M, et al. Association of IRS1 genetic variants with glucose control and insulin resistance in type 2 diabetic patients from Bosnia and Herzegovina. Drug Metab Personal Ther. 2019. https://doi.org/10.1515/dmpt-2018-0031
Besse-Patin A, Jeromson S, Levesque-Damphousse P, Secco B, Laplante M, Estall JL. PGC1A regulates the IRS1: IRS2 ratio during fasting to influence hepatic metabolism downstream of insulin. Proc Natl Acad Sci. 2019;116(10):4285–90.
Mammarella S, Romano F, Di Valerio A, Creati B, Esposito DL, Palmirotta R, et al. Interaction between the G1057D variant of IRS-2 and overweight in the pathogenesis of type 2 diabetes. Hum Mol Genet. 2000;9(17):2517–21.
Vana DR, Adapa D, Prasad V, Choudhury A, Ahuja G. Diabetes mellitus types: key genetic determinants and risk assessment. Genet Mol Res. 2019. https://doi.org/10.4238/gmr16039952
Butt H, Shabana, Hasnain S. The C1431T polymorphism of peroxisome proliferator activated receptor γ (PPARγ) is associated with low risk of diabetes in a Pakistani cohort. Diabetol Metab Syndr. 2016;8(1):67.
Shabana, Shahid SU, Li KW, Acharya J, Cooper J, Hasnain S, et al. Effect of six type 2 diabetes susceptibility loci and an FTO variant on obesity in Pakistani subjects. Eur J Hum Genet 2016;24(6):904–10.
Samdani RA, Mahmood K. Web site of on-line index of PLB. 2003.
Mohan V, Shah S, Saboo B. Current glycemic status and diabetes related complications among type 2 diabetes patients in India: data from the A1chieve study. J Assoc Physicians India. 2013;61(1 Suppl):12–5.
Mardilovich K, Pankratz SL, Shaw LM. Expression and function of the insulin receptor substrate proteins in cancer. Cell Commun Signal. 2009;7(1):14.
de la Monte SM, Wands JR. Alzheimer's disease is type 3 diabetes—evidence reviewed. J Diabetes Sci Technol. 2008;2(6):1101–13.
Brancati FL, Kao WL, Folsom AR, Watson RL, Szklo M. Incident type 2 diabetes mellitus in African American and white adults: the Atherosclerosis Risk in Communities Study. JAMA. 2000;283(17):2253–9.
Thameem F, Puppala S, Schneider J, Bhandari B, Arya R, Arar NH, et al. The Gly (972) Arg variant of human IRS1 gene is associated with variation in glomerular filtration rate likely through impaired insulin receptor signaling. Diabetes. 2012;61(9):2385–93.
Jellema A, Zeegers M, Feskens E, Dagnelie P, Mensink R. Gly972Arg variant in the insulin receptor substrate-1 gene and association with type 2 diabetes: a meta-analysis of 27 studies. Diabetologia. 2003;46(7):990–5.
Burguete-Garcia AI, Cruz-Lopez M, Madrid-Marina V, Lopez-Ridaura R, Hernández-Ávila M, Cortina B, et al. Association of Gly972Arg polymorphism of IRS1 gene with type 2 diabetes mellitus in lean participants of a national health survey in Mexico: a candidate gene study. Metabolism. 2010;59(1):38–45.
Van Dam R, Hoebee B, Seidell J, Schaap M, Blaak E, Feskens E. The insulin receptor substrate-1 Gly972Arg polymorphism is not associated with Type 2 diabetes mellitus in two population-based studies. Diabet Med. 2004;21(7):752–8.
Porzio O, Federici M, Hribal ML, Lauro D, Accili D, Lauro R, et al. The Gly 972→Arg amino acid polymorphism in IRS-1 impairs insulin secretion in pancreatic β cells. J Clin Investig. 1999;104(3):357–64.
Stumvoll M, Fritsche A, Volk A, Stefan N, Madaus A, Maerker E, et al. The Gly972Arg polymorphism in the insulin receptor substrate-1 gene contributes to the variation in insulin secretion in normal glucose-tolerant humans. Diabetes. 2001;50(4):882–5.
Wang H, Rissanen J, Miettinen R, Kärkkäinen P, Kekäläinen P, Kuusisto J, et al. New amino acid substitutions in the IRS-2 gene in Finnish and Chinese subjects with late-onset type 2 diabetes. Diabetes. 2001;50(8):1949–51.
Ouederni TB, Sanchez-Corona J, Martinez SEF, Maiz HB, Skhiri HA, Abid HK, et al. The G1057D polymorphism of IRS-2 gene is not associated with type 2 diabetes and obese patients among ethnic groups in Tunisian population. Clin Biochem. 2009;42(10):1169–73.
Celik SK, Nurcan A, Kuçuk MU, Yalin S, Akbay E. The G1057D polymorphism of insulin receptor substrate-2 (IRS2) gene with type 2 diabetes in the Turkish population. Yeni Tıp Dergisi. 2014;31(3):168.
Withers DJ, Gutierrez JS, Towery H, Burks DJ, Ren J-M, Previs S, et al. Disruption of IRS-2 causes type 2 diabetes in mice. Nature. 1998;391(6670):900–4.
Singh S. The genetics of type 2 diabetes mellitus: a review. J Sci Res. 2011;55:35–48.
Bodhini D, Radha V, Deepa R, Ghosh S, Majumder P, Rao M, et al. The G1057D polymorphism of IRS-2 gene and its relationship with obesity in conferring susceptibility to type 2 diabetes in Asian Indians. Int J Obes. 2007;31(1):97–102.
Haghani K, Bakhtiyari S. The study on the relationship between IRS-1 Gly972Arg and IRS-2 Gly1057Asp polymorphisms and type 2 diabetes in the Kurdish ethnic group in West Iran. Genet Test Mol Biomark. 2012;16(11):1270–6.
Khodaeian M, Enayati S, Tabatabaei-Malazy O, Amoli MM. Association between genetic variants and diabetes mellitus in Iranian populations: a systematic review of observational studies. J Diabetes Res. 2015;2015:585917. https://doi.org/10.1155/2015/585917.
Radha V, Mohan V. Studies on genetics of diabetes in ethnic Indians. Textbook of diabetes mellitus RSSDI, vol. 1. 2nd ed. New Delhi: Jaypee Publications. 2008:349–57.
The authors thank to the University of the Punjab (PU16) for providing financial support.
Ethics approval and consent to participate
The study was approved by the institutional ethics committee (Ethical Committee, School of Biological Sciences, University of the Punjab, Pakistan) and all procedures were carried out in compliance with the Helsinki Declaration.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Figure S1.
Amplification of IRS-1 gene region containing G972R polymorphism. Last well contains a known sized DNA ladder (ThermoScientific SM#0321). Figure S2. Amplification of IRS-2 gene region containing G1057D polymorphism. Last well contains a known sized DNA ladder (ThermoScientific SM#0241).
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Ijaz, A., Babar, S., Sarwar, S. et al. The combined role of allelic variants of IRS-1 and IRS-2 genes in susceptibility to type2 diabetes in the Punjabi Pakistani subjects. Diabetol Metab Syndr 11, 64 (2019). https://doi.org/10.1186/s13098-019-0459-1
- Diabetes mellitus