1. Introduction
The UN and WHO have pointed to diabetes as the greatest health issue from a global epidemic prospect [1]. Global diagnosis criteria have already been synchronised between the WHO, the International Diabetes Federation (IDF), and the American Diabetes Association (ADA)[1]. The 2021 edition of IDF’s atlas estimates that 537 million adults between 20 and 79 years of age live with diabetes (10.5% prevalence), of whom 90% are type 2 patients [2]. Moreover, 240 million people are undiagnosed, meaning that almost half of the world’s diabetic patients are not aware of their ailment [2]. By 2045, about 783 million people will suffer from diabetes (12.2% prevalence).
In Mexico, the Secretary of Health, through the National Institute for Public Health (INSP), has carried out the National Health and Nutrition Survey (ENSANUT) for over 25 years. According to the latest data of ENSANUT 2021 COVID 19 [3], diabetes’ prevalence was 11.1% (IC95%=9.5-12.8) in a country whose population was 128.9 million by 2020.
Glycaemic control is fundamental to treat diabetes. The United Kingdom Prospective Diabetes Study (UKPDS), confirmed that it significantly reduced complication rates among DMT2 patients [4]. This study’s results established that the HbA1c control objective <7% (53 mmol) reported a microvascular complication decrease among both type 1 and type 2 patients, thus DMT2 patient handling should be focused on adequate control of hyperglycaemia and other risk factors [5, 6].
Metformin is classified as a biguanide and it is considered an essential drug by the WHO. It is the first-line therapy because of its efficacy, safety, and cost. Besides, it has a beneficial effect on HbA1c and weight reduction, and it can reduce cardiovascular and death risks [4].
Metformin’s primary effect is the hepatic inhibition of gluconeogenesis, however, its action mechanism is still debated [7]. Patients receiving a daily dose between 500 and 2,550 mg reach a plasmatic concentration ranging from ~10 up to 40 µM.
Metformin’s consensus plasmatic concentration in treated humans ranges between 0.1 mg/L up to 4 mg/L (~10 a 40 µM) [7]. Metformin is a positively charged hydrophilic molecule, it is mainly transported by organic cations 1-3 (OCT1-3), equilibrant nucleoside transporter 4 (ENT4), and multidrug and toxin extrusions 1 and 2k (MATE 1 and MATE 2-k) [7, 8]. In the intestinal lumen, metformin is absorbed by enterocytes through ENT4 and OCT3, then it is released into the intravascular space through OCT1 and distributed across the organism. Hepatocytes engage on metformin through OCT1 and OCT3, excreting it biliary through MATE1, although that is not the only elimination path, as renal tubular cells engage on metformin in the intravascular space by way of OCT2, from where it is subsequently excreted as urine via MATE1, MATE2-k, and OCT1 [8]. Thus, metformin distribution in humans is consistent with these transporters’ expression, where higher concentrations are found in the liver, kidneys, and intestines, whilst lower uptakes are found in peripheral organs [7]. The efficacy and toxicity of any drug are determined by its pharmacokinetic and pharmacodynamic balance. Inherent genetic variations contribute to this variability, since multiple genes codify proteins directly involved in the aforementioned pharmacological balance [9]. Transporters involved in metformin’s pharmacokinetics and pharmacodynamics belong to the SLC (Solute Carrier) family. This superfamily encompasses a few groups such as OCTs. OCT1 is a 554 aminoacid protein with a molecular weight of 61,154 Da; it is codified by gene SLC22A1, located in the long arm of chromosome 6, in region 25.3 (6q25.3), possessing a total of 12 exons [10]. Like many other members of the SLC22 family, it has 12 helixes or transmembrane domains (TMDs), including a large extracellular loop between domains 6 and 7 [11]. SLC22A1 is a highly polymorphic gene, and its polymorphisms may induce an altered OCT1 function, which affects metformin’s pharmacokinetics and the eventual response to it [8, 12]. More than 23 genetic variants or single nucleotide polymorphisms (SNPs) have been identified for SLC22A1. This study focuses on one exonic (due to direct protein changes) and one intronic polymorphism (due to possible changes in the protein’s translation). The exonic polymorphism in question is rs72552763, since it may represent another 3: rs3567514, rs34305973, and rs35191146, located in exon 7, which would be found by TMD 9 in the final protein [11, 13]. The intronic OCT1 polymorphism is rs622342, a variant which affects transcription rate decrease [14]. Studying the frequence and functionality of these two polymorphisms could lead to identifying pharmacogenetic biomarkers of relevance in DMT2 therapeutic response.
2. Results
2.1. Clinical and Biochemical Data of the Studied Population
Out of 103 patients, 59 fulfilled inclusion criteria. As shown in
Table 1 [15], patients were grouped by gender. Statistically significant differences only regarded height (
p<0.001). HbA1c median was found at 6.30 (5.80 – 7.30) and 6.40 (6.10 – 7.85) respectively for men and women. Uncontrolled (HbA1c >7%) proportions were 29.42% for men and 26.19% for women, where no statistical significance regarding gender was found.
2.2. Allelic and Genotypic Frequencies of the Studied Polymorphisms among DMT2 Patients Undergoing Metformin Treatment
Table 2 shows the allelic and genotypic frequencies, as well as the
p value for the Hardy-Weinberg equilibrium (HWE) test, which was > 0.05 across all analysed polymorphisms in every single case.
2.3. HbA1c Control by SNPs in SLC22A1, SLC22A2, and SLC22A3.
We assessed frequency distribution among controlled (HbA1c< 7%) and uncontrolled (HbA1c≥ 7%) patients. We analysed frequency distribution across the two groups based on both polymorphism genotype and dominant genotypic model for every SNP; yet, no statistical significance was found (
Table 3).
2.4. Studied Biomarker Values by Genotypic Model
We conducted a statistical inference across the different SNP genotypes and grouped them according to an individualised genotypic model (
Table 4). There were significant differences in HbA1c levels by rs72552763 (
p=0.0220) and rs622342 (
p=0.0093) in
SLC22A1; however, this result found no echo in fasting glucose levels. The pertaining post hoc trials revealed a statistical difference (
p<0.05) by GAT/del (
Figure 1) and CC in rs72552763 and rs622342 (
Figure 2) respectively.
There was no significant difference by SNPs in
SLC22A2 and
SLC22A3 regarding HbA1c levels across genotypes. Grouping by dominant genotypic model (GAT/GAT vs GAT/del + del/del) revealed statistical significance in both HbA1c (p=0.0058) (
Figure 3) and fasting glucose levels (p=0.0260), whilst the dominant genotypic model of rs622342 (A/A vs A/C+C/C) revealed statistical significance only in HbA1c levels (p=0.0108) (
Figure 4), not in fasting glucose (p=0.1551). Models applied to rs316019 and rs2076828 yielded no significant difference.
2.5. Metformin Plasmatic Concentrations by Genotype and Genotypic Model
Metformin plasmatic concentrations (ng/ml) were assessed by comparing their levels across all four polymorphisms (
Table 5). There were no significant differences. The dominant genotypic model assessment revealed no significant difference either.
3. Discussion
There are official reports which indicate an HbA1c non-control prevalence of about 75% among DMT2 patients [16]; yet, our Mexican-Mestizo sample undergoing metformin monotherapy reported a non-control prevalence of 29.42% among men and 26.19% among women (p= 1.000), where no frequency difference regarding gender was found. This study reported an overall non-control prevalence of 27.11%, a very inferior rate with respect to ENSANUT 2012 [16]. It is important to clarify that ENSANUT reports consider diabetic populations undergoing any kind of treatment, not just metformin monotherapy. Our data suggest that the disease’s non-control proportion in the case of metformin monotherapy is lower than presently reported by ENSANUT. This may be due to our sample’s median diagnosis period (3 years), which could be expressing data on the disease’s first stages. However, according to central tendency measurements revealed by our sample’s description (table 1), at least 50% of the studied population has a BMI higher than 30 Kg/m2, which represents a non-control risk as per previously reported data [17].
OCTs have already been reported to directly affect metformin pharmacokinetics. Moreover, genetic variations caused by polymorphisms codifying SLC22A1, SLC22A2, and SLC22A3, have been proven to alter metformin pharmacokinetics, affecting its hepatic action. Polymorphism rs72552763 in SLC22A1, where 3 bases (ATG) in codon 420 (p.M420del) are deleted, has been studied in Mexican populations by Reséndiz-Abarca [18] and Menjivar et al. [19]. On rs622342 Reséndiz found no significant differences between alleles AA, AC, and the minor allele CC regarding HbA1c among patients treated with metformin for less than 3 years. In a linear model association analysis of rs622342, adjusted by gender, age, disease duration and waist perimeter, Reséndiz did find a significant difference by CC, where HbA1c value (p<0.001) increased after 12 months of observation [18]. Nonetheless, Menjívar’s study of rs72552763 found no significant difference between fasting glucose (mg/dl) (p=0.368) and HbA1c levels (p=0.181).
Even so, Menjívar reported differences in rs72552763 genotypic proportions when comparing controlled and uncontrolled patients depending on HbA1c levels: GAT/GAT= 33/107, GAT/del= 19/39, and del/del=16/16, (p=0.011).
These results differ from those reported by the present work, because we found no significant differences between controlled and uncontrolled patients across rs72552763 genotypes (p= 0.4678). Still, our analysis of HbA1c suggests that its lowest levels are found in GAT/GAT as compared to the minor allele (p= 0.0220).
It shall be mentioned that Menjívar’s study [19] included patients undergoing a combined treatment of metformin and glibenclamide, whereas our study includes only metformin monotherapy, like Reséndiz-Abarca [18]. Our results coincide with Christensen’s work [20], which reported that reduced function alleles in OCT1 SNPs, including rs72552763, is associated with a statistically significant decreasing trend in metformin’s stationary concentration and an increase in HbA1c levels during the cohort’s observation (p= 0.043). We analysed metformin’s plasmatic concentrations, and while its median was higher by GAT/GAT, 504.04 (131.65 – 900.16)(ng/ml) with respect to deletion carriers of s72552763, 236.48 (86.76 – 785.24)(ng/ml), there was no statistical significance (p=0.2384).
The Audit trail in Tayside, Scotland (DARTS) [21], a cohort survey including 1531 patients treated with metformin, reported that genotypes of rs72552763 have no effect on HbA1c levels or the biomarker’s decrease after 42 months; but in that study the proportion of del/del was 3 times less than our sample’s (15.25% vs 4.7%) [22].
On the other hand, in 2015 Umamaheswaran 2015 [22] reported that recently diagnosed DMT2 patients from Southern India, were treated with lifestyle adequation and metformin monotherapy, starting at 500 mg/day, which increased to 2.5 g/day after 12 weeks. When studying metformin response across AA, AC, and CC, in rs622342, a comparison between respondent (HbA1c decrease ≥ 0.5 %) and non-respondent patients (fasting glucose >180 mg; HbA1c decrease <1 % after 12 weeks) by way of a logistic regression based on AA as the reference group found that minor allele C presented a greater non-response risk in genotypes AC (p= 0.011) OR= 3.50 (1.39 - 8.84; 95% CI) and CC (p=0.033) OR= 5.60 (1.24 - 25.80; 95% CI). A dominant genotypic model ( AA vs AC+CC), found that minor allele C is at greater non-response risk (p= 0.003) OR= 3.85 (1.61 - 9.19; 95% CI). These results [22] coincide with ours, as we found higher HbA1c levels among patients carrying the minor allele C of rs622342 in SLC22A1 (p= 0.0108).
4. Materials and Methods
4.1. Study Design and Sample Description (DMT2 Patients)
Study design
This is an observational, transversal, clinical and analytical trial, conducted at HRAEI between April 2018 and April 2019. Patients’ histories were compiled collecting clinical and biochemical data found in their electronic record. These included age, disease duration, HbA1c levels, etc. A peripheral venous blood sample was collected from every patient using two tubes containing ethylenediaminetetraacetic acid (EDTA).
Patient recruitment
Out of a previously studied sample of 103 Mexican-Mestizo DMT2 patients (diagnosed by ADA and WHO criteria) from Hospital Regional de Alta Especialidad de Ixtapaluca (HRAEI) [17], those 59 who fulfilled the metformin monotherapy marker were selected for this new study.
4.2. Clinical Evaluation
The procedure has been previously reported in detail [17] (pp16-17).
4.3. Genotyping Procedure
The procedure has been previously reported in detail [17] (pp17-18).
4.4. Plasmatic Metformin Determination
The procedure has been previously reported in detail [17] (pp18-19).
4.5. Statistical Analyses
Inferential analysis.
The pertaining statistical analyses were conducted using R language (available at
https://www.r-project.org/). Apposite code lines were generated to programme every analysis. Either Kolmogorov-Smirnov or Shapiro-Wilk tests were performed for the population description analysis. According to their distribution, quantitative variables were analysed using either T-Student’s or U Mann-Whitney’s tests, accounting for p<0.05 as a statistically significant value. Qualitative variables were analysed through Pearson’s Chi-square test. The Kruskal-Wallis test was used to conduct comparisons between more than two groups, where p<0.05 was considered a significant value as well.
Genotypic and Allelic Frequency Analysis
Allelic frequencies were counted, and respective expected values were calculated for each genotype. A p value > 0.05 defined the frequencies in Hardy–Weinberg equilibrium as calculated through Pearson’s Chi square test.
5. Conclusion
Summarising, although our sample is small, the obtained data allow us to conclude that polymorphisms of rs72552763 and rs622342 in SLC22A1, seem to relevantly affect metformin response among Mexican-Mestizo DMT2 patients. Still, longitudinal studies in a similar population undergoing metformin monotherapy are necessary to observe this effect’s consistency in time.
Author Contributions
Conceptualization, ADIEL ORTEGA-AYALA, CARLOS BARTOLO-MONTIEL, Fernando De Andrés, Adrián LLerena and Juan Arcadio Molina-Guarneros; Data curation, ADIEL ORTEGA-AYALA, CARLOS BARTOLO-MONTIEL, Fernando De Andrés, Adrián LLerena and Juan Arcadio Molina-Guarneros; Formal analysis, ADIEL ORTEGA-AYALA, CARLOS BARTOLO-MONTIEL, Fernando De Andrés, Adrián LLerena and Juan Arcadio Molina-Guarneros; Funding acquisition, Adrián LLerena and Juan Arcadio Molina-Guarneros; Investigation, Adrián LLerena and Juan Arcadio Molina-Guarneros; Methodology, ADIEL ORTEGA-AYALA, CARLOS BARTOLO-MONTIEL, Fernando De Andrés, Adrián LLerena and Juan Arcadio Molina-Guarneros; Resources, Adrián LLerena and Juan Arcadio Molina-Guarneros; Software, ADIEL ORTEGA-AYALA and CARLOS BARTOLO-MONTIEL; Validation, ADIEL ORTEGA-AYALA, CARLOS BARTOLO-MONTIEL, Fernando De Andrés, Adrián LLerena and Juan Arcadio Molina-Guarneros; Writing – original draft, ADIEL ORTEGA-AYALA, CARLOS BARTOLO-MONTIEL, Fernando De Andrés, Adrián LLerena and Juan Arcadio Molina-Guarneros; Writing – review & editing, ADIEL ORTEGA-AYALA, Fernando De Andrés, Adrián LLerena and Juan Arcadio Molina-Guarneros.
Funding
Supported by the Programa de Apoyo a Proyectos de Investigación e Innovación Tecnológica (Project No. IN205121) of the Dirección General de Apoyo a Personal Académico, Universidad Nacional Autónoma de México. The present study was also partially supported by AEXCID, Junta de Extremadura, Spain (“Hampi-Yachay” AEXCID 22IA002).
Institutional Review Board Statement
The study was conducted in accordance with the Declaration of Helsinki. Protocol NR-004-2016 was approved on the 7th of April, 2016 by the Research and Ethics Commissions of the Hospital Regional de Alta Especialidad de Ixtapaluca, the Secretary of Health, and the Research Division of UNAM’s Faculty of Medicine 001/SR/2016 01/17/2017.
Informed Consent Statement
All patients who agreed to voluntarily be involved in the present study provided written informed consent. No identifying information, including names, initials, date of birth, or hospital numbers, images, or statements are included in the manuscript.
Data Availability Statement
The data presented in this study are available on request via the corresponding author. These data are not publicly available because the patients and researchers are bound to an agreement establishing that only the head of the study and Mexican health authorities shall have access to them, in accordance with the presidential decree of 16th April, 2015, sanctioning the General Law on Transparency and Access to Public Information.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Standl, E., Khunti, K., Hansen, T. B., & Schnell, O. The global epidemics of diabetes in the 21st century: Current situation and perspectives. European Journal of Preventive Cardiology (2019). 26(2_suppl), 7–14. [CrossRef]
- Sun H., Pouya P., Karuranga S., Pinkepank M., Ogurtsova K., Duncan B.B., Stein C., Basit A., Chan J.C.N., Mbanya J.C., Pavkov M. E., Ramachandaran A, Wild S.H., James S., Herman W.H., Zhang P., Bommer C., Kuo S., Boyko E.J., Dianna J., Magliano D.J. IDF DiabetesAtlas:Global, regional and country-level diabetes prevalence estimates for 2021 and projections for 2045. Diabetes Research and Clinical Practice (2022) 183 I09II9, 1-13. [CrossRef]
- Romero-Martínez M, Barrientos-Gutiérrez T, Cuevas-Nasu L, Bautista-Arredondo S, Colchero A, Gaona-Pineda E.B, Lazcano-Ponce E, Martínez-Barnetche J, Alpuche-Aranda C, Rivera-Dommarco J, Shamah-Levy T. Metodología de la Encuesta Nacional de Salud y Nutrición 2020 sobre Covid-19. Salud Publica Mex (2021) May 3;63(3):444-451. [CrossRef]
- Intensive blood-glucose control with sulphonylureas or insulin compared with conventional treatment and risk of complications in patients with type 2 diabetes (UKPDS 33). UK Prospective Diabetes Study (UKPDS) Group. Lancet (1998). 352: 837-853,.
- Diabetes Care Volume 45, Supplement 1, January (2022). Matthew C. Riddle, MD. [CrossRef]
- Diagnóstico y Tratamiento Farmacológico de la Diabetes Mellitus Tipo 2 en el Primer Nivel de Atención. Guía de Evidencias y Recomendaciones: Guía de Práctica Clínica. México, Instituto Mexicano del Seguro Social. (2018). 55. http://imss.gob.mx/profesionales-salud/gpc (accessed on July 18th 2023).
- LaMoia T. and Shulman G.I. Cellular and Molecular Mechanisms of Metformin Action Endocr Rev. 2021 Feb; 42(1): 77–96. [CrossRef]
- Staiger, H., Schaeffeler, E., Schwab, M., & Häring, H. U. Pharmacogenetics: Implications for modern type 2 diabetes therapy. Review of Diabetic Studies, (2015). 12(3–4), 363–376. [CrossRef]
- Huang, C., & Florez, J. C. Pharmacogenetics in type 2 diabetes: Potential implications for clinical practice. Genome Medicine (2011). 3(11). [CrossRef]
- Weizzman Institute of Science. (2022). Gene cards SLC22A1. https://www.genecards.org/cgi-bin/carddisp.pl?gene=SLC22A1(Accessed July 20th 2023).
- Koepsell, H., Lips, K., & Volk, C. Polyspecific organic cation transporters: Structure, function, physiological roles, and biopharmaceutical implications. Pharmaceutical Research (2007) 24(7), 1227–1251. [CrossRef]
- Lai Y. Organic anion, organic cation and zwitterion transporters of the SLC22 and SLC47 superfamily (OATs, OCTs, OCTNs and MATEs). In Transporters in Drug Discovery and Development. (2013). [CrossRef]
- Srijib, G. PharmGKB summary: very important pharmacogene information for SLC22A1. Pharmacogenetics and Genomics (2014). 23(1), 1–7. [CrossRef]
- Becker M.L., Visser L.E., van Schaik, R.H.N., Hofman A., Uitterlinden A.G., Stricker B.H.C. Interaction between polymorphisms in the OCT1 and MATE1 transporter and metformin response. Pharmacogenet. Genom. (2010) 20, 38–44. Available online: https://pubmed-ncbi-nlm-nih-gov.pbidi.unam.mx:2443/19898263/ (accessed on 2 May 2022).
- WHO. Physical Status: The Use of and Interpretation of Anthropometry, Report of a WHO Expert Committee. 1995. Available online: https://apps.who.int/iris/handle/10665/37003 (accessed on July 21th 2023).
- Shamah Levy T, Cuevas Nasu L, Morales Ruan MC, Mundo Rosas V, Méndez Gómez-Humarán I, Villalpando Hernández S. Profile of the Health and Nutritional Status of Older Adults in Mexico. 2012 National Health and Nutrition Survey. J Frailty Aging. 2013;2(4):184-91. [CrossRef]
- Ortega-Ayala A, Rodríguez-Rivera NS, Andrés F. d, LLerena A, Pérez-Silva E, Espinosa-Sánchez AG, Molina-Guarneros JA Pharmacogenetics of Metformin Transporters Suggests No Association with Therapeutic Inefficacy among Diabetes Type 2 Mexican Patients. Pharmaceuticals (2022). 15(7):774. [CrossRef]
- Reséndiz-Abarca C., Flores-Alfaro E., Suárez-Sánchez F., Cruz M., Valladares-Salgado A., Alarcón-Romero L.D.C., Wacher-Rodarte N.A., Gómez-Zamudio J.H. Altered Glycemic Control Associated With Polymorphisms in the SLC22A1 (OCT1) Gene in a Mexican Population With Type 2 Diabetes Mellitus Treated With Metformin: A Cohort Study. J. Clin. Pharmacol. (2019), 59 (10) 1384–1390. Available online: https://pubmed-ncbi-nlm-nih-gov.pbidi.unam.mx:2443/31012983/ (accessed on Jul 18th 2023).
- Menjivar M., Sánchez-Pozos K., Jaimes-Santoyo J., Monroy-Escutia J., Rivera- Santiago C., Granados-Silvestre A., Ortiz-López M.G. Pharmacogenetic Evaluation of Metformin and Sulphonylurea Response in Mexican Mestizos with Type 2 Diabetes. Curr. Drug Metab. (2020), 21, 291–300. Available online: https://pubmed-ncbi-nlm-nih-gov.pbidi.unam.mx: 2443/32407269/ (accessed on May 2nd 2022).
- Christensen MH , Brasch-Andersenb C , Greene H, Nielsen F , Damkierc P , Beck-Nielsen H and Brosena K. The pharmacogenetics of metformin and its impact on plasma metformin steady-state levels and glycosylated hemoglobin A1c. Pharmacogenetics and Genomics (2011) 21(12):p 837-850. [CrossRef]
- Morris A.D., Boyle D.R. , MacAlpine R, Emslie-Smith A, Jung R.T., Newton RW, MacDonald T. M. for the DARTS/MEMO Collaboration. The diabetes audit and research in Tayside Scotland (darts) study: electronic record linkage to create a diabetes register. BMJ (1997) 315:524.
- Zhou, K., Donnelly, L. A., Kimber, C. H., Donnan, P. T., Doney, A. S. F., Leese, G., Hattersley, A. T., McCarthy, M. I., Morris, A. D., Palmer, C. N. A., & Pearson, E. R. Reduced-function SLC22A1 polymorphisms encoding organic cation transporter 1 and glycemic response to metformin: A GoDARTS study. Diabetes, (2009) 58(6), 1434–1439. [CrossRef]
- Umamaheswaran G, Praveen R.G. Damodaran S. E. Das A.K. Adithan C. Influence of SLC22A1 rs622342 genetic polymorphism on metformin response in South Indian type 2 diabetes mellitus patients. Clin Exp Med (2015) 15:511–517.
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).