Allele frequency and genotype distribution of a common variant in the 3´-untranslated region of the SLC22A3 gene in patients with type 2 diabetes: Association with response to metformin
Maryam Ghaffari-Cherati1, Abdolkarim Mahrooz2, Mohammad Bagher Hashemi-Soteh2, Seyyedeh Raheleh Hosseyni-Talei1, Ahad Alizadeh3, Saeid Mofid Nakhaei4
1 Department of Clinical Biochemistry and Genetics, Faculty of Medicine, Mazandaran University of Medical Sciences, Sari, Iran
2 Department of Clinical Biochemistry and Genetics, Faculty of Medicine, Mazandaran University of Medical Sciences; Immunogenetic Research Center, Mazandaran University of Medical Sciences, Sari, Iran
3 Department of Epidemiology and Reproductive Health, Reproductive Epidemiology Research Center, Royan Institute for Reproductive Biomedicine, ACECR, Tehran, Iran
4 Medical Diagnostic Laboratory, Bu.Ali Hospital, Mazandaran University of Medical Sciences, Sari, Iran
|Date of Submission||11-Feb-2016|
|Date of Decision||01-Jun-2016|
|Date of Acceptance||23-Jun-2016|
|Date of Web Publication||18-Oct-2016|
Department of Clinical Biochemistry and Genetics, Faculty of Medicine, Mazandaran University of Medical Sciences, Km 17 Khazarabad Road, Sari
Source of Support: None, Conflict of Interest: None
Background: Organic cation transporter 3 (OCT3) is an excellent transporter for metformin, which is used as first-line therapy for type 2 diabetes (T2D). OCT3 genetic variants may influence the clinical response to metformin. This study aimed to determine the genotype and allele frequency of OCT3-564G>A (rs3088442) variant and its role in the glycemic response to metformin in patients with newly diagnosed T2D. Materials and Methods: Based on the response to metformin, 150 patients were classified into two groups: Sixty-nine responders (decrease in glycated hemoglobin [HbA1c] values by more than 1% from the baseline) and 81 nonresponders (decrease in HbA1c values <1% from the baseline). HbA1c levels were determined by chromatography. The variant OCT3-564G>A was genotyped using polymerase chain reaction - based restriction fragment length polymorphism. Results: The genotypes frequencies were 51.3% GG, 36% AG, and 12.7% AA. Allele frequency of major allele (G) and minor allele (A) in OCT3-564G>A variant was found to be 0.69 and 0.31, respectively. Fasting glucose, HbA1c, body mass index, and lipid profile in both GG genotypes and GA + AA group decreased significantly after 3 months of metformin therapy compared with baseline (P < 0.05). In both responders and nonresponders, HbA1c and fasting glucose levels were lower in patients with the GA + AA genotype than in those with the GG genotype; however, the differences were not statistically significant (P > 0.05). Conclusion: The A allele frequency (which may be a protective allele against coronary heart disease) in the Iranian diabetic patients was lower than Iranian, Caucasian and Japanese healthy populations. Metformin is useful in improving the lipid profile, in addition to its impacts in glycemic control, and these effects are regardless of OCT3-564G>A variant.
Keywords: Genetic variation, glycated hemoglobin, metformin, organic cation transporter 3, type 2 diabetes
|How to cite this article:|
Ghaffari-Cherati M, Mahrooz A, Hashemi-Soteh MB, Hosseyni-Talei SR, Alizadeh A, Nakhaei SM. Allele frequency and genotype distribution of a common variant in the 3´-untranslated region of the SLC22A3 gene in patients with type 2 diabetes: Association with response to metformin. J Res Med Sci 2016;21:92
|How to cite this URL:|
Ghaffari-Cherati M, Mahrooz A, Hashemi-Soteh MB, Hosseyni-Talei SR, Alizadeh A, Nakhaei SM. Allele frequency and genotype distribution of a common variant in the 3´-untranslated region of the SLC22A3 gene in patients with type 2 diabetes: Association with response to metformin. J Res Med Sci [serial online] 2016 [cited 2020 Jan 26];21:92. Available from: http://www.jmsjournal.net/text.asp?2016/21/1/92/192508
| Introduction|| |
Diabetes is a major threat to public health worldwide, and its prevalence has rapidly increased and even has influenced the younger population.,, The World Health Organization estimates that about 366 million people worldwide will develop diabetes by 2030, the majority of who develop type 2 diabetes (T2D). Lifestyle negative changes such as obesity, weight gain, reduced physical activity, high-calorie diet and smoking contribute to the increased prevalence of the disease.,
Metformin is the most widely used antidiabetic medication among different medications administered to treat the disease, and it is considered as the first-line oral treatment for T2D., In general, the mechanisms of metformin action include reduction of hepatic glucose production, increased glucose uptake in skeletal tissue, suppression of the intestinal absorption of glucose, and improvement of insulin sensitivity., Metformin reduces triglyceride (TG) levels, free fatty acids, and low-density lipoprotein (LDL) in plasma and increases high-density lipoprotein (HDL). Due to slight hydrophobic property, metformin needs a transporter to cross plasma membrane into the cell. A family of transporters called organic cation transporters (OCTs) is among transporters with an important role in bioavailability, clearance, and pharmacodynamics of metformin. OCTs have three members: OCT1-3. OCT3 is coded by SLC22A3 gene located on chromosome 6q26-q27. OCT3 has a high-capacity for binding to the substrate and an extensive tissue distribution in the body, particularly in muscles and adipose tissue., This transporter contributes to the transfer of endogenous physiological amino compounds such as catecholamines. Like two other members of OCT family, OCT3 is an excellent transporter for metformin. Various factors such as nature of the disease and organ function affect drug response, but estimates show that genetics can dictate 20–95% of variability in different responses to the same medication. Hence, OCT3 genetic variants may influence the clinical response to metformin. OCT3-564G>A is a noncoding variant in 3´-untranslated region (3´-UTR) of OCT3., The displacement G >A creates a binding site for microRNA called miR-147. The presence of miR-147 leads to negative expression of OCT3 gene. In other words, this variant contributes to the gene expression of OCT3 and affects the final production of mRNA.
So far, few studies have been conducted on the relationship between OCT3 and metformin response. Therefore, similar studies on various variants of OCT3 should be conducted in different ethnic populations. Hence, in this study, we evaluated whether the variant OCT3-564G>A influences metformin response and also the distribution of genotypes of this variant in Iranian patients with T2D.
| Materials and Methods|| |
This observational study conducted on 150 newly diagnosed patients with T2D (according to the WHO criteria) treated with metformin from March 2014 to August 2015. The study was performed in Northern province of Iran, Mazandaran. The patients were followed for 3 months and during this period, they received 1000 mg of metformin (Aria pharmaceutical company, Iran) twice a day. None of the patients were receiving antidiabetic medication prior to their diabetes diagnosis. Information about medical history, personal habits, demographic parameters, and medication use was obtained through a questionnaire. Patients were divided into two groups based on the response to metformin: (1) Responder group with a reduction in HbA1c values by more than 1% compared to before taking metformin and (2) nonresponder group with a reduction in HbA1c values <1% compared to before taking metformin. Subjects with type 1 diabetes, and previous history of renal failure, chronic diseases, autoimmune, and liver diseases was not included in the study. The study protocol was approved by the Ethical Committee at Mazandaran University of Medical Sciences (ethnic number: 93-1303) and informed consent was obtained from all participants.
Sampling and laboratory assays
Blood samples were taken from subjects after an overnight fast (for 10-12 h). Sera were isolated by low-speed centrifugation, and the serum aliquots were stored at −70°C until use. The HbA1c levels were determined by boronate affinity technique (Axis-Shield PoC AS, Oslo, Norway; accuracy, failure <5%). HDL-cholesterol (HDL-C), total cholesterol (TC), triglycerids (TGs), fasting blood glucose (FBG), alanine aminotransferase (ALT) levels were measured using an auto-analyzer (Prestige, Japan). Levels of LDL-cholesterol (LDL-C) were determined by the Friedewald formula. All of the tests were performed before and after 3 months of metformin therapy.
Genotyping of OCT3-564G>A variant
The OCT3-564G>A variant was genotyped using the restricted fragment length polymorphism analysis after polymerase chain reaction amplification. The designed primers used for this variant were: F 5'-AGATTGCATGGAGGATGAAC-3' and R 5'-TGTTACAGAGGAGGTGGACG-3'. DNA was amplified with an initial melting temperature of 93°C for 3 min, followed by 35 cycles of 50 s at 93°C, 40 s at 55°C, and 40 s at 72°C, with a final extension step of 5 min at 72°C. Amplification products from each sample (171 bp) were digested with the restriction enzyme AciI (thermoscientific) at 37°C for 16 h, and resulted in 67 and 104 bp fragments, which were subjected to electrophoresis on a 2.5% agarose gel. Wild-type (GG) patients were identified by the presence of 67 and 104 bp fragments. Heterozygous (GA) patients were identified by the presence of 171, 104, and 67 bp fragments, while the presence of an undigested fragment (171 bp) was the basis for the identification of mutants (AA).
The distribution normality of the variables was checked by the Kolmogorov–Smirnov test. To analyze the differences between the parametric variables, we used independent t-test or paired t-test. Mann–Whitney U-test used to analyze the nonparametric variables. The Chi-square test was used for checking Hardy–Weinberg equilibrium. Deviation from the equilibrium may be indicative of problematic assays. A P < 0.05 was accepted as statistically significant. Statistical analyses were performed by SPSS (version 16.0) software (SPSS Inc, Chicago, IL).
| Results|| |
[Table 1] indicates the changes in study parameters after 3 months of metformin therapy. As can be seen in the table, most of these parameters including systolic and diastolic blood pressure (P < 0.001), body mass index (BMI) (P < 0.001), fasting glucose (P < 0.001), and HbA1c (P < 0.001) were significantly reduced after 3 months of treatment. Compared with baseline levels before treatment, TG (P < 0.001), TC (P = 0.001), and LDL-C (P < 0.001) were improved following 3 months of treatment. HDL-C levels were increased (P = 0.101) after 3 months of treatment compared with baseline levels.
|Table 1: Change in the study variables from baseline to 3 months of metformin treatment (n=150)|
Click here to view
OCT3-564G>A genotypes were in Hardy–Weinberg equilibrium in the study patients. The frequency of GG, AG, and AA genotypes were 51.3%, 36%, and 12.7%, respectively. The allele frequency of G and A was found to be 0.69 and 0.31, respectively.
As shown in [Table 2], the majority of the studied parameters such as HbA1c, fasting glucose, BMI, and lipid profile in both GG genotypes and GA + AA group decreased significantly after 3 months of metformin therapy compared with baseline. Although there was an increase in HDL-C levels in both patients with the GA + AA genotype and GG genotype after 3 months of treatment with metformin, statistically significant differences were not found.
|Table 2: Change in the study parameters from baseline to 3 months of metformin therapy according to the genotypes of organic cation transporter 3-564G>A|
Click here to view
In the current study, subjects were divided into two groups: Metformin responders (n = 69) and metformin nonresponders (n = 81). There were no statistically significant differences between responders and nonresponders with respect to the study parameters, except HbA1c levels before treatment (results not shown). When we analyzed the study parameters with respect to metformin response according to OCT3-564G>A genotypes, our findings showed no statistically significant differences between responders and nonresponders with respect to the majority of the parameters of this study [Table 3] and [Table 4]. In both nonresponders [Table 3] and responders [Table 4], fasting glucose levels were lower in patients with the GA + AA genotype than in those with the GG genotype; however, the differences were not statistically significant. Furthermore, our results showed that in GA + AA group compared with GG group, reduction in average HbA1c values was higher in responders than in nonresponders [Figure 1], although were not found the statistically significant differences among the groups.
|Table 3: Change in the study variables after 3 months of metformin therapy in metformin nonresponders according to the genotypes of organic cation transporter 3-564G>A|
Click here to view
|Table 4: Change in the study variables after 3 months of metformin therapy in metformin responders according to the genotypes of organic cation transporter 3-564G>A|
Click here to view
|Figure 1: Change in glycated hemoglobin levels after 3 months of metformin therapy in metformin responders and nonresponders based on the genotypes of OCT3-564G>A|
Click here to view
| Discussion|| |
To the best of our knowledge, there are few studies on the allele frequency of OCT3-564G>A in different populations. In our study, the major allele frequency of G and minor allele frequency of A were 0.69 and 0.31, respectively, in patients with T2D. Hengen et al. reported that G and A allele frequencies were 0.53 and 0.47, respectively, in healthy Caucasians in the United States. Moreover, in the study of Aoyama et al., G and A allele frequencies were 0.51 and 0.49, respectively, in a healthy Japanese population. It should be noted that according to another study we are doing, the frequency of A allele in a healthy Iranian population (n = 152) was 0.38. It seems that minor allele frequency of A in Iranian patients with T2D and also in healthy Iranian population is lower than that of healthy Caucasian and Japanese populations. According to an investigation recently published by Li et al., AA genotype of OCT3-564G>A variant compared to GG genotype significantly reduces the risk of coronary heart disease (CHD) in a Chinese population. In other words, they reported that this variant can contribute to reducing the risk of CHD through its supporting role against inflammatory responses. Thus, lower frequency of A allele in Iranian patients with T2D may contribute to the risk of CHD in these patients. It should be noted that mortality in diabetic patients is more due to cardiovascular complications.
Interindividual variations are involved in drug disposition and response. For example, a sufficient response is not observed in certain metformin users that is one of the most important predicaments about efficacy of the important antidiabetic drug., One reason for this poor response is genetic variability in metformin transporters. OCT3 such as OCT1 and OCT2 is an important transporter for metformin. There are few studies on the effect of genetic variants of OCT3 on the pharmacokinetics of metformin. Chen et al. showed that some different variants in OCT3 such as rs68187715 and rs8187717 affect metformin uptake in vitro. However, Tzvetkov et al. investigated the effect of some variants of OCT3 such as rs3123634 and rs2292334 on pharmacokinetics in healthy volunteers, but found no statistically significant relationship between these variants and clearance of metformin.
OCT3-564G>A is one of the variants in 3´-UTR region of OCT3 which affect the expression of this transporter. Li et al. showed that G >A displacement in the variant leads to binding miR-147 to mRNA resulting in the decreased expression of OCT3. Moreover, Nies et al. indicated that three variants in the noncoding region including OCT3-564G>A are associated with decreased expression of mRNA in OCT3 in the liver cells. In the present study, given the importance of OCT3-564G>A in OCT3 expression, we investigated whether this variant can affect the response to metformin in patients with T2D. As expected, the results showed that fasting glucose and HbA1c levels were significantly decreased after 3 months of treatment with metformin. When we analyzed these factors according to OCT3-564G>A genotypes, it was found that fasting glucose and HbA1c levels were significantly decreased in both GG and GA + AA groups after 3 months of metformin therapy. In addition, our study are consistent with results obtained by other authors , showed that metformin is effective in improving the lipid profile. The results analysis according to the genotypes of OCT3-564G>A showed that the lipid profile is improved after 3 months of treatment with metformin in both patients with GG and GA + AA genotypes.
Further analyses according to responders and nonresponders to metformin revealed that fasting glucose and HbA1c do not change significantly in responders and nonresponders with respect to the genotypes of study variant. In other words, it seems that OCT3-564G>A is not effective in the glycemic response to metformin. However, it should be noted that metformin can be transported by two other members of the OCT family, i.e., OCT1 and OCT2. This issue is more important in the major organ of metformin activity, i.e., the liver in which OCT1 is the main transporter. Hence, for a more comprehensive review of the role of OCT3 variants in response to metformin, these variants should be studied with functional variants of two other members particularly OCT1.
The relatively small sample size should be considered as a limitation of the study, and thus further research in different populations with a larger sample size is needed to clarify the role of the OCT3-564G>A variant in metformin response.
In general, considering the different consequences of the variations in drug transporters, the effect of OCT3 variations as a high-capacity transporter widely expressed in various tissues cannot be ignored.
| Conclusion|| |
The frequency of A allele (which may be a protective allele against CHD) in the Iranian diabetic patients was lower compared with Caucasian, Japanese, and Iranian healthy populations. Metformin is useful in improving the lipid profile, in addition to its effects in glycemic control, and these impacts are regardless of the variant OCT3-564G>A.
Financial support and sponsorship
This study received funding support from Mazandaran University of Medical Sciences.
Conflicts of interest
There are no conflicts of interest.
| Authors' Contribution|| |
AM contributed in the conception of the work, conducting the study, writing the article, revising the draft, approval of the final version of the manuscript, and agreed for all aspects of the work MG contributed in the conception of the work, study design, samples and data collection, drafting and revising the draft, approval of the final version of the manuscript, and agreed for all aspects of the work MBH contributed in the conception of the work, drafting and revising the draft, approval of the final version of the manuscript, and agreed for all aspects of the work SRH contributed in the conception of the work, drafting and revising the draft, approval of the final version of the manuscript, and agreed for all aspects of the work AA contributed in the conception of the work, data analysis, drafting and revising the draft, approval of the final version of the manuscript, and agreed for all aspects of the work SMN contributed in the conception of the work, drafting and revising the draft, approval of the final version of the manuscript, and agreed for all aspects of the work.
| References|| |
Zinman B, Harris SB, Neuman J, Gerstein HC, Retnakaran RR, Raboud J, et al.
Low-dose combination therapy with rosiglitazone and metformin to prevent type 2 diabetes mellitus (CANOE trial): A double-blind randomised controlled study. Lancet 2010;376:103-11.
Gao Y, Yoon KH, Chuang LM, Mohan V, Ning G, Shah S, et al.
Efficacy and safety of exenatide in patients of Asian descent with type 2 diabetes inadequately controlled with metformin or metformin and a sulphonylurea. Diabetes Res Clin Pract 2009;83:69-76.
Burcelin R, Serino M, Chabo C, Blasco-Baque V, Amar J. Gut microbiota and diabetes: From pathogenesis to therapeutic perspective. Acta Diabetol 2011;48:257-73.
Weber MB, Narayan KM. Preventing type 2 diabetes: Genes or lifestyle? Prim Care Diabetes 2008;2:65-6.
Nolan CJ, Damm P, Prentki M. Type 2 diabetes across generations: From pathophysiology to prevention and management. Lancet 2011;378:169-81.
Rizos CV, Elisaf MS. Metformin and cancer. Eur J Pharmacol 2013;705:96-108.
Kashi Z, Masoumi P, Mahrooz A, Hashemi-Soteh MB, Bahar A, Alizadeh A. The variant organic cation transporter 2 (OCT2)-T201M contribute to changes in insulin resistance in patients with type 2 diabetes treated with metformin. Diabetes Res Clin Pract 2015;108:78-83.
Zhou K, Pearson ER. OCT3 is a major determinant of metformin effect in skeletal muscle. Pharmacogenomics J 2011;12:708-9.
Gong L, Goswami S, Giacomini KM, Altman RB, Klein TE. Metformin pathways: Pharmacokinetics and pharmacodynamics. Pharmacogenet Genomics 2012;22:820-7.
Mahrooz A, Parsanasab H, Hashemi-Soteh MB, Kashi Z, Bahar A, Alizadeh A, et al.
The role of clinical response to metformin in patients newly diagnosed with type 2 diabetes: A monotherapy study. Clin Exp Med 2015;15:159-65.
Viollet B, Foretz M. Revisiting the mechanisms of metformin action in the liver. Ann Endocrinol (Paris) 2013;74:123-9.
Chen EC, Liang X, Yee SW, Geier EG, Stocker SL, Chen L, et al.
Targeted disruption of organic cation transporter 3 attenuates the pharmacologic response to metformin. Mol Pharmacol 2015;88:75-83.
Verhaagh S, Schweifer N, Barlow DP, Zwart R. Cloning of the mouse and human solute carrier 22a3 (Slc22a3/SLC22A3) identifies a conserved cluster of three organic cation transporters on mouse chromosome 17 and human 6q26-q27. Genomics 1999;55:209-18.
DeGorter MK, Kim RB. Hepatic drug transporters, old and new: Pharmacogenomics, drug response, and clinical relevance. Hepatology 2009;50:1014-6.
Bleasby K, Castle JC, Roberts CJ, Cheng C, Bailey WJ, Sina JF, et al.
Expression profiles of 50×enobiotic transporter genes in humans and pre-clinical species: A resource for investigations into drug disposition. Xenobiotica 2006;36:963-88.
Wieland A, Hayer-Zillgen M, Bönisch H, Brüss M. Analysis of the gene structure of the human (SLC22A3) and murine (Slc22a3) extraneuronal monoamine transporter. J Neural Transm (Vienna) 2000;107:1149-57.
Kerb R. Implications of genetic polymorphisms in drug transporters for pharmacotherapy. Cancer Lett 2006;234:4-33.
Nies AT, Koepsell H, Winter S, Burk O, Klein K, Kerb R, et al.
Expression of organic cation transporters OCT1 (SLC22A1) and OCT3 (SLC22A3) is affected by genetic factors and cholestasis in human liver. Hepatology 2009;50:1227-40.
Li L, He M, Zhou L, Miao X, Wu F, Huang S, et al.
A solute carrier family 22 member 3 variant rs3088442 G>A associated with coronary heart disease inhibits lipopolysaccharide-induced inflammatory response. J Biol Chem 2015;290:5328-40.
Friedewald WT, Levy RI, Fredrickson DS. Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin Chem 1972;18:499-502.
Hengen N, Lizer MH, Kidd RS. Evaluation of genetic variations in organic cationic transporter 3 in depressed and nondepressed subjects. ISRN Pharmacol 2011;2011:161740.
Aoyama N, Takahashi N, Kitaichi K, Ishihara R, Saito S, Maeno N, et al.
Association between gene polymorphisms of SLC22A3 and methamphetamine use disorder. Alcohol Clin Exp Res 2006;30:1644-9.
Bell DS. Drugs for cardiovascular risk reduction in the diabetic patient. Curr Diab Rep 2001;1:133-9.
Pacanowski MA, Hopley CW, Aquilante CL. Interindividual variability in oral antidiabetic drug disposition and response: The role of drug transporter polymorphisms. Expert Opin Drug Metab Toxicol 2008;4:529-44.
Takane H, Shikata E, Otsubo K, Higuchi S, Ieiri I. Polymorphism in human organic cation transporters and metformin action. Pharmacogenomics 2008;9:415-22.
Kashi Z, Mahrooz A, Kianmehr A, Alizadeh A. The role of metformin response in lipid metabolism in patients with recent-onset type 2 diabetes: HbA1c level as a criterion for designating patients as responders or nonresponders to metformin. PLoS One 2016;11:e0151543.
Chen L, Pawlikowski B, Schlessinger A, More SS, Stryke D, Johns SJ, et al.
Role of organic cation transporter 3 (SLC22A3) and its missense variants in the pharmacologic action of metformin. Pharmacogenet Genomics 2010;20:687-99.
Tzvetkov MV, Vormfelde SV, Balen D, Meineke I, Schmidt T, Sehrt D, et al.
The effects of genetic polymorphisms in the organic cation transporters OCT1, OCT2, and OCT3 on the renal clearance of metformin. Clin Pharmacol Ther 2009;86:299-306.
Bosi E. Metformin – The gold standard in type 2 diabetes: What does the evidence tell us? Diabetes Obes Metab 2009;11 Suppl 2:3-8.
Abdul-Ghani MA, Puckett C, Triplitt C, Maggs D, Adams J, Cersosimo E, et al.
Initial combination therapy with metformin, pioglitazone and exenatide is more effective than sequential add-on therapy in subjects with new-onset diabetes. Results from the efficacy and durability of initial combination therapy for type 2 diabetes (EDICT): A randomized trial. Diabetes Obes Metab 2015;17:268-75.
[Table 1], [Table 2], [Table 3], [Table 4]
|This article has been cited by|
||Clinical and genetic predictors of diabetes drug’s response
| ||Adriana Fodor,Angela Cozma,Ramona Suharoschi,Adela Sitar-Taut,Gabriela Roman |
| ||Drug Metabolism Reviews. 2019; : 1 |
|[Pubmed] | [DOI]|
||The Role of Histamine in the Pathophysiology of Asthma and the Clinical Efficacy of Antihistamines in Asthma Therapy
| ||Kohei Yamauchi,Masahito Ogasawara |
| ||International Journal of Molecular Sciences. 2019; 20(7): 1733 |
|[Pubmed] | [DOI]|
||OCT3 promoter haplotype is associated with metformin pharmacokinetics in Koreans
| ||Eun Young Kwon,Jae-Yong Chung,Hyo Jin Park,Bo Min Kim,Minsuk Kim,Ji Ha Choi |
| ||Scientific Reports. 2018; 8(1) |
|[Pubmed] | [DOI]|
||Physiological serum 25-hydroxyvitamin D concentrations are associated with improved thyroid function—observations from a community-based program
| ||Naghmeh Mirhosseini,Ludovic Brunel,Giovanna Muscogiuri,Samantha Kimball |
| ||Endocrine. 2017; |
|[Pubmed] | [DOI]|