Implications of genetic variations, differential gene expression, and allele-specific expression on metformin response in drug-naive type 2 diabetes

Purpose Metformin is widely used to treat type 2 diabetes mellitus (T2DM) individuals. Clinically, inter-individual variability of metformin response is of significant concern and is under interrogation. In this study, a targeted exome and whole transcriptome analysis were performed to identify predictive biomarkers of metformin response in drug-naive T2DM individuals. Methods The study followed a prospective study design. Drug-naive T2DM individuals (n = 192) and controls (n = 223) were enrolled. T2DM individuals were administered with metformin monotherapy and defined as responders and non-responders based on their glycated haemoglobin change over three months. 146 T2DM individuals were used for the final analysis and remaining samples were lost during the follow-up. Target exome sequencing and RNA-seq was performed to analyze genetic and transcriptome profile. The selected SNPs were validated by genotyping and allele specific gene expression using the TaqMan assay. The gene prioritization, enrichment analysis, drug-gene interactions, disease-gene association, and correlation analysis were performed using various tools and databases. Results rs1050152 and rs272893 in SLC22A4 were associated with improved response to metformin. The copy number loss was observed in PPARGC1A in the non-responders. The expression analysis highlighted potential differentially expressed targets for predicting metformin response (n = 35) and T2DM (n = 14). The expression of GDF15, TWISTNB, and RPL36A genes showed a maximum correlation with the change in HbA1c levels. The disease-gene association analysis highlighted MAGI2 rs113805659 to be linked with T2DM. Conclusion The results provide evidence for the genetic variations, perturbed transcriptome, allele-specific gene expression, and pathways associated with metformin drug response in T2DM.