Genetic Polymorphisms Associated with Insulin Resistance Risk in Normal BMI Indians

1. Study Participants

The research enrolled 191 participants (90 men and 101 women) of Indian ethnicity residing in India, aged 18-65 years, between January 2022 and December 2023. HOMA2-IR was used to categorise participants into cases and controls. Individuals with HOMA2-IR > 2 were classified as cases (57 participants), and those with HOMA2-IR < 2 (134 participants) were categorised as controls. The HOMA2-IR was calculated using the online tool provided by the Medical Science Division of The University of Oxford¹⁸. An electronic registration form was used to gather the demographic data of the participants, including their self-reported weight (kg) and height (cm). The formula for calculating body mass index (BMI) was weight (kg)/height (m^2)¹⁹. For adult males, the formula for Body Fat Percentage (BFP) was 1.20 x BMI + 0.23 x Age - 16.2, for an adult female, the formulas for Body Fat Percentage (BFP) were 1.20 x BMI + 0.23 x Age - 5.420²⁰ and Fat Mass was BFP x Weight x 0.0120²⁰. To ensure participant safety and ethical consideration, the study included only individuals free from: cancer history, cardiovascular or renal failure, mental illnesses, and pregnancy or lactation. Every participant gave written, informed consent in accordance with the guidelines in the 2008 updated Declaration of Helsinki²¹. The study protocol received approval from the Answer Genomics Ethical Review Committee. Figure 1 provides a detailed graphical representation of the study’s structure and key components, outlining the research methodology, participant selection, and analysis framework.

• Download figure
• Open in new tab

Figure 1: Study Flowchart

2. Laboratory Measurements

Following a 12-hour fast, participants’ blood was drawn for analysis of important metabolic markers, including, high-density lipoprotein cholesterol (LDL-C), total cholesterol (TC), triglycerides (TG), fasting plasma glucose (FPG), and low-density lipoprotein cholesterol (LDL-C),. The assessments were conducted using Beckman DxC 700 AU for all markers except HbA1c, which was analyzed with Tosoh G-8, and fasting insulin levels measured by Beckman UniCel DxI 800, adhering to the protocols provided by the manufacturers. The healthy reference ranges for these markers were established as follows: FPG (70-100 mg/dl), TC (0-100 mg/dl), TG (0-150 mg/dl), HDL-C (40-60 mg/dl), and LDL-C (0-100 mg/dl). To determine insulin sensitivity, the Homeostasis Model Assessment of Insulin Resistance (HOMA2-IR) was applied, calculated from fasting plasma glucose (in mmol/l) and fasting insulin (in mU/l) using the formula: fasting glucose times fasting insulin/22.5¹⁰. Participants with a HOMA2-IR score above 2 were categorized as insulin resistant, while those with scores below 2 were considered insulin sensitive. The clinical characteristics of the study subjects and the distribution of these groupings are detailed in Table 2.

3. Genotyping and Single-nucleotide polymorphism(SNP) selection

This study investigated potential genetic influences on insulin resistance by analyzing participants’ DNA. The extraction of DNA was conducted on blood specimens utilizing the Qiagen blood extraction kit, acclaimed for its capability to isolate genomic DNA of superior quality from blood samples. The process of genotyping SNPs was executed with the aid of the Illumina Infinium Global Screening Array (GSA) V3 platform, recognized for its extensive representation of genetic variance throughout the genome²². The Illumina iScan system was used for genotyping, while Genome Studio V2 software was used for interpreting the raw intensity data, quality control, and data export²³.

The study focused on four genes (INTS10, LINC01427 - LINC00261, SENP2 and SLC22A11) involved in metabolism and potentially influencing insulin resistance, as indicated in Table 1. The selection of SNPs was based on a comprehensive literature review of scientific literature identifying genetic markers associated with metabolic traits.

4. Statistical analysis

The evaluation of statistical differences commenced by examining the variations in age, glucose, and lipid metrics such as LDL cholesterol, HDL cholesterol, total cholesterol, triglycerides, and HbA1c between two groups. Given the data’s non-parametric characteristics, the Mann-Whitney U test (utilising the Scipy library) was utilised for these analyses¹⁷. Gender distribution across the groups was evaluated using the Chi-square test, with detailed results presented in Table 2.

Utilising PLINK software²⁴ for quality control, the study filtered genotype data by removing samples with low genotyping rates, excluding SNPs with significant missing data, and eliminating rare variants for their minimal statistical power. Heterozygosity checks were also conducted to ensure accurate genotype distribution. Subsequently, a study was carried out using SNPstats software²⁵ to ascertain the correlation between the described four SNPs and insulin resistance. The analysis yielded odds ratios (ORs) together with their 95% confidence intervals (CIs). Table 3 summarises the results, which show significant SNP markers under different genetic inheritance models.

In the control group, Hardy-Weinberg equilibrium was maintained with a significance level set at 0.05. The assessment of risk was conducted through the calculation of odds ratios (OR) along with their 95% confidence intervals (95% CI). The study also evaluated the linkage disequilibrium (LD) between the SNPs, with the LD coefficient D being determined via the use of Snpstats software ²⁵. Additionally, power analyses were conducted based on power and sample size calculations²⁶. Each SNP underwent a Bonferroni adjustment was carried out for each SNP to rule out the possibility of a type I error brought on by multiple tests.

Subject characteristics

In this study, all 191 participants were genotyped, with their clinical features and lipid metabolism indicators detailed in Table 2. The investigation revealed no meaningful statistical disparities in age or gender between the case and control groups, nor among males and females within those groups (p>0.05). Yet, notable distinctions were observed in BMI (0.01), Fat Mass (0.007), HBA1C (0.01), Glucose Fasting (0.001), and Serum Triglycerides (0.02), with the exceptions of BFP, Total Cholesterol, and Serum HDL Cholesterol when comparing cases to controls.

Association of the 4 SNPs with IR

In this study, the relationship between IR and four SNPs in the Indian population’s normal weight individuals was examined. To ensure the study had sufficient statistical power, a sample size calculation was performed. This calculation aimed to detect significant associations, considering a significance level (alpha) of 0.05.

The frequency of genotypes for all four SNPs among both cases and controls was evaluated to examine their compliance with Hardy-Weinberg equilibrium (HWE) through a Fisher exact test. None of the SNPs demonstrated a significant deviation from HWE (p > 0.05).suggesting the absence of genotyping errors and population stratification, the results are mentioned in Table 3.

We performed chi-square tests to evaluate the association between genotypes and allele frequencies of each SNP with phenotype (IR). After Bonferroni correction for multiple testing (p < 0.0125), all four SNPs displayed significant associations with IR:

Chi-squared analyses were performed to investigate the relationship between the genotypic and allelic variations of four SNPs (rs920590, rs7274134, rs6762208, and rs2078267) within two groups distinguished by their HOMA2 IR levels (>2 and <2).The obtained p-values indicate significant associations between SNP genotypes and insulin resistance status.The allelic frequency difference between the cases and control were also found to be significant (p< 0.05)

• For rs920590 (INTS10), individuals carrying the T/T genotype had an elevated risk of IR compared to those with C/C-C/T genotypes [Odds Ratio (OR) = 4.01; 95% Confidence Interval (CI): 1.55-10.34; p < 0.0014] under the recessive model.

• rs7274134 (LINC01427 - LINC00261), individuals carrying the C/T-T/T genotypes had an elevated risk of IR [Odds Ratio (OR) = 2.60, 95% CI: 1.37-4.96, p < 0.0032] according to the dominant model.

• rs6762208 (SENP2), individuals carrying the A/C-C/C genotypes were associated with an increased risk of IR [Odds Ratio (OR) = 3.11, 95% CI: 1.55-6.24, p < 0.0013] under the dominant model.

• rs2078267 (SLC22A11), individuals carrying the T/T genotype showed a higher risk of IR [Odds Ratio (OR) = 3.72, 95% CI: 1.22-11.28, p < 0.0091] under the recessive model.

INTS10

Protein coding gene INTS10 (Integrator Complex Subunit 10) is a part of the Integrator (INT) complex, which is a complex engaged in the transcription pathway leading to gene expression²⁷. In a diabetic retinopathy study done on mice, it was found that the transcription of INTS10 in the inner retinal cells of diabetes induced mice was upregulated early, suggesting that INST10 contributes to diabetic retina due to high insulin resistance²⁸. In another GWAS study carried out on the African population, a polymorphism was reported to be associated with BMI that lies between INTS10 and LPL genes, which the study suggests is a BMI loci²⁹. LPL gene has Lipoprotein lipase activity that is crucial in amassing the triglycerides from the blood and storing them in the adipocytes, which is seen more in women than in men³⁰. Serum LPL mass is linked to metabolic syndrome, coronary heart disease, and ultimately insulin sensitivity, according to a study by Onat et al.³¹.

LINC01427 - LINC00261

With more than 200 nucleotides, long non-coding RNAs (lncRNAs) are a type of significant regulatory RNAs that are involved in a variety of cellular functions in both healthy and pathological conditions.^{32, 33}. They affect translation processes, RNA splicing, stability, chromatin structure and function, and gene transcription. Additionally, they play a role in the formation and regulation of organelles and nuclear condensates³⁴. LncRNAs are classified by their genomic locations. LINC01427 and LINC00261 are long intergenic ncRNAs (lincRNAs), which are located between protein-coding genes³⁵.

It has been demonstrated that lncRNA dysregulation occurs in both animal and human pancreatic islets. Numerous long noncoding RNAs (lncRNAs) are connected to the development of insulin resistance (IR) and are engaged in multiple stages of the insulin manufacturing process^{36, 37}. Dysregulated levels of certain lncRNAs in peripheral blood were positively connected with IR, impaired glucose management, inflammation, and transcriptional markers of senescence in a study comparing type 2 diabetes (T2D) patients with healthy controls. Even after controlling for confounding variables, there was a substantial correlation established between this dysregulation and type 2 diabetes³⁸. A study by Junpei et al. revealed that a dysregulated lncRNA expression profile is associated with beta cell dysfunction in humans³⁹. Beta cell dysfunction leads to impaired insulin secretion, while diminished responsiveness of target tissues to normally secreted insulin leads to insulin resistance. The failure of beta cell function exacerbates insulin resistance, thereby accelerating the progression of type 2 diabetes^{40, 33}.

Research on mice has revealed that a specific long non-coding RNA (lncRNA), called lncSHGL, suppresses gluconeogenesis and lipogenesis in the liver³⁹. When lncSHGL is restored in mice, it improves hyperglycemia, insulin resistance, and fat accumulation in obese diabetic mice. Conversely, inhibiting lncSHGL in the liver leads to increased fasting blood sugar and fat buildup in normal mice. These results imply that long non-coding RNAs regulate inflammation and the synthesis of fat, which are critical factors in the initiation and development of insulin resistance and glucose regulation⁴¹.

SENP2

SENP2 belongs to the family of proteolytic enzymes SUMO-specific proteases (SENPs) that reverse the effects of sumoylation, which is a post-translational modification that controls the activity and viability of proteins⁴². Koo YD et al., investigated the role of SENPs in energy metabolism. It was found that when a specific muscle cell line was treated with saturated fatty acids, it led to the increase of SENP2 expression, thereby activating the fatty acid oxidation. Overexpression of SENP2 in skeletal muscle increased high-fat diet-induced obesity and IR, suggesting a potential therapeutic target for IR-linked metabolic disorders⁴². The study carried out by Chung SS et al reveals that SENP2 is crucial in adipogenesis regulation. Its expression increases upon adipocyte differentiation and is dependent on protein kinase A activation⁴³. Knockdown of SENP2 in mice leads to the reduction of adipogenesis, PPARgamma, and C/EBPalpha mRNA levels. Sumoylation of C/EBPbeta reverses this effect. Overexpression of C/EBPbeta overcomes knockdown’s inhibitory effect on adipogenesis⁴³. Krapf SA et al., studied how (SENP2) can influence fatty acid and glucose metabolism in primary adipocytes. The researchers found that silencing of SENP2 reduced glucose uptake, oxidation and lipogenesis, while increasing lipid oxidation.This suggests that SENP2 plays a crucial control on energy metabolism in primary human adipocytes⁴⁴.

SLC22A11

Studies suggest that elevated uric acid levels can lead to metabolic syndrome and insulin resistance, with fructose raising uric acid levels⁴⁵. The protein-coding gene SLC22A11, involved in transporting and excreting uric acid, is linked to nephrolithiasis, uric acid, and gout. Its genes OAT4 and URAT1 encode renal urate transporters, significantly affecting Serum Urate⁴⁶. Insulin regulates cellular processes through various signaling pathways, with serum- and glucocorticoid-inducible kinases (sgk) acting as mediators⁴⁷. Sgk increases transport activity and OAT4 expression by overcoming Nedd4–2 inhibitory effects⁴⁸. Through the action of ubiquitin ligase Nedd4–2, insulin increases OAT4 expression and transport activity. Nedd4–2 specific siRNA is used to knock down Nedd4–2, weakening regulation. Because insulin competes with sgk2 rather than working through sgk2 to regulate OAT4, the effects of insulin and sgk2 are cumulative⁴⁷. IGF-1, a liver hormone, is crucial for metabolism, growth, and development. It increases OAT4 transport activity and SUMOylation in kidney-derived cells, but PKB-specific inhibitors block its regulation⁴⁹.

The investigation of optimal inheritance models reveals varying modes of genetic inheritance contributing to IR susceptibility. While the recessive model is suitable for rs920590 and rs2078267, the dominant model is favoured for rs7274134 and rs6762208. These findings highlight the complexity of genetic factors influencing IR and the need to consider diverse genetic inheritance patterns in understanding the genetic architecture of metabolic disorders.

This study also emphasizes the significance of clinical parameters in understanding IR development. Serum LDL cholesterol, serum triglycerides, and fat mass exhibit significant associations with IR, indicative of the multifaceted nature of metabolic dysfunction. Combining genomic and clinical data highlights the interaction between hereditary and environmental factors in metabolic health and offers a thorough framework for understanding the pathophysiology of IR⁵⁰.

The validation of genetic markers associated with IR offers insights for personalized therapeutic strategies. By understanding the genetic determinants underlying IR, clinicians can tailor treatment approaches to individual genetic profiles, enhancing the management of metabolic disorders. Identification of novel genetic targets opens avenues for targeted intervention and prevention strategies, paving the way for precision medicine in metabolic health.

While this study provides significant insights into IR’s genetic underpinnings, further investigations are warranted to unravel the intricate interplay between genetic and environmental factors contributing to metabolic dysfunction. Longitudinal studies incorporating larger cohorts and diverse populations are needed to validate and extend these findings. Functional studies are required to understand the molecular mechanisms underlying the observed associations and can provide deeper insights into IR’s pathogenesis. Integrating multi-omics approaches offers a comprehensive understanding of the complex interactions driving metabolic health and disease, facilitating the development of targeted therapeutic interventions and precision medicine strategies aimed at mitigating the burden of metabolic disorders.