Association of the Reduced Function Met420del Polymorphism of SLC22A1 with Metformin-Induced Gastrointestinal Intolerance in Ethiopian Patients with Type 2 Diabetes Mellitus

Introduction

Metformin is indicated as a first-line glucose lowering agent for patients with type 2 diabetes mellitus (T2DM),¹ though it has been frequently associated with gastrointestinal symptoms than most other oral glucose lowering agents.^2,3 Moreover, although metformin is still recommended as a drug of choice for patients with T2DM,^4–6 there are increasing reports demanding a more personalized approaches in its use.^7,8 Previous studies reported that the variability to metformin treatment associated incidence of side effects is related to differences in individual genetic profile.^9,10

Metformin diffuses poorly across biological membranes as it is a hydrophilic molecule at physiological pH.¹¹ Variations in the gene encoding for one of its main transporters (the organic cation transporter-1(OCT1)),^12,13 SLC22A1, could thus potentially affect metformin associated gastrointestinal adverse effects.¹⁴ Furthermore, several studies postulated that decreased metformin transport through OCT1 that is expressed in enterocytes might raise local metformin concentration in the intestine, leading to metformin’s gastrointestinal side effects.^15–17

In addition, in recent years, the alteration in the gut microbiome composition is increasingly recognized as a potential mediator of metformin-associated gastrointestinal adverse effects.¹⁸ To this effect, a comparative study on the gut microbiota profile of metformin-tolerant and metformin-intolerant patients with T2DM revealed that metformin-induced shift in the composition and function of the gut microbiota potentially mediate the drug intolerance.¹⁹ Moreover, recent studies showed that the common metformin induced gastrointestinal adverse effects have been attributed to gases (such as CO2 and H2S) produced by the gut microbiota, seemingly through metformin-induced selection and metabolic modelling of Escherichia spp. and Akkermansia muciniphila.^20–22 Thus, the decreased metformin transport through OCT1 associated with the reduced functional variants of SLC22A1 might increase metformin concentration in the gut, inducing selective alteration in the metabolic profile of the human gut microbiota and a subsequent shift to potentially mediate metformin associated gastrointestinal adverse effect.

However, there are inconsistent reports on the association between the role of individual SLC22A1 variants and gastrointestinal side effects.¹⁴ To address this issue, the reduced function methionine deletion at codon 420 (Met420del) variant of SLC22A1 (rs72552763) was selected for this study since its allele frequency is higher in Ethiopian population as compared to other African population.^23–25 Furthermore, Met420del is the most studied functional variant of SLC22A1 gene for its influence on metformin pharmacokinetics though there were inconsistent reports based an extensive survey of recent literature. In addition, this variant is generally observed in combination with a rare loss of function polymorphism, Cys88Arg (rs55918055), that causes improper membrane localization of OCT1, resulting in metformin gastrointestinal intolerance.¹⁵ Thus, in this study, we investigated the association of an individual Met420del polymorphism of SLC22A1 with metformin induced gastrointestinal intolerance in Ethiopian patients with T2DM.

Materials and Method

Study Design and Setting

The study setting was described in our recent study, which reported the association of Met420del polymorphism of SLC22A1 with metformin response.²³ In this study, a retrospective observational study was conducted on T2DM patients on metformin treatment for <3 years to assess the association of SLC22A1 (rs72552763) polymorphism with metformin induced gastrointestinal intolerance.

Selection of the Study Participants

All patients were diagnosed according to the WHO criteria²⁶ and information on comorbidities, medication use, clinical data, biochemical test parameters, and metformin gastrointestinal intolerance were obtained from clinical records. Socio-demographic data, metformin gastrointestinal intolerance, life style, and medication adherence were obtained from questionnaire interview. Information on gastrointestinal intolerance was obtained both from medical charts and through interview as the information was not available in the charts for some participants and also to confirm for those available.

Eligibility Criteria

Inclusion Criteria

Age 18–65, recently diagnosed (less than 3 years) unrelated T2DM patients on metformin monotherapy within the first 6 months of treatment (for intolerant group) and on metformin treatment for more than 6 months (tolerant group), and willing to give consent to participate in the study were included in the study.

Exclusion Criteria

Pregnant or lactating women, patients with renal or hepatic insufficiency (serum creatinine >1.2mg/dl for renal insufficiency; serum glutamic-pyruvic transaminase or serum glutamic-oxaloacetic transaminase >40IU/L for hepatic insufficiency), active smokers (an adult who smoked even one puff during the past seven days), patients with malignancies or thyroid disorders or depression depending on the documented medical history, patients not adherent to medications, patients with chronic gastrointestinal diseases (including chronic liver disease, cholelithiasis, chronic pancreatitis, inflammatory bowel disease and gastroduodenal ulcer), and patients unable to recall being on metformin or having the gastrointestinal adverse effects were excluded.

Sample Size Estimation and sampling

Comparison of two proportion formula²⁷ was used to determine whether the rs72552763 polymorphism of SLC22A1 exists between metformin intolerant and metformin tolerant groups, considering the following assumptions: p₁ = 0.30,²⁸ 95% CI, 5% margin of error and a power of 80%. Accordingly, the calculated sample size was set at 21 for each group in the present study. Adding a 10% to account for potential dropouts, the final sample size was 47 (24 for tolerant and 23 for the intolerant group). Purposive sampling technique was used to enroll study subjects.

Collection of Data

In this study, initial preliminary identification of potentially eligible study participants was done based on the attending physician’s review of the medical chart and/or active patient compliant to the attending physician working at the diabetic follow-up clinic. Subsequently, a written consent was obtained from each potential study participant before collection of data. In order to reduce potential source of bias, data were collected by trained nurses working on diabetic follow-up clinic after the participants agreed to written informed consent. An interviewer administered questionnaire and a chart review were used for data collection.

Thus, relevant information was collected from their medical charts and confirmed through oral interview whether gastrointestinal intolerance was a reason for discontinuation and/or not increasing metformin dose over the first 6 months of therapy. However, for those participants who did not have a recorded metformin-associated gastrointestinal adverse effects in their medical charts but actively complaining about the adverse effects to their attending physician, we confirmed through oral interview whether gastrointestinal intolerance was a reason for discontinuation and/or not increasing metformin dose. Accordingly, patients were categorized as metformin tolerant and intolerant groups. Patients meeting the following were taken as intolerant: i) Patients who were on metformin monotherapy and stopped receiving metformin immediate release (IR) formulation within the first 6 months of treatment and switched to another oral hypoglycemic agent (including metformin slow-release formulations) within 6 months of the last metformin IR prescription and reported gastrointestinal adverse effects as the reason for switching;^15,29 ii) Patients who were on metformin monotherapy and could not increase their metformin IR dose >500 mg daily despite an HbA1c >7.0% because of the reported gastrointestinal adverse effects.³⁰

Metformin-tolerant individuals were defined as those treated with >2000 mg of metformin daily for more than 6 months (excluding modified-release formulations of metformin) and reported no adverse effects.¹⁵ A daily dose of metformin was defined as the last prescribed dose for intolerant patients and as an average dose in the first 6 months (after satisfying tolerance definition) of metformin treatment for tolerant patients.

A four-item Self-Reported Measure of Adherence scale were used to assess medication adherence and individuals who responded “NO’ for all the four questions were considered adherent to medication.³¹ Furthermore, individuals who modified their lifestyle as recommended for greater than 3 days in the last seven successive days were considered adherent to diet.³² Likewise, study participants who accomplished at least 150 min per week (3 days) of moderate-intensity physical exercise were considered as active or inactive if perform any lesser.³³

Genotyping

Blood collection, DNA extraction and genotyping of the rs72552763 polymorphism were described in our recent study.²³ Briefly, genotyping of the rs72552763 polymorphism was performed using the TaqMan^® Drug Metabolism Enzyme Genotyping Assay (Assay ID: C__34211613_10) from Applied Biosystems (Carlsbad, CA, USA) and ChamQ Geno-SNP Probe Master Mix (Vazyme Biotechnology, Singapore). The qPCR was conducted using the StepOnePlus™ Real-Time PCR System. The genotyping quality was assessed by looking at the allelic discrimination plots for the color and number of discrete clusters as the software makes genotyping call based on VIC/FAM dye ratios and cluster orientation (Figure 1). The quality value was assigned manually in each well for each genotype call to be 100%. Furthermore, a no template control sample was also run in each assay to confirm whether proper genotype call was made or not.

Figure 1 Allelic discrimination plot for the delGAT/GAT polymorphic SNP rs72552763.

Statistical Analysis

Statistical Package for Social Sciences (SPSS) Version 26 for Windows (IBM Corps, Armonk, New York, NY, USA) was utilized for statistical analyses. Categorical variables were described as percentages. Continuous variables were described as medians and interquartile range (IQR). The normality of continuous variables was determined by Kolmogorov–Smirnov and Shapiro–Wilk tests and all the data were not normally distributed. For continuous dependent variables, associations were assessed by Independent-Samples Mann–Whitney U-test. The associations between categorical variables were assessed by a Chi-square test and 95% confidence interval. Logistic regression model analysis was also performed for categorical dependent variables. The Hosmer–Lemeshow goodness-of-fit test was done to assess the model fit (Hosmer–Lemeshow statistic > 0.05). A p-value less than 0.05 was considered statistically significant. In order to evaluate the association of rs72552763 polymorphism with metformin induced gastrointestinal adverse effects, the number of metformin-intolerant patients was assessed. In this study, the Hardy–Weinberg equilibrium (HWE) test was performed using Gene-Calc software, which uses the Chi-square test to determine whether or not observed genotype frequencies are consistent with Hardy-Weinberg expectations.³⁴

Ethics Approval and Informed Consent

The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board (IRB) of the College of Health Science, Addis Ababa University (protocol number:118/20/SoP) and IRB of St. Paul’s Hospital Millennium Medical College (RN: DM23/576). National Ethical approval was granted by the National Research Ethics Review Committee, Ministry of Education (MOE), Ethiopia (RN:02/246/572/22). All participants received information about the process and purpose of the study in their local language (Amharic) and provided with a written informed consent. All collected data were anonymized and held confidential.

Results

Characteristics of the Study Participants

In this study, 24 metformin tolerant and 23 metformin intolerant individuals with T2DM were recruited. The Independent-Samples Mann–Whitney U-test showed no significant difference in the median values of BMI and age at diagnosis between the two groups. Obviously, there was a significant difference in the median metformin daily doses between the metformin tolerant and intolerant groups. The median daily dose of metformin was 1000 mg in metformin intolerant group while it was 2000 mg in the tolerant group. Furthermore, the Independent-Samples Mann–Whitney U-Test showed no significant difference in the median values of baseline HbA1c and FBG between metformin tolerant and intolerant study participants (Figure 2). However, there were significant differences in sex and level of physical activity between metformin tolerant and intolerant groups in a chi-square test analysis (Table 1).

Table 1 General Characteristics of the Metformin Intolerance Study Participants (n = 47)

Figure 2 Independent-Samples Mann–Whitney U-test between metformin tolerant and intolerant groups on their baseline FBG (panel (A) and HbA1c levels (panel (B).

Genotype Frequency

The minimum allele frequency (MAF) of the 3-base pair (GAT) deletion mutation at rs72552763 was 9.6%, whereas that of the wild-type GAT (G) allele was 90.4% in the current study. No patient was homozygous for the rs72552763 deletion (del) allele in this study (Table 2). The SNP rs72552763 was in line with the principles of the HWE with p-value of 0.984.

Table 2 Genotype and Allele Frequency Distribution of the rs72552763 Polymorphism Across Intolerance Study Participants (n = 47)

Association of rs72552763 with Metformin Intolerance

In the logistic regression analyses, our study found that there was no significant association between metformin intolerance and genotypes of the SNP rs72552763 (Table 3). In the chi-square test analyses, the alleles of rs72552763 polymorphism were also not significantly associated with metformin intolerance (Table 4). Furthermore, no significant intergroup difference in age were observed between the study participants who were metformin tolerant and intolerant in a logistic regression analysis (Table 3). By contrast, the female gender (OR=4.019, 95% CI (1.048–15.416), p<0.05) and physical inactivity (OR=5.455, 95% CI (1.414–21.035), p<0.05) were found to be significantly associated with increased metformin intolerance (Table 3).

Table 3 Logistic Regression Analysis on the Predictors of Metformin Intolerance Among Study Participants (n = 47)

Table 4 Chi-Square Test Analysis on the Effect of the Alleles of the rs72552763 Polymorphism on Metformin Intolerance Among Study Participants (n = 47)

Discussion

In this study, metformin gastrointestinal intolerance was defined in 47 patients with T2DM. Since we involved an individual variant of SLC22A1, the Met420del polymorphism (rs72552763), to associate with metformin gastrointestinal intolerance in our study, intolerance phenotype was defined to include common and/or severe gastrointestinal intolerance, as both share a common underlying mechanism despite the difference in the intensity and duration of the symptoms.²⁹ Moreover, as no study participant was homozygous for the deletion allele in this study, comparison was made between the patient homozygous for the wild-type GAT allele (G_G genotype) and the patient heterozygous for the deletion allele in the association study.

In our study, no significant association was observed between the Met420del variant of SLC22A1 (rs72552763) and metformin intolerance. Indeed, this result is in accordance with earlier finding, which showed that Met420del variant had no association with metformin side-effects, although the same study reported the association of common gastrointestinal side effects of metformin therapy in an individual with reduced function variants of SLC22A1 (rs628031 and rs36056065).¹⁴ The findings from two previous studies also showed that individual variants (including rs72552763) were not associated with the occurrence of side effects, though the number of SLC22A1 reduced-function alleles was associated with common gastrointestinal adverse effects to metformin²⁹ and individuals carrying two reduced function SLC22A1 alleles had higher chances of severe intolerance compared with individuals with one or no reduced function alleles.¹⁵ In addition, in the same study¹⁵ individuals with two SLC22A1 reduced-function alleles receiving treatment with known OCT1 inhibitors exhibited a nearly four times more gastrointestinal intolerance. Collectively, these findings indicate that inadequate SLC22A1 transport appears to be a major contributor to the metformin gastrointestinal intolerance, although there is lack of consistency in the association between individual reduced function SLC22A1 variants and metformin gastrointestinal intolerance. Thus, it is plausible to assume that the combined effect of the reduced function variants of SLC22A1 is valuable to predict metformin induced gastrointestinal adverse effects rather than the individual variants. However, as our study involved only a single reduced functional variant of SLC22A1, Met420del variant, this assumption needs to be confirmed by further studies involving more than one reduced functional variants of the gene. Besides, we assume that the absence of association between Met420del polymorphism and metformin induced gastrointestinal intolerance in our study might be related to the absence of participant who is homozygous for the deletion allele and, most importantly, due to the smaller sample size. Thus, a more complete investigation of SLC22A1 variants would be required to fully assess the effect of the gene on metformin induced gastrointestinal intolerance as several variants with a more severe loss of function have been described.

Furthermore, our study found that the female gender was associated with more than 4 times greater chance of being intolerant to metformin. This finding adds to the previous knowledge that indicated sex-related differences in metformin intolerance.^14,15,29 However, further studies are needed to assess the underlying factors for the potential sex differences in metformin intolerance. Moreover, there was a significant difference in the level of physical activity between metformin tolerant and intolerant groups. It is not clear why physical inactivity is related to metformin gastrointestinal intolerance. Further prospective studies are required to confirm this observation. Age appeared to have influence on the occurrence of side effects depending on the time of T2DM diagnosis. Evidence for this assertion comes from the observation that whilst age had significant effect on the occurrence of side effects in patients with a longer duration of T2DM,^14,15 it was found to have no association in patients with newly diagnosed T2DM,²⁹ which is concordant with our findings. This might have to do with the fact that both metformin tolerant and intolerant patients had similar age in the current study, as all the study participants were recently diagnosed T2DM individuals.

One of the major limitations of this study is the relatively small sample size. We used the minimum sample size estimated by comparison of two proportion formula. However, more reliable results could have been obtained if we expanded the sample size. Thus, the findings of this study need to be cautiously interpreted and require validation in larger cohort studies involving multiple variants of SLC22A1 gene reported to affect metformin induced gastrointestinal intolerance. We excluded patients on medications, such as proton pump inhibitors and cimetidine that could potentially affect OCT1 in the present study. As a result, how such medications affect the outcome is lacking. Future studies including such factors into account may help to fully establish the effect of the gene on metformin induced gastrointestinal intolerance.

Conclusions

In this study, rs72552763 polymorphism was not significantly associated with metformin induced gastrointestinal intolerance in Ethiopian patients with T2DM. The present study is the first to investigate the association of rs72552763 with metformin intolerance in the Ethiopian population with T2DM. Thus, the finding of this study might provide additional important data that contributes to knowledge synthesis for the Ethiopian and African genetic diversity in terms of metformin associated adverse effects. In our study, we also found that the female gender and physical inactivity were risk factors for metformin induced gastrointestinal intolerance.