Skip to main content

ORIGINAL RESEARCH article

Front. Endocrinol., 25 July 2024
Sec. Clinical Diabetes

Sex-specific association of circulating Isthmin-1 with isolated post-challenge hyperglycemia

Jiahua Fan,,&#x;Jiahua Fan1,2,3†Jialin He,&#x;Jialin He2,3†Jiangyuan Zhu,Jiangyuan Zhu2,3Jialu Yang,Jialu Yang2,3Jingmeng Ju,Jingmeng Ju2,3Jingyi Huang,Jingyi Huang2,3Zhihao HuangZhihao Huang2Zhuoyu Zhang,Zhuoyu Zhang2,3Wenkang LiWenkang Li2Min Xia,Min Xia2,3Yan Liu,*Yan Liu2,3*
  • 1Department of Clinical Nutrition, Guangzhou Chest Hospital, Guangzhou, Guangdong, China
  • 2Guangdong Provincial Key Laboratory of Food, Nutrition and Health, Guangzhou, Guangdong, China
  • 3Department of Nutrition, School of Public Health, Sun Yat-sen University, Guangzhou, Guangdong, China

Introduction: To explore the distribution of Isthmin-1 (ISM1) level and its association with isolated post-challenge hyperglycemia (IPH).

Methods: A total of 522 participants without a history of diabetes were invited to attend a standard 75g 2-h oral glucose tolerance test (OGTT), and 71 subjects were further invited for a 3-h oral minimal model test. Insulin sensitivity and β-cell function were evaluated using both HOMA and estimated from OGTT. Circulating ISM1 levels were determined by a commercially available ELISA kit.

Results: A total of 76 (14.6%) participants were diagnosed as IPH, accounting for 61.3% of the newly diagnosed diabetes. ISM1 levels were significantly higher in men than in women (1.74 ng/mL versus 0.88 ng/mL). The inverse correlation between ISM1 and β-cell function and IPH was only significant in men. After multivariate adjustment, per unit increment in ISM1 was associated with 0.68-fold (95% CI: 0.49-0.90) reduced odds ratio (OR) of IPH in men. Compared to men with the lowest ISM1 levels, the adjusted OR of IPH with the highest ISM1 levels decreased by 73% (95% CI: 0.11-0.61). Moreover, incorporation of ISM1 into the New Chinese Diabetes Risk Score (NCDRS) model yielded a substantial improvement in net reclassification improvement of 58% (95% CI: 27%-89%) and integrated discrimination improvement of 6.4% (95% CI: 2.7%-10.2%) for IPH.

Conclusions: ISM1 was significantly and independently associated with IPH, and serves as a feasible biomarker for the early identification of men with high risk of IPH.

GRAPHICAL ABSTRACT
www.frontiersin.org

Graphical Abstract ISM1 was significantly and independently associated with IPH, and might be serves as a feasible biomarker for the early identification of men with high risk of IPH.

Highlights

● Isthmin-1 (ISM1) is a newly identified adipokine which promotes glucose uptake and improves glucose tolerance in murine models. However, its distribution and association with glycemic control in general population remain largely unknown.

● To explore the sex-specific distribution of ISM1 and its association with post-challenge hyperglycemia.

● Circulating ISM1 is substantially higher in men than women. Significant association of ISM1 with post-challenge hyperglycemia exists only in men.

● Circulating ISM1 serves as a potential biomarker for post-challenge hyperglycemia in men.

Introduction

Isolated post-challenge hyperglycemia (IPH), typically defined as 2-h glucose level higher than 11.1 mol/L with normal fasting glucose, is an established feature of diabetes and has been epidemiologically associated with worsening cardiometabolic health (13) and increased risk of mortality and morbidity (46). It is estimated that 46.6% of adults with undiagnosed diabetes have IPH in China (7), and the prevalence of IPH has been reported to account for up to half of the newly diagnosed diabetes in Asian populations (8). Timely pharmacological or lifestyle interventions to decrease the long-term risk in individuals with IPH is therefore mandatory. To date, however, IPH can only be detected by 2-h OGTT, which is time-consuming, inconvenient and not generally included in regular physical examinations, thus leaving a large proportion of patients with IPH undiagnosed and untreated. Clearly, considering the high and increasing prevalence of diabetes worldwide (9), identifying novel biomarkers for easy and early identification of subjects at high risk of IPH is of significance.

Isthmin-1 (ISM1), was first identified as a gene expressed in the Xenopus midbrain-hindbrain organizer called isthmus, with a role in early brain development (10). Recently, it was found to be a secreted polypeptide hormone that regulates adipose tissue glucose uptake (11), which plays an essential role in the maintenance of post-challenge glucose homeostasis (11). Ablation of ISM1 impaired both basal and insulin-stimulated glucose uptake by adipocytes, whereas, therapeutic administration of recombinant ISM1 improved diabetes in diet-induced obese mice (11). Moreover, despite a small sample size, ISM1 was proven to be detectable in the circulation and was reported to be lower in patients with type 2 diabetes than in those with normoglycemic control (12). However, as a newly identified adipokine, clinical investigation of circulating ISM1 in the general population is sparse and its association with IPH remains largely unknown.

In this study, we assessed the distribution of circulating ISM1 in a deeply phenotyped cohort of community residents and determined its association with IPH.

Research design and methods

Study population

This study was approved by the Ethics Committee of the School of Public Health, Sun Yat-sen University (2019-127) and was conducted in accordance with the principles of the Declaration of Helsinki. Written informed consent was obtained from each individual. Local residents, who had lived in Dongguan City, Guangdong, China for over 5 years were invited for diabetes screening through flyers and advertisements. The detailed inclusion and exclusion criteria were as follows:

Inclusion criteria: 1) aged 30-65 years old; 2) stable body weight (<10% of body weight change) for three months before the study; and 3) absence of diagnosed diabetes.

Exclusion criteria: 1) history of malignant tumor, severe liver or kidney disease, chronic hepatitis, thyroid disorders, autoimmune diseases, cardiovascular or cerebrovascular diseases; and 2) women during pregnancy or lactation.

A total of 535 participants were consecutively recruited, and 13 participants with outlier ISM1, defined as an extreme value of the upper 2.5% of the distribution (13) were excluded. Finally, a total of 522 participants were included in the current analysis.

Collection and assessment of covariates

Demographic characteristics, medication, and lifestyles were collected through structured questionnaires by trained staffs. Current smokers were defined as those who had smoked at least one cigarette per day for more than 6 months (14). Current drinking was defined as taking any alcoholic beverage at least once a week in the past 6 months (14). Physical activity was assessed by the short form of the International Physical Activity Questionnaire (IPAQ). Participants were categorized into low, moderate and high levels according to the guidelines for data processing and analysis of the IPAQ and subjects with moderate or higher intensity activity was defined as having an exercise behavior (14).

Physical examinations were carried out by trained personnel. Blood pressure was measured on the right upper arm in the sitting position after at least 10-15 minutes of rest using a validated digital automatic analyzer (Omron HEM-7136). Hypertension was defined as systolic blood pressure (SBP) ≥140 mmHg, diastolic blood pressure (DBP) ≥90 mmHg, self-reported hypertension previously diagnosed by physician or current use of antihypertensive medications (15, 16). Overweight was defined as BMI ≥25 kg/m2, and central obesity was defined as waist circumference ≥80 cm for women and 94 cm for men (17). Metabolic syndrome was defined according to the National Cholesterol Education Programme Adult Treatment Panel III (NCEP ATP III) criteria (18).

Biochemical measurements

Venous blood was taken after overnight fasting for 10-12 hours and subjected to laboratory analysis at local clinical laboratory. Glucose (fasting, 1-h and 2-hpost 75g glucose challenge), triglyceride (TG), total cholesterol (TC), HDL cholesterol (HDL-c), LDL cholesterol (LDL-c) and serum creatinine (Cr) concentrations were measured using an autoanalyzer (BS-800 Biochemical Autoanalyzer; Mindray; China). Estimate glomerular filtration rate (eGFR) was calculated using the Chronic Kidney Disease Epidemiology Collaboration equation (CKD-EPI) (19). Glycated hemoglobin (HbA1c) was measured by a designated high-performance liquid chromatography method (HA-8180, Arkray, Shiga, Japan). Insulin and C-peptide levels were determined with commercial ELISA kits (Mercodia, Uppsala, Sweden).

Ascertainment of glucose metabolism

All eligible participants were asked to fast overnight and were required to attend the 2-h standard oral glucose tolerance test (OGTT), with blood sampling at 0, 60 min and 120 min after oral loading of 75 g glucose. Among them, 71 subjects (38 men and 33 women) further participated in the 3-h oral minimal model tests with blood sampling at 0, 10, 20, 30, 60, 90, 120, and 180 min after a 75-g oral glucose challenge (20, 21). According to the American Diabetes Association (ADA) thresholds of glycemic categories (22), newly diagnosed diabetes was defined as fasting blood glucose (FBG) ≥ 7.0 mmol/L (126 mg/dL) or 2-h glucose ≥ 11.1 mmol/L (200 mg/dL); post-challenge hyperglycemia (PH) was defined as 2-h glucose ≥11.1 mmol/L (200 mg/dL); IPH was defined as FPG < 7.0 mmol/L (<126 mg/dL) and 2-h glucose ≥11.1 mmol/L (200 mg/dL); prediabetes was defined FBG of 5.6-6.9 mmol/L (100-125 mg/dL) or 2-h glucose of 7.8-11.0 mmol/L (140-199 mg/dL); impaired glucose tolerance (IGT) was defined as 2-h glucose of 7.8-11.0 mmol/L (140-199 mg/dL); isolated impaired glucose tolerance (iIGT) was defined as FBG < 5.6 mmol/L (< 100 mg/dL) and 2-h glucose of 7.8-11.0 mmol/L (140-199 mg/dL). Detailed definition and classification criteria of glucose status based on the 2-h OGTT are shown in Figure 1.

Figure 1
www.frontiersin.org

Figure 1 Diagram showing the definition of glycemic status. The definition and diagnostic criteria of different glycemic status was based on 2-h OGTT according to the guideline from American Diabetes Association.

Determination of ISM1

Fasting serum ISM1 levels were determined by a commercial ELISA kit (Biomatlk, Canada) according to the manufacturer’s instructions and compared with purified human ISM1 standards. The mean intra-assay and inter-assay coefficients of variation were 2.9-8.5% and 4.3-14.9%, respectively.

Measurements of insulin sensitivity and β-cell function

The area under the curve (AUC) of glucose (AUCglu) and insulin (AUCins) over the 2 hours was calculated with the trapezoidal method (23). Insulin sensitivity was estimated by the following indices: 1) the homeostatic model assessment of insulin resistance (HOMA-IR) = FBG (mmol/L) × fasting insulin (mU/L)/22.5; 2) TyG index = log [TG (mg/dL) × FBG (mg/dL)/2] (24) and 3) Matsuda index = 10000/√[Glu0 (mg/dL) × Ins0 ((mU/L)] × [Glumean (mg/dL) × Insmean ((mU/L)] (25). β-Cell function was evaluated by the following 5 indices: 1) HOMA-β = fasting insulin (mU/L) × 20/[FBG (mmol/L) - 3.5]; 2) Stumvoll first phase insulin secretion = 2,032 + 4.681 × Ins0 (pmmol/L) - 135.0 × Glu120 (mmol/L) + 0.995 × Ins120 + 27.99 × BMI - 269.1 × Glu0 (mmol/L) (26); 3) Stumvoll second phase insulin secretion = 277 + 0.800 × Ins0 (pmmol/L) - 42.79 × Glu120 (mmol/L) + 0.321 × Ins120 (pmmol/L) + 5.338 × BMI; 4) The insulinogenic index (IGI60) = [Ins60 (mU/L) - Ins0 (mU/L)]/[Glu60 (mmol/L) - Glu0 (mmol/L)] (27); and 5) The oral disposition index (DI) was evaluated as AUCins/AUCglu multiple by Matsuda index (27). Ins0, Ins60, Ins120, Glu0, Glu60 and Glu120 were insulin and glucose during 2-h OGTT at 0, 60 and 120 min, respectively, and the units of insulin and glucose were shown in the above formula. Glumean and Insmean referred to the mean values of glucose and insulin at 0 and 120 min in 2-h OGTT.

Oral minimal model-derived measurements

During the 3-h OGTT, glucose, insulin and C-peptide concentrations were measured in mmol/L, mU/L and nmol/L, respectively. The insulinogenic index (IGI) = (Ins30 - Ins0)/(Glu30 - Glu0) and C-peptide index (CPI) = (Cpep30 - Cpep0)/(Glu30 - Glu0). Oral minimal model-derived measurements, including the insulin sensitivity index (SI), dynamic insulin secretion (PhiD), static insulin secretion (PhiS) and total insulin secretion (PhiT), were calculated. Insulin sensitivity was estimated by the insulin/glucose minimal model, insulin secretion measures were derived from the glucose/C-peptide minimal model with SAAM II 2.3 software (The Epsilon Group, Charlottesville, VA), and the glucose derivatives were determined by MATLAB R2016a software (MathWorks, Natick, MA), as previously described (28).

Statistical analysis

The normal distribution assumption was tested with the Shapiro-Wilk test. Data are presented as the mean ± SD for normally distributed variables and median with interquartile range for non-normally distributed parameters. Categorical variables were compared by χ2 test, and continuous variables were compared with independent Student’s t test or the Mann-Whitney U test as appropriate. To assess sex differences in ISM1, a least-squares regression model adjusted for age, history of hypertension, BMI, waist circumference, TG, HDL-c, eGFR and HOMA-IR was performed. In addition, differences in ISM1 levels according to different metabolic statuses were also determined by the least-squares regression model adjusted for age, and the box plots were based on the residuals from ordinary regression models of each response on age. Associations of ISM1 with various metabolic parameters, and indices related to insulin sensitivity and β-cell function were determined by age-adjusted partial Spearman correlation analysis. Linear regression models were performed to analyze the association of ISM1 with oral minimal model-derived measurements.

To evaluate the independent association of ISM1 with IPH in men and women separately, both linear and non-linear associations of ISM1 and glycemic status were first examined with restricted cubic spline (RCS) with 3 knots at the 10th, 50th and 90th centiles fitted in the multivariable logistic regression model. Age, history of hypertension, TG, HDL-c, eGFR, waist circumference and HOMA-IR were included in the RCS model as covariates. Furthermore, as a categorical variable, the independent associations of ISM1 with odds ratios (ORs) of IPH in men and women were determined with multivariable logistic regression model adjusted for potential confounders as specified below. Model 1: adjusted for age, history of hypertension, TG, HDL-c and eGFR. Model 2: Model 1 plus waist circumference. Model 3: Model 2 plus HOMA-IR. Moreover, to assess the additive value of ISM1 to the New Chinese Diabetes Risk Score (NCDRS), an established model for screening undiagnosed diabetes in Chinese individuals (29), the C-statistic, continuous net reclassification improvement (NRI) (30) and the integrated discrimination improvement (IDI) were employed. The likelihood ratio test was used to compare these models, and the calibration of these models was evaluated by the Hosmer-Lemeshow statistic. The confidence intervals (CIs) for the NRI values were computed using 1000 bootstrap samples.

Additionally, to enhance the robustness of our findings, several sensitivity analyses were performed. First, we restricted the analysis to participants who had complete data on 2-h OGTT tests and HbA1c results (n=428), where IPH was determined on both the OGTT test and HbA1c following the ADA guidelines (22). Second, ISM1 outliers were included in the analysis. Third, we expanded subjects with IPH into subjects with PH, defined as individuals with 2-h post-challenge hyperglycemia regardless of the level of fasting glucose. Finally, the E-value methodology was employed to assess the risk of residual confounding (31). We computed an E-value for the primary analysis using an online platform (32). Two-sided P values <0.05 were considered statistically significant and all statistical analyses were performed with R software, version 4.1.2.

Results

Sex-specific distribution of ISM1

The demographic and clinical characteristics of the study participants stratified by sex are presented in Table 1. Overall, the mean age was 47 years, and 48.7% were women. A total of 76 (14.6%) individuals were classified as IPH, accounting for 61.3% of the newly diagnosed diabetes in our cohort. Although there was no obvious difference in age, men were more likely to have unhealthy lifestyles, such as smoking and drinking, a higher prevalence of metabolic disorders, including diabetes, hypertension, dyslipidemia and metabolic syndrome, and worse metabolic profiles in terms of waist circumference, blood pressure and eGFR (P <0.05). Moreover, compared to women, both insulin sensitivity and β-cell function were markedly worse in men (Table 1).

Table 1
www.frontiersin.org

Table 1 Participant characteristics stratified by sex.

Of note, there was a significant difference in the distribution of circulating ISM1 by sex, with a median of 1.74 (1.02-2.63) ng/mL in men and 0.88 (0.22-1.72) ng/mL in women (P <0.001). The higher level of ISM1 in men remained significant even after adjustment for age, history of hypertension, TG, HDL-c, eGFR, BMI, waist circumference and HOMA-IR (P <0.05, Figure 2).

Figure 2
www.frontiersin.org

Figure 2 Distribution of Isthmin-1 according to sex. The P value indicated comparison between men and women after adjustment for age, history of hypertension, TG, HDL-c, eGFR, BMI, waist circumference and HOMA-IR. There was a significant difference in the distribution of circulating ISM1 by sex, with a median of 1.74 (1.02-2.63) ng/mL in men and 0.88 (0.22-1.72) ng/mL in women (P <0.001).

Sex-specific association of ISM1 with metabolic parameters, insulin sensitivity and β-cell function

In women, no significant difference in metabolic parameters, except for slightly decreased eGFR and improved HOMA-β was found with increasing tertiles of ISM1 (P >0.05, Supplementary Table S1). In contrast, in men, stumvoll 1 and stumvoll 2, two widely used measurements for insulin secretion capacity (26), were substantially impaired in men in the lowest tertile of ISM1. Consistently, the prevalence of diabetes, PH and IPH were also found to be remarkably higher in men with the lowest tertile of ISM1 expression (P <0.05, Supplementary Table S2). In a similar fashion, no obvious difference in ISM1 levels between females with different glycemic statuses or metabolic comorbidities was found (all P values >0.05, Figure 3A), while significantly lower ISM1s were observed in males with diabetes, PH and IPH (P <0.05, Figure 3B).

Figure 3
www.frontiersin.org

Figure 3 Distribution of ISM1 levels according to metabolic status in men and women. The box plots were based on the residuals from ordinary regression models of each response on age in (A) women and (B) men, respectively, according to different metabolic statuses. The P values were for comparisons between different metabolic statuses (yes and no) and adjusted for age. PH, post-challenge hyperglycemia; IPH, isolated post-challenge hyperglycemia; IGT, impaired glucose tolerance; iIGT, isolated impaired glucose tolerance.

Moreover, a distinct correlation between ISM1 and various metabolic parameters in men and women was found. Specifically, no significant correlation between ISM1 and any clinical features or metabolic indices related to glycemic control was found in women (all P >0.05, Figure 4A). On the contrary, ISM1 demonstrated a positive correlation with waist circumference and indices related to β-cell function, including Stumvoll 2 and disposition index (DI), in men after adjustment for age (P <0.05, Figure 4B).

Figure 4
www.frontiersin.org

Figure 4 Association of ISM1 with metabolic parameters, insulin sensitivity and β-cell function. Heatmap showing the age-adjusted Spearman correlation coefficient between ISM1 and clinical parameters in (A) women and (B) men, respectively. Insulin sensitivity and β-cell function indices were based on the 2-h OGTT test. The symbol × indicates insignificant correlation. All P values were adjusted for age. WC, waist circumference.

Notably, significant associations between ISM1 and β-cell function, in terms of IGI, PhiS and PhiT (oral minimal model-derived measurements), were also found in the subgroup of men who attended the 3 h-OGTT test (Figure 5A), whereas, no significant associations were observed between ISM1 and parameters related to insulin sensitivity in men and women (Figure 5B).

Figure 5
www.frontiersin.org

Figure 5 Association of ISM1 with oral minimal model-derived measurements. Linear regression models analyze the association of ISM1 with parameters related to insulin sensitivity and β-cell function in women (A) and men (B) with 3-h oral minimal model tests (38 men and 33 women). All P values were adjusted for age.

Independent association of circulating ISM1 with IPH in men

We then employed logistic regression models with cubic natural spline analyses to explore the association of serum ISM1 levels with ORs of IPH in men and women separately. In women, no association of ISM1 with the risk of IPH was found (Figures 6A, C), regardless of the menopause status (Supplementary Figure S1). However, in men, an increasing level of ISM1 was found to be significantly associated with a reduced risk of IPH (P=0.006, Figure 6B) after correcting for age, history of hypertension, waist circumference, TG, HDL-c, eGFR and HOMA-IR. Furthermore, when categorized into three groups according to the tertiles of ISM1 in our cohort, per unit increment in ISM1 was associated with a 32% (OR=0.68, 95% CI: 0.49-0.90) reduced risk of IPH in men. Compared to males in the lowest tertile of ISM1, those in the highest tertile of ISM1 experienced a greater than 70% (OR=0.26, 95% CI: 0.11-0.58) reduction in the risk of IPH after extensive adjustment for potential confounders (both P<0.05, Figure 6C).

Figure 6
www.frontiersin.org

Figure 6 Isthmin-1 and Odds ratios of isolated post-challenge hyperglycemia. The association of serum ISM1 levels (as a continuous variable) with ORs of IPH in (A) women and (B) men, respectively. ORs (solid yellow and blue lines) and 95% CIs (light yellow and light blue areas) were derived from multivariable logistic regression using cubic natural spline. The P values were adjusted for age, history of hypertension, waist circumference, TG, HDL-c, eGFR and HOMA-IR. (C) Association of ISM1 (as a categorical variable) and the ORs of IPH in women (up) men (down). Model 1: adjusted for age, history of hypertension, TG, HDL-c and eGFR. Model 2: Model 1 plus waist circumference. Model 3: Model 2 plus HOMA-IR. IPH, isolated post-challenge hyperglycemia.

Of note, the independent association of ISM1 with IPH in men remained robust in several sensitivity analyses, for instance, restricting the analysis to subjects with complete data of both 2 h-OGTT and HbA1c (n=428, 211 women and 217 men, Supplementary Figure S2) or expanding the analysis to the samples with ISM1 outliers (n=535, Supplementary Figure S3). Moreover, when using a broader definition of 2 h post-challenge hyperglycemia without considering the level of fasting glucose, the inverse association of ISM1 and risk for PH remained significant only in men (n=522, Supplementary Figure S4). Additionally, considering that most observed relationships between lifestyle and genetic factors and type 2 diabetes rarely exceeded 6.87 (the E-value for the adjusted OR of IPH with the highest tertile of ISM1) (9, 33), the E-value analysis suggested that no residual confounding factors remained in our study (Supplementary Table S3). Taken together, the above findings suggest that circulating ISM1 serves as an independent biomarker for IPH in men.

Incremental discriminative value of ISM1 for IPH in men

To further examine the clinical relevance of ISM1, we then evaluated the additive value of ISM1 for the identification of diabetes and IPH on top of NCDRS, a widely used tool for the screening of undiagnosed diabetes in Chinese individuals (29, 34). Incorporation of ISM1 led to a 6% increase in the global performance for IPH differentiation, a modest and only marginally significant improvement (P=0.056, Table 2). Regarding the reclassification potential, which conveys more clinically important information (35), additional inclusion of ISM1 to the NCDRS resulted in 14% and 44% of the population being correctly reclassified as cases and non-cases, respectively, thus yielding an NRI of 58% (95% CI: 27%-89%, P<0.001) and IDI of 6.4% (95% CI: 2.7%-10.2%, P=0.001) for the total population (Table 2). Of note, compared to traditional clinical indices, such as fasting glucose and HbA1c, the incorporation of ISM1 into the NCDRS demonstrated a significantly better discrimination ability for IPH (Table 2). Moreover, such a substantial improvement in risk discrimination and reclassification was not materially changed in all sensitivity analyses, as evidenced by an increase in NRI from 36% to 65% and IDI ranging from 2.7% to 6.2% (all P<0.05,Supplementary Tables S4–S6). All models had acceptable calibration (Hosmer-Lemeshow statistic χ2 range from 4.21 to 13.36, all P >0.10).

Table 2
www.frontiersin.org

Table 2 Additive differentiation and reclassification value of ISM1 for isolated post-challenge hyperglycemia on top of classical model.

Conclusions

In this deeply phenotyped cohort of medication-naïve community residents, we demonstrated an obvious sexual dimorphism in the distribution of circulating ISM1 and its association with β-cell function. Elevated ISM1 was independently associated with a decreased risk of IPH exclusively in men, but not in women. Previous studies have reported that ISM1 levels improved glucose tolerance and were associated with a lower risk of diabetes (11, 12). Conversely, other investigations found that serum ISM1 was significantly increased in patients with T2DM and correlated with the severity of albuminuria (36, 37). Such an ambiguous correlation might be mainly due to a limited sample size, methodology, and population characteristics, especially an incomplete adjustment for potential confounders. By selecting subjects without diagnosed diabetes and medication-naive subjects, extensively adjusting for potential confounders, and stratifying potential risk factors that may affect the ISM1-diabetes relationship, such as sex, we demonstrated a sexual difference of ISM1 with a decreased risk of IPH. Moreover, the incorporation of ISM1 into classical models for IPH screening yielded a significant improvement in the reclassification of Chinese subjects with IPH.

As an important subtype of diabetes, individuals with IPH have been reported to be at an increased risk for both incident cardiovascular disease and all-cause mortality compared to those with normal glucose levels (38, 39). Moreover, they also have an equal risk as those who have previously been diagnosed with diabetes mellitus (40). Additionally, incident diabetes and prediabetes in the elderly always manifest as post-challenge rather than fasting hyperglycemia (8). In our study, IPH accounted for up to 61.3% of newly diagnosed diabetes cases. Concordant with this result, the high prevalence of IPH of approximately 50-60% in subjects with newly diagnosed diabetes has also been described in previous population studies in Europe, Asia and China (7, 8, 41). Despite such a high burden, 2 h-OGTT remains an indispensable tool for the identification of IPH in current clinical practice, leaving an unmet need for novel biomarkers of IPH. In this regard, our findings that ISM1 is independently associated with β-cell function and serves as a sensitive biomarker for IPH in men are of major clinical importance. After adjustment for age, history of hypertension, waist circumference, TG, HDL-c, eGFR, and HOMA-IR, decreased ISM1 remained an independent risk factor for IPH in men. Notably, such an inverse association between circulating ISM1 and increased 2-h glucose persisted regardless of the fasting glucose level. Moreover, incorporation of ISM1 into the classical model for diabetes screening led to a net reclassification of 58% for the total population in our study, which means that approximately three in five subjects with undiagnosed diabetes could be identified by additional inclusion of ISM1, and appropriate interventions can thus be implemented in time to avoid severe and potentially irreversible complications of their unmanaged hyperglycemia.

As a newly identified adipokine, the circulating level of ISM1 and its distribution characteristics have been poorly investigated. Overall, the median level of ISM1 was 1.3 ng/L in our study, slightly lower than the level of 3.72 ng/mL in subjects with a significantly higher BMI (27 in the previous study vs. 25 in our cohort) (12), but comparable to the level of approximately 0.9-1.1 ng/mL found in individuals with diabetes and a similar BMI range to ours (36). Although it was reported to be positively associated with BMI in females (11), no significant associations with BMI were found in our study. To our knowledge, we are the first to report a sexual difference in the distribution of ISM1. In extensively adjusted models, we observed a significant difference in ISM1, with higher levels in men. Sexual differences have also been reported for several adipokines, such as adiponectin, leptin, and vaspin levels (4244). Our finding of sexual dimorphism of ISM1 was further supported by a recent report that ISM1 was only significantly higher in pubertal boys with obesity (45). Such a sex difference in the distribution of ISM1 and other adipokines may partially account for the difference in insulin sensitivity and β-cell function observed between men and women.

Mechanistically, IPH was reported primarily due to deterioration of insulin secretion and, to a lesser extent, to insulin resistance (46, 47). Consistently, no obvious association between circulating ISM1 and insulin sensitivity indices, either based on the basal model or derived from glucose stimulation, was found in our study or the previous study (11, 12, 37, 45). In contrast, a significant positive association of ISM1 levels and β-cell function was found in men by several methodologies, including the Stumvoll indices and oral minimal model test-derived indices, which quantitatively simulate the complex process of glucose metabolism with a mathematical model and analog β-cell responsibility (20, 21). These data were well in line with the previous animal study revealed that ISM1 served as a secreted factor that exogenously and endogenously controls glucose uptake in vitro, which could activate the PI3K-AKT signaling pathway independently of insulin, and insulin-like growth factor receptors improved glucose tolerance (11).Thus, the protective role of ISM1 and IPH might be related to insulin secretion ability. However, more in-depth studies considering the role of ISM1 in insulin secretion are needed in the future.

Our study has several strengths, including the employment of more accurate estimates of pancreatic β-cell function derived from both the 2 h-OGTT test and oral minimal model, such as Stumvoll first and second phase secretion index, PhiD, PhiS and PhiT, other than HOMA-β which only takes into account the basal glucose and insulin concentrations. Such indices based on mathematical models upon glucose stimulation are more sensitive for probing mild defects in β-cell function and provide a more accurate evaluation of the complex process of glucose metabolism (20, 21). Moreover, we also included various sensitivity analyses to ensure the robustness of the findings. Nevertheless, several limitations of this study should be acknowledged. First, our participants were only Chinese recruited from a single center, and thus, further validation in larger and more diverse ethnicities is warranted. Second, despite an extensive adjustment, residual confounding was still possible due to the nature of observational studies. However, the results remained basically the same in various sensitivity analyses. Moreover, considering the large E value that a potentially unmeasured confounding factor would need to deviate from our findings, the E value analysis suggested that no residual confounding factors remained in our study.

In conclusion, our study provides evidence that there is sexual dimorphism in the distribution of circulating ISM1. ISM1 independently associates with β-cell function in men and functions as an easy-to-use tool for the identification of male subjects with high risk of IPH in clinical practice.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

This study complied with the Declaration of Helsinki and was approved by the School of Public Health Institutional Review Board of Sun Yat-sen University. All participants provided written informed consent before data collection. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

JF: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Writing – original draft, Writing – review & editing. JHe: Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing. JZ: Data curation, Investigation, Methodology, Writing – review & editing. JY: Data curation, Investigation, Writing – review & editing. JJ: Data curation, Investigation, Writing – review & editing. JHu: Data curation, Investigation, Writing – review & editing. ZH: Data curation, Investigation, Writing – review & editing. ZZ: Data curation, Investigation, Writing – review & editing. WL: Data curation, Investigation, Writing – review & editing. MX: Resources, Writing – review & editing. YL: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was financially supported by the Natural Science Foundation of Guangdong Province (2021A1515010841), Guangdong Basic and Applied Basic Research Fund Project (Regional Joint Fund-Youth Fund Project, 2019A1515110583); Young Scientists Fund of the National Natural Science Foundation of China (82003455); The Fundamental Research Funds for the Central Universities, Sun Yat-sen University (22qntd4201; 22hytd03). The funders had no role in the preparation, review, or approval of the manuscript or the decision to submit the manuscript for publication.

Acknowledgments

The authors gratefully acknowledge the support of the study participants, study staff, and partner organizations participating in the study survey.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2024.1394190/full#supplementary-material

References

1. Sourij H, Saely CH, Schmid F, Zweiker R, Marte T, Wascher TC, et al. Post-challenge hyperglycaemia is strongly associated with future macrovascular events and total mortality in angiographied coronary patients. Eur Heart J. (2010) 31:1583–90. doi: 10.1093/eurheartj/ehq099

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Barzin M, Hosseinpanah F, Malboosbaf R, Hajsheikholeslami F, Azizi F. Isolated post-challenge hyperglycaemia and risk of cardiovascular events: Tehran Lipid and Glucose Study. Diabetes Vasc Dis Res. (2013) 10:324–9. doi: 10.1177/1479164112469639

CrossRef Full Text | Google Scholar

3. Crandall JP, Shamoon H, Cohen HW, Reid M, Gajavelli S, Trandafirescu G, et al. Post-challenge hyperglycemia in older adults is associated with increased cardiovascular risk profile. J Clin Endocrinol Metab. (2009) 94:1595–601. doi: 10.1210/jc.2008-1829

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Shaw JE, Hodge AM, de Courten M, Chitson P, Zimmet PZ. Isolated post-challenge hyperglycaemia confirmed as a risk factor for mortality. Diabetologia. (1999) 42:1050–4. doi: 10.1007/s001250051269

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Qiao Q, Dekker JM, de Vegt F, Nijpels G, Nissinen A, Stehouwer CD, et al. Two prospective studies found that elevated 2-hr glucose predicted male mortality independent of fasting glucose and HbA1c. J Clin Epidemiol. (2004) 57:590–6. doi: 10.1016/j.jclinepi.2003.10.007

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Silbernagel G, Sourij H, Grammer TB, Kleber ME, Hartaigh B, Winkelmann BR, et al. Isolated post-challenge hyperglycaemia predicts increased cardiovascular mortality. Atherosclerosis. (2012) 225:194–9. doi: 10.1016/j.atherosclerosis.2012.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Yang W, Lu J, Weng J, Jia W, Ji L, Xiao J, et al. Prevalence of diabetes among men and women in China. N Engl J Med. (2010) 362:1090–101. doi: 10.1056/NEJMoa0908292

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Qiao Q, Nakagami T, Tuomilehto J, Borch-Johnsen K, Balkau B, Iwamoto Y, et al. Comparison of the fasting and the 2-h glucose criteria for diabetes in different Asian cohorts. Diabetologia. (2000) 43:1470–5. doi: 10.1007/s001250051557

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Zheng Y, Ley SH, Hu FB. Global aetiology and epidemiology of type 2 diabetes mellitus and its complications. Nat Rev Endocrinol. (2018) 14:88–98. doi: 10.1038/nrendo.2017.151

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Pera EM, Kim JI, Martinez SL, Brechner M, Li SY, Wessely O, et al. Isthmin is a novel secreted protein expressed as part of the Fgf-8 synexpression group in the Xenopus midbrain-hindbrain organizer. Mech Dev. (2002) 116:169–72. doi: 10.1016/S0925-4773(02)00123-5

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Jiang Z, Zhao M, Voilquin L, Jung Y, Aikio MA, Sahai T, et al. Isthmin-1 is an adipokine that promotes glucose uptake and improves glucose tolerance and hepatic steatosis. Cell Metab. (2021) 33:1836–52.e11. doi: 10.1016/j.cmet.2021.07.010

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Wang J, Du J, Ge X, Peng W, Guo X, Li W, et al. Circulating Ism1 Reduces the Risk of Type 2 Diabetes but not Diabetes-Associated NAFLD. Front Endocrinol (Lausanne). (2022) 13:890332. doi: 10.3389/fendo.2022.890332

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Kizer JR, Arnold AM, Benkeser D, Ix JH, Djousse L, Zieman SJ, et al. Total and high-molecular-weight adiponectin and risk of incident diabetes in older people. Diabetes Care. (2012) 35:415–23. doi: 10.2337/dc11-1519

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Luo S, Zhao Y, Zhu S, Liu L, Cheng K, Ye B, et al. Flavonifractor plautii protects against elevated arterial stiffness. Circ Res. (2023) 132:167–81. doi: 10.1161/CIRCRESAHA.122.321975

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Wang Z, Chen Z, Zhang L, Wang X, Hao G, Zhang Z, et al. Status of hypertension in China: Results from the China hypertension survey, 2012-2015. Circulation. (2018) 137:2344–56. doi: 10.1161/CIRCULATIONAHA.117.032380

PubMed Abstract | CrossRef Full Text | Google Scholar

16. NCD Risk Factor Collaboration (NCD-RisC). Worldwide trends in hypertension prevalence and progress in treatment and control from 1990 to 2019: a pooled analysis of 1201 population-representative studies with 104 million participants. Lancet. (2021) 398:957–80. doi: 10.1016/S0140-6736(21)01330-1

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Fan J, Yin S, Lin D, Liu Y, Chen N, Bai X, et al. Association of serum retinol-binding protein 4 levels and the risk of incident type 2 diabetes in subjects with prediabetes. Diabetes Care. (2019) 42:1574–81. doi: 10.2337/dc19-0265

PubMed Abstract | CrossRef Full Text | Google Scholar

18. National Cholesterol Education Program (NCEP) Expert Panel on Detection, Evaluation, and Treatment of High Blood Cholesterol in Adults (Adult Treatment Panel III). Third Report of the National Cholesterol Education Program (NCEP) Expert Panel on Detection, Evaluation, and Treatment of High Blood Cholesterol in Adults (Adult Treatment Panel III) final report. Circulation. (2002) 106:3143–421. doi: 10.1161/circ.106.25.3143

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Levey AS, Stevens LA, Schmid CH, Zhang YL, Castro AF 3rd, Feldman HI, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. (2009) 150:604–12. doi: 10.7326/0003-4819-150-9-200905050-00006

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Cobelli C, Dalla Man C, Toffolo G, Basu R, Vella A, Rizza R. The oral minimal model method. Diabetes. (2014) 63:1203–13. doi: 10.2337/db13-1198

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Chen ME, Chandramouli AG, Considine RV, Hannon TS, Mather KJ. Comparison of β-cell function between overweight/obese adults and adolescents across the spectrum of glycemia. Diabetes Care. (2018) 41:318–25. doi: 10.2337/dc17-1373

PubMed Abstract | CrossRef Full Text | Google Scholar

22. ElSayed NA, Aleppo G, Aroda VR, Bannuru RR, Brown FM, Bruemmer D, et al. 2. Classification and diagnosis of diabetes: Standards of care in diabetes-2023. Diabetes Care. (2023) 46:S19–s40. doi: 10.2337/dc23-S002

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Garaulet M, Lopez-Minguez J, Dashti HS, Vetter C, Hernández-Martínez AM, Pérez-Ayala M, et al. Interplay of dinner timing and MTNR1B type 2 diabetes risk variant on glucose tolerance and insulin secretion: A randomized crossover trial. Diabetes Care. (2022) 45:512–9. doi: 10.2337/dc21-1314

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Wu H, Tremaroli V, Schmidt C, Lundqvist A, Olsson LM, Krämer M, et al. The gut microbiota in prediabetes and diabetes: A population-based cross-sectional study. Cell Metab. (2020) 32:379–90.e3. doi: 10.1016/j.cmet.2020.06.011

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Masafumi Matsuda M, DeFronzo RA. Insulin sensitivity indices obtained from oral glucose tolerance testing,comparison with the euglycemic insulin clamp. Diabetes Care. (1999) 22(9):1462–70. doi: 10.2337/diacare.22.9.1462

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Stumvoll M, Van Haeften T, Fritsche A, Gerich J. Oral glucose tolerance test indexes for insulin sensitivity and secretion based on various availabilities of sampling times. Diabetes Care. (2001) 24:796–7. doi: 10.2337/diacare.24.4.796

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Ohn JH, Kwak SH, Cho YM, Lim S, Jang HC, Park KS, et al. 10-year trajectory of β-cell function and insulin sensitivity in the development of type 2 diabetes: a community-based prospective cohort study. Lancet Diabetes Endocrinol. (2016) 4:27–34. doi: 10.1016/S2213-8587(15)00336-8

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Huang R, Yin S, Ye Y, Chen N, Luo S, Xia M, et al. Circulating retinol-binding protein 4 is inversely associated with pancreatic β-cell function across the spectrum of glycemia. Diabetes Care. (2020) 43:1258–65. doi: 10.2337/dc19-2432

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Zhou X, Qiao Q, Ji L, Ning F, Yang W, Weng J, et al. Nonlaboratory-based risk assessment algorithm for undiagnosed type 2 diabetes developed on a nation-wide diabetes survey. Diabetes Care. (2013) 36:3944–52. doi: 10.2337/dc13-0593

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Pencina MJ, D'Agostino RB Sr, Steyerberg EW. Extensions of net reclassification improvement calculations to measure usefulness of new biomarkers. Stat Med. (2011) 30:11–21. doi: 10.1002/sim.4085

PubMed Abstract | CrossRef Full Text | Google Scholar

31. VanderWeele TJ, Ding P. Sensitivity analysis in observational research: Introducing the E-value. Ann Intern Med. (2017) 167:268–74. doi: 10.7326/M16-2607

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Mathur MB, Ding P, Riddell CA, VanderWeele TJ. Web site and R package for computing E-values. Epidemiology. (2018) 29:e45–e7. doi: 10.1097/EDE.0000000000000864

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Spracklen CN, Horikoshi M, Kim YJ, Lin K, Bragg F, Moon S, et al. Identification of type 2 diabetes loci in 433,540 East Asian individuals. Nature. (2020) 582:240–5. doi: 10.1038/s41586-020-2263-3

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Li W, Xie B, Qiu S, Huang X, Chen J, Wang X, et al. Non-lab and semi-lab algorithms for screening undiagnosed diabetes: A cross-sectional study. EBioMedicine. (2018) 35:307–16. doi: 10.1016/j.ebiom.2018.08.009

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Pencina MJ, D'Agostino RB Sr., D'Agostino RB Jr., Vasan RS. Evaluating the added predictive ability of a new marker: from area under the ROC curve to reclassification and beyond. Stat Med. (2008) 27:157–72; discussion 207-12. doi: 10.1002/sim.2929

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Wang C, Xu M, Feng R, Zhang L, Yin X, Feng R, et al. Serum isthmin-1 levels are positively and independently correlated with albuminuria in patients with type 2 diabetes mellitus. BMJ Open Diabetes Res Care. (2022) 10(5):e002972. doi: 10.1136/bmjdrc-2022-002972

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Liao J, Li Y, Gui X, Zhang Y, Hu X, Cheng L, et al. Serum isthmin-1 was increased in type 2 diabetic patients but not in diabetic sensorimotor peripheral neuropathy. Diabetes Metab Syndr Obes. (2023) 16:2013–24. doi: 10.2147/DMSO.S411127

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Hanefeld M, Koehler C, Henkel E, Fuecker K, Schaper F, Temelkova-Kurktschiev T. Post-challenge hyperglycaemia relates more strongly than fasting hyperglycaemia with carotid intima-media thickness: the RIAD Study. Risk Factors in Impaired Glucose Tolerance for Atherosclerosis and Diabetes. . Diabetes Med. (2000) 17:835–40. doi: 10.1046/j.1464-5491.2000.00408.x

CrossRef Full Text | Google Scholar

39. DECODE Study Group. Age- and sex-specific prevalences of diabetes and impaired glucose regulation in 13 European cohorts. Diabetes Care. (2003) 26:61–9. doi: 10.2337/diacare.26.1.61

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Adam JM, Josten D. Isolated post-challenge hyperglycemia: concept and clinical significance. Acta Med Indones. (2008) 40:171–5.

PubMed Abstract | Google Scholar

41. Barrett-Connor E, Ferrara A. Isolated postchallenge hyperglycemia and the risk of fatal cardiovascular disease in older women and men. The Rancho Bernardo Study. Diabetes Care. (1998) 21:1236–9. doi: 10.2337/diacare.21.8.1236

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Delaney KZ, Santosa S. Sex differences in regional adipose tissue depots pose different threats for the development of Type 2 diabetes in males and females. Obes Rev. (2022) 23:e13393. doi: 10.1111/obr.13393

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Youn BS, Klöting N, Kratzsch J, Lee N, Park JW, Song ES, et al. Serum vaspin concentrations in human obesity and type 2 diabetes. Diabetes. (2008) 57:372–7. doi: 10.2337/db07-1045

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Lew J, Sanghavi M, Ayers CR, McGuire DK, Omland T, Atzler D, et al. Sex-based differences in cardiometabolic biomarkers. Circulation. (2017) 135:544–55. doi: 10.1161/CIRCULATIONAHA.116.023005

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Ruiz-Ojeda FJ, Anguita-Ruiz A, Rico MC, Leis R, Bueno G, Moreno LA, et al. Serum levels of the novel adipokine isthmin-1 are associated with obesity in pubertal boys. World J Pediatr. (2023) 19(9):864–72. doi: 10.1007/s12519-022-00665-8

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Fukushima M, Suzuki H, Seino Y. Insulin secretion capacity in the development from normal glucose tolerance to type 2 diabetes. Diabetes Res Clin Pract. (2004) 66 Suppl 1:S37–43. doi: 10.1016/j.diabres.2003.11.024

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Abdul-Ghani MA, Tripathy D, DeFronzo RA. Contributions of beta-cell dysfunction and insulin resistance to the pathogenesis of impaired glucose tolerance and impaired fasting glucose. Diabetes Care. (2006) 29:1130–9. doi: 10.2337/diacare.2951130

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Isthmin-1, adipokine, isolated post-challenge hyperglycemia, diabetes, sexual dimorphism

Citation: Fan J, He J, Zhu J, Yang J, Ju J, Huang J, Huang Z, Zhang Z, Li W, Xia M and Liu Y (2024) Sex-specific association of circulating Isthmin-1 with isolated post-challenge hyperglycemia. Front. Endocrinol. 15:1394190. doi: 10.3389/fendo.2024.1394190

Received: 01 March 2024; Accepted: 11 July 2024;
Published: 25 July 2024.

Edited by:

Elliot Rayfield, Icahn School of Medicine at Mount Sinai, United States

Reviewed by:

Emily Jane Gallagher, Icahn School of Medicine at Mount Sinai, United States
Tingting Geng, Fudan University, China

Copyright © 2024 Fan, He, Zhu, Yang, Ju, Huang, Huang, Zhang, Li, Xia and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yan Liu, liuyan215@mail.sysu.edu.cn

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.