Literature DB >> 23867834

Homeostatic model assessment indices in evaluation of insulin resistance and secretion in hemodialysis patients.

Stanisław Niemczyk1, Katarzyna Szamotulska, Kinga Giers, Mariusz Jasik, Zbigniew Bartoszewicz, Katarzyna Romejko-Ciepielewska, Ewa Paklerska, Małgorzata Gomółka, Joanna Matuszkiewicz-Rowińska.   

Abstract

BACKGROUND: Some previous observations suggest that insulin resistance and glucose metabolism disturbances are frequent complications of chronic kidney disease. However, there are no conclusive studies on other indices of the effectiveness of insulin action in end-stage renal disease (ESRD) patients, including chronically hemodialysed (HD) ones.
MATERIAL AND METHODS: The groups comprised 33 non-diabetic ESRD hemodialysed patients and 33 healthy controls matched for age, sex, and body mass index (BMI). In both groups, HOMA-%B, HOMA-%S, HOMA-IR indices, and DI were calculated using HOMA1 and HOMA2 as measures of insulin resistance. The indices were also assessed in subgroups divided according to BMI.
RESULTS: Mean fasting plasma glucose concentrations were lower in ESRD patients than in healthy persons (82.4±10.4 vs. 93.9±11.6, p=0.001). Fasting serum insulin concentrations were similar in both groups (median 6.8 vs. 6.0 mU/l, p=0.698). HOMA1-%B values were higher in ESRD patients than controls (median 137.1 vs. 81.6, p=0.002). HOMA1-%S (median 75.6 vs. 71.5) and HOMA1-IR (median 1.3 vs. 1.4) values were not significantly different (p=0.264 and p=0.189, respectively). DI1 levels were higher for HD patients than for healthy subjects (median 1.16 vs. 0.53, p<0.001). In subgroup analysis, all statistically significant differences were restricted mainly to persons with BMI <25 kg/m2. Similar results as for the HOMA1 model were obtained for HOMA2.
CONCLUSIONS: 1. HOMA beta-cell function is strongly correlated with HOMA insulin resistance in HD patients. 2. In non-diabetic ESRD hemodialysed patients, the HOMA indices and DI may be useful and important models in interpretation of glucose metabolism disturbances.

Entities:  

Mesh:

Substances:

Year:  2013        PMID: 23867834      PMCID: PMC3724569          DOI: 10.12659/MSM.883978

Source DB:  PubMed          Journal:  Med Sci Monit        ISSN: 1234-1010


Background

Insulin resistance (IR) is a frequent complication in the end-stage renal disease (ESRD), including hemodialysed patients (HD), and is a significant risk factor for cardiovascular diseases in this group of patients. Insulin resistance is already present in the early stages of chronic kidney disease (CKD), even in patients with a normal glomerular filtration rate (GFR) [1]. GFR in the first period of CKD depends on the severity of proteinuria [2]. Disorders of carbohydrate metabolism are also common in patients with ESRD. Some scientists have evaluated IR in patients in various stages of CKD and proved that greater insulin resistance is associated with lower renal function [3]. Renal replacement therapy, which corrects the disorders that affect IR, may reduce IR intensity [4,5]. Some studies on IR concerning patients with ESRD before renal replacement therapy confirm the existence of increased values of indices of IR [6,7]. However, the study of Bastürk and Unsal proved that there was no significant difference in GFR at the end of the study between patients who had or did not have increased IR. Furthermore, HOMA-IR was not significantly different in patients with or without renal failure [8]. There are also no conclusive data on insulin secretion in ESRD. The method of establishing glucose metabolism disturbances in ESRD patients is not clearly defined. The objective of the present study was to investigate the usefulness of homeostatic model assessment indices of insulin action (HOMA-%B, HOMA-%S, HOMA-IR) as well as disposition index (DI) calculated with 2 methods in ESRD patients without diabetes in comparison to healthy control subjects.

Material and Methods

Analytical procedures and data analysis

The study included 66 patients, aged 50 to 80 years old, hospitalized in the Nephrology Department at the Medical University of Warsaw. Group I comprised 33 non-diabetic ESRD subjects undergoing haemodialysis (25 men and 8 women) and group II included 33 healthy controls matched for age, sex, and body mass index (BMI). The groups were subdivided by BMI into the following subgroups: <25 (n=15 pairs), 25–30 (n=9 pairs), and >30 (n=9 pairs) kg/m2. The cause of ESRD included idiopathic chronic glomerulonephritis in non-diabetic individuals. Fasting plasma glucose concentrations were assayed by glucose hexokinase (Roche Cobas Integra 800 analyzer). Serum insulin concentrations were determined by IRMA (Immunotech). HOMA-IR, HOMA-%B, HOMA-%S, and disposition index were calculated by 2 methods: HOMA1 and HOMA2. For HOMA1, the following formulae were used. HOMA-IR2, HOMA-%B2 and HOMA-%S2 were derived using the HOMA calculator () and disposition index 2 results from multiplication of HOMA-%S2/100 by HOMA-%B2/100. The study was approved by the Ethics Committee at the Medical University of Warsaw on 20 June 2006 with the amendment made on the 21 October 2009. Informed consent was given by every patient.

Statistical methods

All statistical data are presented as means, standard deviations (SD), and medians. Wilcoxon exact test was used to compare matched groups. The correlation between metabolic variables was assessed by Spearman’s rho. A p value of <0.05 was considered statistically significant. The data were analysed using the program SPSS v.18.

Results

The results of fasting glycemia and fasting serum insulin are shown in Table 1 and 2, and the results of HOMA parameters and DI in the ESRD patients and control subjects are presented in Tables 3–5.
Table 1

Fasting plasma glucose concentration in non-diabetic ESRD patients undergoing haemodialysis and healthy control group.

Fasting plasma glucose (mg/dl)HDControl groupp-value*
BMI <25.0Mean ±SD80.9±9.693.9±8.50.011
Median78.095.0
BMI 25.0–30.0Mean ±SD77.6±8.792.3±15.60.078
Median79.089.0
BMI ≥30.0Mean ±SD88.8±11.295.4±12.60.359
Median92.093.0
TotalMean ±SD82.4±10.493.9±11.60.001
Median79.093.0

Wilcoxon exact test

Table 2

Fasting serum insulin concentration in non-diabetic ESRD patients undergoing haemodialysis and healthy control group.

Fasting plasma insulin (mU/l)HDControl groupp-value*
TotalMean ±SD9.0±11.68.4±5.90.698
Median6.86.0

Wilcoxon exact test

Table 3

Homeostatic model assessment of β-cell function (HOMA1 and 2-%B) in non-diabetic ESRD patients undergoing haemodialysis and healthy control group.

BMIHOMA1-βC (33 pairs)HOMA2-βC (29 pairs)
HDControl groupp-value*HDControl groupp-value*
<25.0Mean ±SD161.4±138.788.8±108.50.022114.0±42.382.5±49.00.042
Median132.956.0102.267.38
25.0–30.0Mean ±SD276.4±265.6112.0±70.00.027126.1±47.693.0±40.30.027
Median164.883.5113.881.2
≥30.0Mean ±SD202.4±188.7135.7±54.90.496104.9±30.5112.2±32.20.945
Median130.8109.799.3100.8
TotalMean ±SD203.9±193.1107.9±86.70.002115.2±40.694.0±42.70.013
Median137.181.6109.381.2

Wilcoxon exact test.

Table 5

Disposition index 1 (DI1) in non-diabetic ESRD patients undergoing haemodialysis and healthy control group.

BMIDI1 (33 pairs)DI2 (29 pairs)
HDControl groupp-value*HDControl groupp-value*
<25.0Mean ±SD1.37±0.880.57±0.270.0051.29±0.370.99±0.170.052
Median1.250.481.320.94
25.0–30.0Mean ±SD1.99±1.780.73±0.430.0981.42±0.441.08±0.460.164
Median1.160.631.331.13
≥30.0Mean ±SD0.91±0.740.58±0.320.3591.10±0.320.84±0.330.078
Median0.550.520.980.74
TotalMean ±SD1.41±1.200.62±0.33<0.0011.28±0.390.98±0.330.006
Median1.160.531.290.96

Wilcoxon exact test.

Fasting plasma glucose and serum insulin concentrations

Glucose concentrations were significantly lower in ESRD patients compared with healthy subjects (p=0.001). However, this difference was higher in persons with BMI <25 and 25–30 kg/m2 (p=0.011 and p=0.078, respectively) than in obese subjects (p=0.359) (Table 1). There was no significant difference between insulin concentrations in ESRD patients and healthy subjects (median 6.8 vs. 6.0 mU/l, p=0.698), irrespective of BMI (Table 2).

Homeostatic model assessment of β-cell function (HOMA1-%B)

The homeostatic model assessment of β-cell function (HOMA1-%B) values were significantly higher in ESRD patients compared with healthy subjects (median 137.1 vs. 81.6, p=0.002). In subgroups with BMI <25.0 and 25.0–30.0 kg/m2, HOMA1-%B values were significantly higher in HD patients compared with healthy subjects (p=0.022 and p=0.027, respectively). In the group with BMI >30 kg/m2 there was no significant difference between HD patients and controls (p=0.496) (Table 3). Similarly, the HOMA2-%B values were significantly higher in ESRD patients compared with healthy subjects (median 109.3 vs. 81.2, p=0.013).

Homeostatic model assessment of insulin sensitivity (HOMA1-%S) and insulin resistance (HOMA1-IR)

The values of homeostatic model assessment of insulin sensitivity HOMA1-%S (median 75.6 vs. 71.5, p=0.264) and insulin resistance HOMA1-IR (median 1.3 vs. 1.4, p=0.189) were not significantly different in ESRD patients compared with healthy controls (Table 4).
Table 4

Homeostatic model assessment insulin sensitivity (HOMA1-%S) and homeostatic model assessment insulin resistance (HOMA1-IR) in non-diabetic ESRD patients undergoing haemodialysis and healthy control group.

BMIHOMA1 (33 pairs)HOMA2 (29 pairs)
HDControl groupp-value*HDControl groupp-value*
HOMA-IRMean ±SD2.0±3.02.0±1.50.1891.0±0.41.1±0.80.559
Median1.31.40.90.8
HOMA-S%Mean ±SD91.5±66.372.9±40.80.264121.9±53.4121.6±59.90.189
Median75.671.5107.4126.4

Wilcoxon exact test.

No statistically significant difference was found for HOMA2-%S (p=0.189) and HOMA2-IR (p=0.559).

Homeostatic model assessment of disposition index

Disposition index 1 (DI1) (Figure 1) was higher for the whole HD group than for controls (median 1.16 vs. 0.53, p<0.001). However, after subgroup analysis, the difference was significant in the group with BMI <25 kg/m2 only (median 1.25 vs. 0.48 p=0.005) (Table 5).
Figure 1

Insulin sensitivity and secretion in HOMA 1 (A) and HOMA 2 (B) models in non-diabetic ESRD patients undergoing haemodialysis and healthy control group.

Disposition index 2 (DI2) was also significantly higher for the HD patients (median 1.29 vs. 0.96, p=0.006).

Discussion

IR significantly contributes to the development of carbohydrate metabolism disorders in many diseases, including ESRD patients. Clinically, it is characterized by normal serum insulin concentration associated with abnormal glucose response [9]. IR is associated with prevalent CKD and fast decrease in renal function in elderly patients, whereas co-existing metabolic syndrome predicts the risks of prevalent and incident CKD [10]. DeFronzo et al. claimed that in IR accompanying uremia, the suppressive effect of insulin on gluconeogenesis in the liver or stimulated glucose uptake by hepatocytes is normal [11]. Not all researchers confirm the increase of insulin resistance in patients with CRF. It is indicated in our study as well (Table 2) [12]. The gold standard in evaluating insulin resistance is believed to be the euglycemic clamp method described by DeFronzo et al. [11], which is considered the best technique for insulin resistance assessment because it provides a direct measurement of the whole-body sensitivity to insulin, particularly in skeletal muscle. This technique differentiates between peripheral and hepatic insulin resistance because of a direct and accurate measurement [13]. The HOMA-IR test evaluates hepatic rather than peripheral insulin resistance [13]. In CKD, insulin resistance exists mostly as a peripheral mechanism. Some researchers believe that HOMA-IR measurement cannot be the gold standard and is not an accurate method for evaluation of insulin resistance [11,13]. However, a 2000 study by Bonora et al. showed an excellent correlation between euglycemic hyperinsulinemic clamp and HOMA-IR measurements in patients with various degrees of glucose tolerance and insulin sensitivity [14,15]. Shoji et al. showed that HOMA-IR can be alternative technique to assess resistance to insulin in patients with and without renal failure [8,16]. The cause of IR and accompanying chronic renal failure seem to be multifactorial and very complex. The post-receptor signaling pathways of insulin seem to be essential [17]. IR is a derivative of disorders connected with the uremic environment; hence the influence of chronic inflammation, anemia, secondary hyperparathyroidism, and chronic acidosis is stressed. A significant role is attributed to the function of adipose tissue, especially visceral tissue, playing a more important part in secretion of many adipocytokines, which concentrations in renal failure are significantly increased [18]. The following factors determine the existence of IR in CFR patients: excess visceral fat, adipokine deregulation and accumulation, metabolic acidosis, oxidative stress, vitamin D deficiency, decreased physical activity, and accumulation of uremic toxins [19]. Insulin resistance is correlated with serum leptin concentrations and it is strongly correlated with the plasma leptin: adiponectin ratio (LAR) [20,21]. It is widely accepted that ESRD patients may show glucose intolerance due to IR, as the effect of their reduced peripheral sensitivity to the hypoglycaemic activity of insulin. Hemodialysis improves IR observed in uremic patients and, secondarily to this, glycemia could be significantly lowered in ESRD patients compared to healthy controls, as in our results (Table 1) [5]. In our study, there was no significant difference in fasting serum insulin concentrations between hemodialysed ESRD patients and healthy subjects (Table 2). Long-term consequences of increased IR include development of diabetes type 2, many vascular disorders, and malignancies. At the moment no single, commonly accepted test for establishing IR in clinical practice is available [22]. Relatively simple, increasingly common, non-invasive alternatives to the clamp technique have been proposed, such as homeostasis model assessment (HOMA) [23]. There is an inverse linear correlation between IR measured by euglycemic hyperinsulinemic clamp and log-transformed IR calculated by HOMA. IR measured by HOMA is considered an alternative method to assess IR in patients with renal failure [13,24]. The results of our study seem to confirm the usefulness of HOMA. The most popular, commonly accepted method of IR measurement is the HOMA-IR index [25]. There are 2 main methods of calculating HOMA-IR: from the formula and with the calculator (the Oxford method used to calculate HOMA2 is considered to be more accurate [26]). This also applies to other examined indices. HOMA-%S indicates insulin sensitivity and has various capabilities for assessment. This is the inverse of HOMA-IR multiplied by 100. The HOMA-%B index is a method for assessing β-cell function from basal glucose and insulin concentrations. The disposition index is an effect of insulin sensitivity and insulin secretion. It is generally constant for a patient; a change in value seems to be a very sensitive parameter of disturbances in glucose metabolism, as found in our study [27,28]. Regardless of the method of calculation, in our material, the DI in patients on dialysis is higher than in patients without renal failure (Table 5). HOMA analysis is important in both clinical and epidemiological studies. However, the HOMA-IR (or HOMA-%S) index is usually examined and HOMA-%B and DI are calculated less often. For example, the study by Heald et al. confirmed the usefulness of the HOMA-%B index in a general population of migrants from India [29]. Matthews et al. demonstrated the usefulness of this test for the determination of pancreatic beta-cell function both in diabetic and non-diabetic patients [30]. HOMA-%B has also been used to assess pancreatic function in many other populations [31-34]. It is also useful in predicting the function of beta-cells in time [35]. A correlation of 0.7 was demonstrated between HOMA-%B and metabolic clamp study in assessing the function of beta-cells in patients with diabetes and a correlation of 0.5 to 0.8 with the intravenous glucose tolerance test in non-diabetic and diabetic patients [30,33]. Hermans et al. demonstrated the usefulness of the beta-cell function index in the dynamic CIGMA test (CIGMA-%B) with continuous infusion of glucose [36]. CIGMA-%B provides better discrimination for distinguishing beta-cell function than other tests in patients with normal or impaired glucose tolerance and type 2 diabetes. The intravenous glucose tolerance test did not demonstrate any advantage in quantifying insulin secretion in these patients compared with the model-based HOMA and CIGMA tests. It has also been demonstrated that the areas under the curve of insulin concentrations, C-peptide, and glucose obtained in the course of the oral glucose tolerance test have greater usefulness for scientists in the evaluation of insulin secretion by pancreatic beta-cells than HOMA-%B [37,38]. The interpretation of HOMA-%B should be simultaneous with the evaluation of HOMA-%S [28]. In our study, HOMA-%B was higher in HD patients than in healthy subjects, regardless of the method of counting. This proves undisturbed pancreatic function in patients with ESRD. On the other hand, there was no difference in HOMA-IR and HOMA-%S between HD and the control group, which may suggest low sensitivity of the examined indices, the actual lack of a difference, or low usefulness of examined tests carried out on an empty stomach (Tables 3 and 4, Figure 1). The role of overstimulation of the pancreas to produce insulin in CRF is not obvious. Tuzcu et al. analyzed HOMA-%B in HD and continuous ambulatory peritoneal dialysis (CAPD) patients and found higher values in CAPD than HD or the control group, but similar in HD and the control group [39]. Shehab-Eldin et al. also did not confirm beta-cell dysfunction, despite the fact that the results of this study were similar to ours [36]. Pancreatic beta-cell function was also established in other studies, in non-diabetic and diabetic patients with normal and impaired renal function using HOMA-%B, plasma glucose, and insulin concentrations obtained at fasting or during 75-gram oral glucose tolerance test (insulin sensitivity index), insulinogenic index, first-phase insulin secretion index, and area under the curve of plasma insulin [12]. In non-diabetic patients, there was no difference in insulin secretion between patients with normal and impaired renal function. The homeostatic responsivity index HOMA-B, derived from basal measurements of insulin and glucose, is a relatively easy and common method of assessment of beta-cell function. The use of C-peptide instead of insulin in HOMA-B has been encouraged to avoid the confounding effect of hepatic insulin extraction. Although HOMA-B is widely used because of its simplicity, it has its limitations because it is assessed under non-stimulated conditions [28]. Dynamic tests seem to be necessary for a precise evaluation of disorders of carbohydrate metabolism in patients with CRF, as well as in elderly patients [41]. Attempts to determine the usefulness and accuracy of HOMA methods are important and were presented in this article. It is significant because these methods are safe, simple, and non-aggravating for the patient [13]. The results of our study can be helpful in choosing the method of assessment of insulin resistance in HD patients. Homeostatic model assessment beta-cell function (HOMA-%B) is strongly correlated with insulin resistance level (HOMA IR) in HD patients. In non-diabetic and end-stage renal disease patients undergoing hemodialysis, the homeostatic model assessment indices and disposition index may be useful and important models in interpretation of glucose metabolism disturbances.
  41 in total

1.  The homeostasis model in the San Antonio Heart Study.

Authors:  S M Haffner; H Miettinen; M P Stern
Journal:  Diabetes Care       Date:  1997-07       Impact factor: 19.112

2.  Correct homeostasis model assessment (HOMA) evaluation uses the computer program.

Authors:  J C Levy; D R Matthews; M P Hermans
Journal:  Diabetes Care       Date:  1998-12       Impact factor: 19.112

3.  Insulin resistance and hyperinsulinemia are already present in patients with incipient renal disease.

Authors:  D Fliser; G Pacini; R Engelleiter; A Kautzky-Willer; R Prager; E Franek; E Ritz
Journal:  Kidney Int       Date:  1998-05       Impact factor: 10.612

4.  A prospective analysis of the HOMA model. The Mexico City Diabetes Study.

Authors:  S M Haffner; E Kennedy; C Gonzalez; M P Stern; H Miettinen
Journal:  Diabetes Care       Date:  1996-10       Impact factor: 19.112

5.  Insulin resistance in uremia mediated by postbinding defects.

Authors:  D Smith; R A DeFronzo
Journal:  Kidney Int       Date:  1982-07       Impact factor: 10.612

6.  Insulin resistance in patients with chronic kidney disease.

Authors:  Shuzo Kobayashi; Kyoko Maesato; Hidekazu Moriya; Takayasu Ohtake; Toshio Ikeda
Journal:  Am J Kidney Dis       Date:  2005-02       Impact factor: 8.860

7.  Comparison of tests of beta-cell function across a range of glucose tolerance from normal to diabetes.

Authors:  M P Hermans; J C Levy; R J Morris; R C Turner
Journal:  Diabetes       Date:  1999-09       Impact factor: 9.461

8.  Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man.

Authors:  D R Matthews; J P Hosker; A S Rudenski; B A Naylor; D F Treacher; R C Turner
Journal:  Diabetologia       Date:  1985-07       Impact factor: 10.122

9.  High prevalence of abnormal glucose tolerance and metabolic disturbances in first degree relatives of NIDDM patients. A study in Catalonia, a mediterranean community.

Authors:  A Costa; M Rios; R Casamitjana; R Gomis; I Conget
Journal:  Diabetes Res Clin Pract       Date:  1998-09       Impact factor: 5.602

10.  Serum leptin concentrations correlate to plasma insulin concentrations independent of body fat content in chronic renal failure.

Authors:  P Stenvinkel; O Heimbürger; F Lönnqvist
Journal:  Nephrol Dial Transplant       Date:  1997-07       Impact factor: 5.992

View more
  8 in total

1.  Kidney function, β-cell function and glucose tolerance in older men.

Authors:  Ting Jia; Ulf Risérus; Hong Xu; Bengt Lindholm; Johan Ärnlöv; Per Sjögren; Tommy Cederholm; Tobias E Larsson; Talat Alp Ikizler; Juan J Carrero
Journal:  J Clin Endocrinol Metab       Date:  2014-11-27       Impact factor: 5.958

2.  Visceral adipose tissue is associated with insulin resistance in hemodialyzed patients.

Authors:  Kinga Giers; Stanisław Niemczyk; Katarzyna Szamotulska; Katarzyna Romejko-Ciepielewska; Ewa Paklerska; Zbigniew Bartoszewicz; Ryszard Pacho; Mariusz Jasik; Joanna Matuszkiewicz-Rowińska
Journal:  Med Sci Monit       Date:  2015-02-20

3.  Preoperative oral carbohydrate improved postoperative insulin resistance in rats through the PI3K/AKT/mTOR pathway.

Authors:  Zhiguo Wang; Yiqing Liu; Qi Li; Canping Ruan; Bin Wu; Qiang Wang; Zhiqian Hu; Huanlong Qin
Journal:  Med Sci Monit       Date:  2015-01-01

4.  A Prospective Study of Renal Transplant Recipients: A Fall in Insulin Secretion Underpins Dysglycemia After Renal Transplantation.

Authors:  David Langsford; Varuni Obeyesekere; Sara Vogrin; Jessie Teng; Richard J MacIsaac; Glenn Ward; Frank Alford; Karen M Dwyer
Journal:  Transplant Direct       Date:  2016-10-07

5.  A Pilot Study on the Metabolic Impact of Mediterranean Diet in Type 2 Diabetes: Is Gut Microbiota the Key?

Authors:  Shámila Ismael; Marta P Silvestre; Miguel Vasques; João R Araújo; Juliana Morais; Maria Inês Duarte; Diogo Pestana; Ana Faria; José B Pereira-Leal; Joana Vaz; Pedro Ribeiro; Diana Teixeira; Cláudia Marques; Conceição Calhau
Journal:  Nutrients       Date:  2021-04-08       Impact factor: 5.717

6.  Chronic activation of cardiac Atg-5 and pancreatic Atg-7 by intermittent fasting alleviates acute myocardial infarction in old rats.

Authors:  Sahar Mohamed El Agaty; Noha A Nassef; Doaa A Abou-Bakr; Aya A Hanafy
Journal:  Egypt Heart J       Date:  2022-04-13

7.  Effect of Date Fruit Consumption on the Glycemic Control of Patients with Type 2 Diabetes: A Randomized Clinical Trial.

Authors:  Alexandra E Butler; Jenan Obaid; Pearl Wasif; Jean V Varghese; Rawan Abdulrahman; Dalal Alromaihi; Stephen L Atkin; Naji Alamuddin
Journal:  Nutrients       Date:  2022-08-25       Impact factor: 6.706

8.  One Year of Yoga Training Alters Ghrelin Axis in Centrally Obese Adults With Metabolic Syndrome.

Authors:  Angus P Yu; Felix N Ugwu; Bjorn T Tam; Paul H Lee; Christopher W Lai; Cesar S C Wong; Wendy W Lam; Sinead Sheridan; Parco M Siu
Journal:  Front Physiol       Date:  2018-09-20       Impact factor: 4.566

  8 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.