- 1Department of Geriatric Medicine, Qilu Hospital, Cheeloo College of Medicine, Shandong University, Jinan, China
- 2Key Laboratory of Cardiovascular Proteomics of Shandong Province, Qilu Hospital, Cheeloo College of Medicine, Shandong University, Jinan, China
- 3Jinan Clinical Research Center for Geriatric Medicine (202132001), Jinan, China
- 4Jinan Aixinzhuoer Medical Laboratory, Jinan, China
Introduction: Diabetic sarcopenia (DS) is characterized by muscle atrophy, slower nerve conduction, reduced maximum tension generated by skeletal muscle contraction, and slower contraction rate. Hence, DS can cause limb movement degeneration, slow movement, reduced balance, reduced metabolic rate, falls, fractures, etc. Moreover, the relevant early biological metabolites and their pathophysiological mechanism have yet to be characterized.
Method: The current cross-sectional study employed serum metabolomics analysis to screen potential noninvasive biomarkers in patients with diabetic sarcopenia. A total of 280 diabetic patients were enrolled in the study (n = 39 sarcopenia [DS], n = 241 without sarcopenia [DM]). Ten patients were randomly selected from both groups. Non-targeted metabolomic analysis was performed by ultra-high-performance liquid chromatography-electrospray ionization tandem mass spectrometry.
Results: A total of 632 differential metabolites were identified, including 82 that were significantly differentially abundant (P < 0.05, VIP > 1, FC > 1.2 or FC < 0.8). Compared with the DM group, the contents of pentadecanoic acid, 5'-methylthioadenosine (5'-MTA), N,N-dimethylarginine (asymmetric dimethylarginine, ADMA), and glutamine in the DS group were significantly increased, while that of isoxanthohumol was decreased.
Discussion: Based on receiver operating characteristic curve analysis, pentadecanoic acid, 5'-MTA, ADMA, and glutamine may serve as potential biomarkers of DS. Moreover, ATP-binding cassette (ABC) transporters and the mammalian target of the rapamycin signaling pathway were found to potentially have important regulatory roles in the occurrence and development of DS (P < 0.05). Collectively, the differential metabolites identified in this study provide new insights into the underlying pathophysiology of DS and serve as a basis for therapeutic interventions.
1 Introduction
Type 2 diabetes mellitus (T2DM) is a chronic multifactor metabolic disease. The associated chronic hyperglycemia and insulin resistance can induce downstream metabolic changes, resulting in the imbalance of protein, lipid and carbohydrate metabolism, resulting in cardiovascular and cerebrovascular diseases (1), neurodegenerative diseases (2), muscle atrophy (3) and other systemic multi-organ chronic damage and even failure. Consequent detrimental effects occur in muscle mass and muscle strength, thus, significantly impacting the musculoskeletal system and leading to the development of sarcopenia (4, 5). Indeed, the prevalence of sarcopenia is high in patients with T2DM, ranging from 7% to 29.3% (6). However, the pathogenesis and etiology of diabetic sarcopenia (DS) are complex (7), with nutritional deficiency (8), muscle fat accumulation (9, 10), inflammation (11), oxidative stress (9, 12), lack of exercise (13), and diabetic complications such as diabetic microangiopathy (14), diabetic neuropathy (15), and diabetic foot disease (16) reported as potential etiologic factors. DS is characterized by muscle atrophy, slower nerve conduction, reduced maximum tension generated by skeletal muscle contraction, and gradually slower contraction rates. Moreover, with an increase in age, middle-aged and elderly individuals with DS further experience limb movement degeneration, reduced balance, and metabolic rate, as well as increased risk of falls, fractures, etc. Furthermore, although DS is a complication of T2DM with high incidence and different treatment methods (11), diagnostic criteria have only been established for sarcopenia, while that for DS remains unclear.
Metabolomics analysis effectively detects and quantifies small molecules, including metabolites, within samples to identify differentially abundant metabolites across various disease states or physiological periods (17). In particular, it can be applied to detect early markers of disease, discover drug targets, elucidate disease mechanisms, and facilitate improved disease diagnosis. Currently, the diagnostic criteria for DS primarily focus on assessing skeletal muscle strength, mass, etc., with no established serological indicators to guide diagnosis. In this study, metabolomic analysis is employed to screen serum metabolomic markers of DS and elucidate potential underlying mechanisms associated with DS development and pathogenesis. Hence, the primary aim of this study is to identify targets for early clinical diagnosis and treatment of DS, as well as potential markers capable of predicting disease prognosis and development.
2 Materials and methods
2.1 Study population
A total of 280 patients aged > 50 years were consecutively recruited from the Department of Geriatric Medicine, Qilu Hospital, Shandong University, Shandong Province, China, comprising diabetic patients with sarcopenia (DS, n =39) and without sarcopenia (DM, n = 241) were enrolled. From this patient population, we randomly selected ten patients from the DS and ten from the DM groups. T2DM was defined based on the standard criteria set by the American Diabetes Association of the Endocrinology Department (18). Sarcopenia was defined, based on the criteria recommended by the Asian Working Group for Sarcopenia (AWGS) (19), as low muscle mass with no effect on muscle strength or physical performance. Appendicular skeletal muscle mass (ASM) index (ASM/height2) < 7.0 kg/m2 in men and ASMI < 5.7 kg/m2 in women was defined as low muscle mass, and handgrip strength < 26 kg in men and < 18 kg in women was defined as low muscle strength. Low physical performance was defined as a gait speed < 0.8 m/s. Exclusion criteria included a history of stroke, stents, artificial pacemaker, serious bone, joint, or neuromuscular diseases, mental diseases, metal implants, or a history of malignant tumors, liver disease, end-stage kidney disease, thyroid dysfunction, arthritis, carpal tunnel syndrome, serious diabetic complications (kidney failure, diabetic retinopathy, diabetic foot, etc.), bariatric surgery, or lower limb surgery. Trained interviewers collected all questionnaire data.
This study was approved by the Medical Ethics Committee of the Qilu Hospital of Shandong University (Approval No: 2021274). All methods and analyses were carried out according to approved protocols and guidelines. All study participants provided written informed consent.
A general questionnaire was used to collect data related to age, sex, and medical history. Standardized equipment was used to measure patient height and weight. Participants wore light clothing and removed their shoes during measurements. Body mass index (BMI) was calculated by dividing the weight (kg) by height squared (m)2. The patient was placed in the supine position with an appropriate pressure cuff placed above the anterior cubital fossa to measure systolic blood pressure (BP) and diastolic BP in both arms; the lower right and left systolic BP and diastolic BP were measured using a sphygmomanometer in the supine position. Systolic BP and diastolic BP was measured at the levels of the posterior tibial and tibial arteries. The ankle brachial index (ABI) for both lower limbs was calculated.
2.2 Measurement of skeletal muscle mass, muscle strength, and physical performance
The regional body composition of patients with diabetes was determined using INBODY720. Bioimpedance analysis (BIA) was employed to screen whole-body muscle mass. Body composition was analyzed using a multi-frequency bioelectrical impedance analysis device (InBody (720) body composition analyzer, Biospace Co., Ltd., Seoul, South Korea); participants fasted the night prior and sat quietly within the device for 5 min before measurements were made. Segmental resistance was measured at six frequencies (1, 5, 50, 250, and 500 Hz, and 1 MHz); output was recorded. The ASMI was calculated by dividing the ASM by height squared in meters.
Handgrip strength was used to assess muscle strength. Grip strength was measured using The Jamar dynamometer (Sammons Preston, Inc., Bolingbrook, IL, USA) in the standard postures recommended by the American Association of Hand Therapists (ASHT)(20). Grip strength measurements were recorded for both arms of each patient and repeated thrice until the maximum load was reached.
2.3 Metabolomic analysis
2.3.1 Sample collection
Fasting blood samples were collected between 7:00 and 8:00 a.m. from the antecubital vein and stored in a Progel tube in a cryogenic rack (2°C–8°C). The whole blood was immediately centrifuged at 1000 × g for 10 min at room temperature to separate the serum. Serum samples were immediately stored at -80°C for metabolomic analysis.
2.3.2 Metabolomic analyses
After thawing at 4°C, 100 μL of each sample was extracted. The mixture was prepared using 800 µL of pre-cooled methanol/acetonitrile (1:1, v/v) and 100 µL pre-cooled water, subjected to 1 h ultrasonic treatment in an ice bath, and then incubated with the sample at -20 °C for 2 h. After centrifugation at 4°C and 16,000 × g for 20 min, the supernatant was collected and drained using a high-speed vacuum centrifuge. When mass spectrometry was required, the sample was redissolved with 100 µL of methanol-aqueous solution (1:1, v/v), and the supernatant was collected for sample analysis following centrifugation at 20,000 × g and 4°C for 20 min.
Non-targeted metabolomics analysis was performed using the DS. SHIMADZU-LC30 UHPLC system with an ACQUITY UPLC® HSS T3 column (2.1 × 100 mm, 1.8 µm; Waters, Milford, MA, USA). The samples were separated by UPLC and analyzed by mass spectrometry using a QE Plus mass spectrometer (Thermo Scientific). The HESI source was employed for ionization. Metabolite structure was identified by accurate mass number matching (mass tolerance < 20 ppm) and secondary spectrogram matching (mass tolerance < 0.02 Da). The HMDB, MassBank, and other public databases, as well as the local metabolite standard database, were searched, and the MSDIAL software (RIkagaku KENkyusho/Institute of Physical and Chemical Research, Japan) was employed to conduct retention time correction, peak alignment, and peak area extraction. For the extracted data, the missing values in the group were deleted; 50% ion peaks were not included in subsequent statistical analysis. The total peak area of the positive and negative ion data was normalized, and the peaks were integrated; R software was used for pattern recognition. After the data were subjected to unit variance scaling (UV), subsequent data analysis was conducted.
2.3.3 KEGG enrichment analysis
The Kyoto Encyclopedia of Genes and Genomes (KEGG; www.genome.jp/kegg/pathway.html) database was used to analyze the metabolite data, and Fisher’s exact test was used for enrichment analysis. After performing false discovery rate (FDR) correction for multiple comparisons, the significantly enriched pathways (P < 0.05) were screened.
2.3.4 Differential metabolite identification
All multivariate data analysis and modeling used R (version: 4.0.3) and R packages. The data was centered on the mean, and Pareto scaling was adopted. Principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (OPLS-DA) were used to establish the model. All evaluation models were overfitted via permutation test. At the univariate analysis level, the metabolites were identified based on the statistical significance threshold for the variable influence of the projected (VIP) values obtained in the OPLS-DA model and by a two-tailed Student’s t-test (P-value) for normalized raw data. In multigroup analysis, one-way analysis of variance (ANOVA) was used to calculate P values. Metabolites with VIP values > 1.0, FC > 1.2 or < 0.8, and P values < 0.05 were considered statistically significant. All identified differential metabolites were applied for cluster analysis using R package.
2.4 Laboratory measurements
Fasting blood samples were collected and centrifuged to separate serum (12,000 × g, 10 min, 4°C). Glucose(Glu), triglyceride (TG), total cholesterol (TC), high-density lipoprotein (HDL) cholesterol, low-density lipoprotein (LDL) cholesterol (chemiluminescent immunoassay, Roche), and glycosylated hemoglobin (HbA1c; high-performance liquid chromatography assay, HPLC-873G8, Tosoh) were measured by the laboratory of Qilu Hospital of Shandong University.
2.5 ELISA detection
A total of 280 serum samples were collected according to section 2.3.1 (39 from the DS group and 241 from the DM group) for ELISA using a Human 5′-Methylthioadenosine (MTA) ELISA Kit (MM-61971H, FEIYA Biotechnology, Jiangsu, China) and Human Asymmetric Dimethylarginine (ADMA) ELISA Kit (MM-51745H1, FEIYA Biotechnology). The assay was carried out according to the manufacturer’s instructions. The absorbance (OD value) was measured at a wavelength of 450 nm with a microplate reader, and the concentration of the target substance was calculated using a standard curve.
2.6 Statistical analysis
In this study, SPSS 23.0 software was used for statistical analysis, and the S-W method was used to assess data normality. Measurement data conforming to normal distribution were expressed in the form of “mean ± standard deviation”, and an independent sample t-test was used to analyze differences between the two groups. Count data were expressed as a “composition ratio,” and the Chi-square test was used to analyze the difference between the two groups. All P < 0.05 were considered statistically significant.
3 Results
3.1 Characteristics of study participants
The clinical characteristics of diabetic patients with DS are presented in Table 1. No significant differences were observed in sex, however, significant differences were detected in muscle mass (P < 0.05), grip strength (P < 0.01), and calf/waist/hip circumference (P < 0.05) between the DS and DM groups. Moreover, compared with the DM control group, DS patients were older and had a longer course of T2DM (P < 0.05). DS patients also had lower BMI and DBP, as well as higher pulse pressure (P < 0.05), serum triglyceride concentrations, plasma glucose, and glycosylated hemoglobin levels (P < 0.05 for each). Moreover, serum HDL levels were lower in DS patients compared with those in the DM controls (P < 0.05). Meanwhile, no significant differences were observed in uric acid concentration.
Table 1 Clinical characteristics of participants in the diabetic patients with sarcopenia (DS) and diabetic patients without sarcopenia (DM) control groups.
3.2 Serum metabolic profiling using ultra high-performance liquid chromatography coupled with quadrupole-exactive mass spectrometry
The overall distribution trend of all samples was observed through PCA analysis (Figures 1A, B). The total variance of the data represented by the first two principal components in positive ion mode was 17.14% (Figures 1A). In negative ion mode, the first two principal components represented 21.73% of the total variance in the data (Figure 1B). The PCA score of the positive-ion group showed a significant separation tendency, whereas an insignificant trend was observed in the separation of PCA scores between negative ion groups. Considering the blank supervision of the PCA model, OPLS-DA was subsequently performed to further filter the differential metabolites between groups (Figures 1C, D). The OPLS-DA model (R2Y, Q2) in positive ion mode (R2Y = 0.968, Q2 = 0.729) and negative ion mode (R2Y = 0.992, Q2 = 0.527) revealed that the R2Y and Q2 were ≥ 0.5, indicating that the model was stable and reliable with good predictive capacity. The permutation test of the OPLS-DA model established in positive and negative ion modes is shown in Figures 1E, F. OPLS-DA model overfitting did not occur and exhibited significant differences between sample distributions, indicating the presence of significant differences in serum metabolites between the two groups. From all PCA and OPLS-DA scores, it can be speculated that the occurrence of DS correlates with changes in the levels of certain metabolites.
Figure 1 Multivariate statistical analysis (A) Principal component analysis (PCA) scatter plot of DS (green circle) vs. DM (red circle) in positive ion mode. (B) PCA scatter plot of DS (green circle) vs. DM (red circle) in the negative ion mode. (C) Orthogonal partial least squares discriminant analysis (OPLS-DA) in positive ion mode of DS vs. DM. (D) OPLS–DA in negative ion mode of DS vs. DM. (E) OPLS–DA permutation test in positive ion mode. (F) OPLS–DA permutation test in negative ion mode.
A total of 632 differential metabolites were detected (Figures 2A, B), including 82 significant differential metabolites (P < 0.05, VIP > 1), primarily classified as alkaloids and their derivatives, benzenoids, lignans, neolignans, and related compounds; lipids and lipid-like molecules; nucleosides, nucleotides, and analogs; organic acids and their derivatives; organic nitrogen compounds; organic oxygen compounds; organoheterocyclic compounds; and phenylpropanoids and polyketides. The predominant differential metabolites included D-glucosamine 6-phosphate, D-glucuronic acid, melibiose, N,N-dimethylarginine (asymmetric dimethylarginine, ADMA), etc. Qualitatively significant differentially abundant metabolites were then applied to construct a hierarchical cluster map for each group of samples (Figure 3). The clustering results were then used to assess the stability of the expression of selected target metabolites within each group. Meanwhile, the metabolites within similar clusters were deemed to have identical expression patterns, which might indicate a corresponding reaction step in the metabolic process.
Figure 2 Volcano plot of differentially abundant metabolites in diabetic patients with sarcopenia (DS) DS vs. diabetic patients without sarcopenia (DM). (A) DS vs. DM in positive ion mode. (B) DS vs. DM in negative ion mode. The abscissa is the multiple change value of the expression difference in metabolites between the two groups, and the ordinate is the statistical test value of the expression difference in metabolites, i.e., P-value. The abscissa and ordinate values were processed logarithmically. Each point represents a specific metabolite.
Figure 3 Heat map analysis of 82 differential metabolites between the DS and DM groups. The values of differential metabolites were normalized and shown on a color scale. The high and low metabolite levels are represented as red and blue scales, respectively.
3.3 ROC analysis of potential DS diagnostic biomarkers
The area under the ROC curve (AUC) for d-glucuronic acid and pentadecanoic acid was 0.94 and 0.700, respectively, whereas 5’-methylthioadenosine (5’-MTA) was 0.83. The concentration of ADMA was 0.86, glutamine was 0.75, and isoxanthohumol (IX) was 0.76 (P < 0.05; Table 2 and Figures 4A-E).
Figure 4 Receiver operating characteristic curve (ROC) of differential metabolites. (A–E) ROC curves of five metabolites: pentadecanoic acid, 5’-methylthioadenosine (MTA), N,N-dimethylarginine (asymmetric dimethylarginine, ADMA), glutamine, and isoxanthohumol (IX). AUC represents the recognition performance. The 95% confidence interval of the AUC was calculated based on the nonparametric resampling method. The points on the curve represent the optimal thresholds determined by the ROC curve to distinguish the two groups. The ordinate is the enrichment rate, which is the ratio of the amount of metabolites enriched in the pathway to the amount of metabolites annotated in the pathway.
3.4 Detection and identification of differential metabolites
Compared with the DM group, the concentrations of pentadecanoic acid, 5’-MTA, ADMA, glutamine and gluconic acid in the DS group were significantly increased, while that of IX was decreased(P < 0.05 for each, Figures 5A-E). We further determined the concentration of ADMA and 5’-MTA in the serum of 280 patients using ELISA. ELISA results showed that the concentrations of 5’-MTA and ADMA were significantly increased in the DS group (P < 0.01, Figure 6), which were consistent with the results of metabolomics analysis.
Figure 5 Analysis of differential metabolite abundance. (A) Quantitative differences between DS and DM controls for plasma pentadecanoic acid concentrations; (B) 5’-methylthioadenosine concentrations; (C) N,N-dimethylarginine concentrations; and (D) glutamine concentrations. (E) Isoxanthohumol concentration. Welch t-test was performed to compare the means of each metabolite in participants with DS and DM controls. * P < 0.05 vs. control. Compared with the DM group, the concentrations of hexadecanoic acid, methionine, N,N-dimethylarginine, glutamine, and gluconic acid in the DS group increased, and the concentration of isoxanthohumol decreased.
Figure 6 ELISA analysis of the concentration of serum ADMA and 5’-MTA in DS and DM group. ** P < 0.01 vs. DM group.
3.5 KEGG pathway enrichment analysis
We used KEGG enrichment analysis to identify the relevant metabolic pathways involved in DS, and the results showed that the mammalian target of rapamycin (mTOR) signaling pathway, ABC transporters, amyotrophic lateral sclerosis, amoebiasis, pyrimidine metabolism, as well as D-arginine, and D-ornithine metabolism pathways were the primary pathways associated with DS (P < 0.05, Figure 7). KEGG enrichment analyses were carried out with the Fisher’s exact test, and FDR correction for multiple testing was performed.KEGG analysis further identified ABC transporters (P = 0.005789, Corrected = 0.116732)and the mTOR signaling pathway (P = 0.011146, Corrected = 0.116732) as the most significantly different between the two groups.
Figure 7 Functional analysis of differential metabolites. (A) KEGG enrichment statistics. The x-axis indicates the negative log transformation of the P-value, and the y-axis indicates the name of the KEGG metabolic pathway. The size and color of the bubbles represent the number and degree of enrichment of different metabolites, respectively. (B) Sankey diagram: the first column on the left represents the differential metabolites that were upregulated or downregulated, and the height of the box represents the amount of data flow, that is, the number of metabolites annotated to pathways. The more corresponding pathways, the higher the box. The red streamline represents the flow direction of upregulated (up) metabolites, whereas the blue streamline represents the flow direction of downregulated (down) metabolites. The second, third, and fourth columns represent the pathways, and the hierarchy increases accordingly. M: metabolism, E: environmental information processing, H: human diseases, C00015: uridine 5’-diphosphate, C00062: L-(+)-arginine, C00526 deoxyuridine, C05402: melibiose; BH: biological process.
4 Discussion
In this study, we randomly selected ten DS and ten DM patients for serum untargeted metabolomics analysis. The concentrations of alkanoic acid, 5 ‘-MTA, ADMA, and glutamine were significantly increased in the DS group compared with the DM group, while that of IX was significantly decreased. KEGG analysis further identified ABC transporters and the mTOR signaling pathway as significantly different between the two groups. To our knowledge, this is the first study to report biomarkers for DS using a metabolomics approach. To identify serum metabolomic markers of DS and explore the related mechanisms, muscle testing was performed on enrolled patients and clinical characteristics were statistically analyzed. The results showed that, compared with the DM group, DS patients were older, had lower BMI, and significantly lower grip strength and muscle mass. Moreover, serum laboratory tests showed significant differences in HbA1c, blood glucose, and blood lipid levels between the two groups, which is consistent with previously reported results (21).
Additionally, higher levels of pentadecanoic acid were detected in the DS group compared with the DM group. As a long-chain saturated fatty acid (LCFA), pentadecanoic acid is a potent biological marker of dairy fat intake (22). As LCFAs accumulate, they become deposited in muscle, consequently compromising muscle integrity. In fact, within muscle tissue, continuous accumulation of LCFAs can lead to disordered protein synthesis and lipotoxicity, activation of immune cells, and induction of muscle cell inflammation, resulting in muscle atrophy and apoptosis of muscle proteins (23). Meanwhile, in non-adipose tissue, ectopic accumulation of LCFAs leads to metabolic dysfunction; increased lipid content causes increased lipid load, muscle atrophy, and a corresponding decrease in muscle content and mass (24). It has been suggested that damage to muscle integrity may be related to intramuscular fat (peroxisome proliferator-activated receptor γ)(25) and intramuscular components associated with the resulting muscle inflammation (cyclo-oxygen-ase-2, inducible nitric oxide synthase) (26). This further suggests a relationship between ectopic intramuscular adipose and muscle loss. In addition, LCFAs can lead to insulin resistance, which is one of the primary etiologic factors of sarcopenia (27). Hence, the LCFA pentadecanoic acid accumulates in skeletal muscle, where it promotes the production of inflammatory mediators, insulin resistance, skeletal muscle protein consumption, and muscle fiber atrophy, and decreases skeletal muscle mass and content, which reflects the etiology and pathogenesis of sarcopenia. Accordingly, pentadecanoic acid may represent a predictive biomarker for DS.
5’-MTA levels were also elevated in the DS group compared to the DM group. MTA is a naturally occurring hydrophobic thioadenine nucleoside (28), which is an intermediate metabolite or by-product of the methionine (MET) cycle and polyamine synthetase, both of which are intricately associated with key cellular pathways, including gene expression regulation, proliferation, differentiation and apoptosis (29). Although MET is a substrate for homocysteine, elevated levels of MET can cause hyperhomocysteinemia (30), which has been shown to reduce skeletal muscle cell viability and produce an energy imbalance (31–33). Moreover, elevated plasma homocysteine can lead to different systemic or local neurological diseases and promote inflammation (34). Hence, hyperhomocysteinemia can cause skeletal muscle injury and dysfunction, as well as reduce skeletal muscle development and growth (35). Therefore, an increase in MET increases sarcopenia to a certain extent, which is consistent with our experimental results.
In our study, ADMA concentrations were higher in DS patients than in DM controls. ADMA, a naturally occurring chemical found in the plasma, is a metabolic by-product of the protein modification process in the human cytoplasm and is closely related to L-arginine (L-Arg) (36). L-Arg is a conditionally essential amino acid and precursor of nitric oxide (NO) (37). As a key endogenous regulator of vascular tone and endothelial function, NO plays an important role in the regulation of sympathetic vasoconstriction in the microcirculation of exercise muscles (38). ADMA has also been shown to inhibit NO production by L-arginine via inhibition of NO synthase production (39), thus leading to an imbalance in vascular contraction and diastolic function of skeletal muscle, blood hypoperfusion, tissue ischemia, and hypoxia. Consequently, skeletal muscle protein formation is negatively impacted, cell and skeletal muscle cells begin to atrophy, and the occurrence and development of sarcopenia is exacerbated. In addition, ADMA can promote the development of chronic inflammation, which can cause damage to skeletal muscle cells. Meanwhile, inflammatory factors also impact ADMA homeostasis (40). Regarding DS, in particular, the enhanced inflammatory response and increased activity of inflammatory factors lead to a decrease in the Arg/ADMA ratio (41), which increases ADMA synthesis and exacerbates the development of diabetes-related muscle mass, strength, and functional decline.
We also found that glutamine concentrations were higher in DS patients than in DM controls. Glutamine is the most abundant amino acid in the blood. As a non-essential amino acid, it is mainly catalyzed by glutamine synthase to synthesize ammonia (NH3) and glutamate in skeletal muscle and is involved in myriad metabolic pathways, including signal transduction, nitrogen transport, energy metabolism, etc. (42, 43). Under stimulation by sarcopenia and hyperglycemia, the body enters a state of long-term, low-level inflammation, causing the abundance of leukocytes, such as lymphocytes, to increase. Unlike skeletal muscle, white blood cells do not have glutamine synthase and cannot synthesize glutamine. In addition, macrophages, lymphocytes, and other inflammatory cells will consume glutamine as a substrate. Therefore, to meet the demand for glutamine anabolism under inflammatory conditions, skeletal muscle cells increase their glutamine production, which becomes released into the blood, causing an increase in plasma glutamine concentration. Indeed, Ilaiwy et al. observed that in vitro muscle cell atrophy is related to an increase in glutamine concentration in culture media (44). Moreover, considering that the DS group was older than the DM group and aging is defined as a disease of anabolic resistance to nutrients (such as amino acids), these patients may have exhibited impaired protein synthesis of glutamine, resulting in reduced muscle protein synthesis and development of sarcopenia. Based on these collective study results, we postulate that glutamine concentration is elevated in the DS group and may serve as an effective metabolic marker of DS.
In addition, we found that IX levels were lower in DS patients than in the DM group. IX is a prenyl flavonoid that inhibits the production of secondary bile acids. Secondary bile acids are cytotoxic DNA-damaging metabolites converted from primary bile acids by 7-α-dehydroxylase (45) that interfere with proper metabolism of sugars and lipids, causing lipids to accumulate in muscle tissues, induction of oxidative stress and the inflammatory response, mitochondrial dysfunction, insulin resistance, and damage to skeletal muscle cells, thereby increasing the risk of developing DS (46). Moreover, xanthohumol, a precursor of IX, has been shown to have anti-obesity effects in animal studies (47). Therefore, IX might improve glucose and lipid metabolism to a certain extent, thereby improving lipid accumulation in diabetes and reducing the possibility of DS.
We also found that 3-methylxanthine was highly expressed in the DS group. 3-Methylxanthine is nephrotoxic in that it impacts kidney function, causing kidney damage and increasing the risk of end-stage renal disease (48). Decreased renal function, as an independent risk factor for sarcopenia, suggests that 3-methylxanthine appears to be a predictive metabolite for diabetic myopathy (49).
The level of L-(+)-arginine was also significantly higher in the DS group than in the DM control group. L-(+)-arginine is an inducer of mTOR, which induces muscle protein synthesis primarily through the mTOR pathway (50). However, studies have shown that when sarcopenia occurs, the L-(+)-arginine concentration in muscles begins to decrease (51). Hence, we postulate that our findings may have resulted from protein decomposition in plasma and a subsequent compensatory increase in L-(+)-arginine concentration, reflecting metabolic changes secondary to sarcopenia.
We also observed a significant increase in D-gluconic acid levels within the DS group. D-Gluconic acid is involved in the antioxidative stress response, thereby upregulating glycolysis and downregulating the TCA cycle, leading to mitochondrial dysfunction and muscle ATP deficiency (52, 53). However, due to the lack of further evidence related to sarcopenia, more studies are needed to verify its applicability as a metabolic marker for DS.
We further identified metabolic pathways associated with DS using KEGG enrichment analysis. Among the enriched pathways were ABC transporter superfamily, mTOR signaling pathway, Chagas disease, amoebiasis, amyotrophic lateral sclerosis, pyrimidine metabolism, D-arginine and D-ornithine metabolism, as well as seven other metabolic pathways. Of these pathways, ABC transporters and the mTOR signaling pathway are the most influential and may represent potential target pathways for DS. ABC transporters, one of the largest families of membrane proteins in most organisms, are involved in myriad key physiological processes and serve as pathogenic factors in many diseases. In fact, mutations in ABC-related genes cause neurological disorders and defects in cholesterol and bile transport (54). Cholesterol transport defects affect lipid metabolism and lead to fat deposition in skeletal muscle, thus affecting the activity of skeletal muscle cells and causing skeletal muscle atrophy. Meanwhile, neurological diseases directly affect the activity of innervated muscles, thus impacting communication between muscles and nerves, as well as the types of muscle fibers and expression of myosin (55). This causes muscle weakness and atrophy, which increases the risk of developing sarcopenia. The mTOR signaling pathway plays a key regulatory role in mitochondrial autophagy in skeletal muscles. Mitochondrial quality in skeletal muscle cells is closely associated with sarcopenia; reduced quality deteriorates the muscle microenvironment, leading to the progression of sarcopenia in naturally or rapidly aging mice (56). In addition, L-(+)-arginine also has an important role in the mTOR signaling pathway. Hence, we purport that ABC transporters and the mTOR signaling pathway may contribute to DS development and, therefore, may represent target pathways for the prevention and treatment of DS.
Certain limitations were noted in this study. First, the relatively small sample size limited the statistical ability to definitively identify metabolites associated with DS. Second, due to the cross-sectional design, we were unable to predict DS prognosis or determine the causal relationship between serum differential metabolites and DS. Further studies are warranted to elucidate the clinical correlation between these metabolites and DS, as well as the related mechanisms. Nevertheless, the advantage of our study lies in the application of metabolomics analysis to identify potential metabolites related to the development and progression of DS. These findings may inform the development of effective preventative and therapeutic modalities for DS.
5 Conclusion
In conclusion, this study explored potential biomarkers of DS using a metabolomic strategy. Pentadecanoic acid, 5’-MTA, ADMA, and glutamine were identified as potential biomarkers of DS, while ABC transporters and the mTOR signaling pathway likely make significant contributions toward DS development. The potential relationships between differential metabolites and DS provide novel insights into the mechanisms underlying DS and potential strategies for future treatments.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by the Medical Ethics Committee of the Qilu Hospital of Shandong University. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author contributions
YT Conceptualization, Methodology, Resource, Investigation, Formal analysis, Data curation, Visualization, Writing - Original draft and Writing - Review & Editing. XL Data curation and Formal analysis. YY Data curation and Formal analysis. BL Conceptualization, Validation and Writing - Review & Editing. YF Data curation and Formal analysis. WZ Data curation. CF Data curation. XY Data curation. ZH Formal analysis and Funding acquisition. MC Supervision, Funding acquisition, Resources and Project administration. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Ministry of Science and Technology of the People’s Republic of China (2015FY111600) Science and Technology Development Plan of Shandong Province (2017G006041).
Acknowledgments
We would like to thank Dr. Haoxin Shi and Dr. Wei Ren.
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.
Abbreviations
ABC transporters, ATP-binding cassette transporter; ADA, American Diabetes Association; ADMA, N,N-dimethylarginine/asymmetric dimethylarginine; ANOVA, analysis of variance; ASM, appendicular skeletal muscle mass ; ASMI, appendicular skeletal muscle mass index; AUC, area under the ROC curve; AWGS, Asian Working Group For Sarcopenia; BIA, bioimpedance analysis; BMI, body mass index; BP, blood pressure; CRP, C-reactive protein; DM, diabetes mellitus; DS, diabetic sarcopenia; FC , fold change; 5’-MTA, 5’-methylthioadenosine; HbA1c, glycosylated hemoglobin ; Hcy, homocysteine; HFS, high-fat high-sugar; IGF-1, insulin-like growth factor 1; IL-6, interleukin-6; IX, isoxanthohumol; KEGG, Kyoto Encyclopedia of Genes and Genomes; LCFA, long-chain saturated fatty acid; MET, methionine; mTOR signaling pathway, mammalian target of rapamycin signaling pathway; NO, nitric oxide; OPLS-DA, orthogonal partial least squares discriminant analysis; PCA, principal component analysis; ROC, receiver operating characteristic; ROS, reactive oxygen species; T2DM, type 2 diabetes mellitus; TNF-α, tumor necrosis factor-α; UHPLC-ESI-MS/MS, ultra-high-performance liquid chromatography-electrospray ionization tandem mass spectrometry; UV, unit variance; VIP, variable importance for the projection.
References
1. Muscoli S, Barillà F, Tajmir R, Meloni M, Della Morte D, Bellia A, et al. The new role of SGLT2 inhibitors in the management of heart failure: Current evidence and future perspective. Pharmaceutics (2022) 14(8):1730. doi: 10.3390/pharmaceutics14081730
2. Rea S, Della-Morte D, Pacifici F, Capuani B, Pastore D, Coppola A, et al. Insulin and exendin-4 reduced mutated huntingtin accumulation in neuronal cells. Front Pharmacol (2020) 11:779. doi: 10.3389/fphar.2020.00779
3. Perry BD, Caldow MK, Brennan-Speranza TC, Sbaraglia M, Jerums G, Garnham A, et al. Muscle atrophy in patients with type 2 diabetes mellitus: roles of inflammatory pathways, physical activity and exercise. Exerc Immunol Rev (2016) 22:94–109.
4. Pereira S, Marliss EB, Morais JA, Chevalier S, Gougeon R. Insulin resistance of protein metabolism in type 2 diabetes. Diabetes (2008) 57:56–63. doi: 10.2337/db07-0887
5. Tabara Y, Ikezoe T, Yamanaka M, Setoh K, Segawa H, Kawaguchi T, et al. Advanced glycation end product accumulation is associated with low skeletal muscle mass, weak muscle strength, and reduced bone density: The nagahama study. J Gerontol A Biol Sci Med Sci (2019) 74:1446–53. doi: 10.1093/gerona/gly233
6. Izzo A, Massimino E, Riccardi G, Della Pepa G. A narrative review on sarcopenia in type 2 diabetes mellitus: Prevalence and associated factors. Nutrients (2021) 13(1):183. doi: 10.3390/nu13010183
7. Cruz-Jentoft AJ, Sayer AA. Sarcopenia. Lancet (2019) 393:2636–46. doi: 10.1016/s0140-6736(19)31138-9
8. Freijer K, Nuijten MJ, Schols JM. The budget impact of oral nutritional supplements for disease related malnutrition in elderly in the community setting. Front Pharmacol (2012) 3:78. doi: 10.3389/fphar.2012.00078
9. Kalinkovich A, Livshits G. Sarcopenic obesity or obese sarcopenia: A cross talk between age-associated adipose tissue and skeletal muscle inflammation as a main mechanism of the pathogenesis. Ageing Res Rev (2017) 35:200–21. doi: 10.1016/j.arr.2016.09.008
10. Rhee EJ. Nonalcoholic fatty liver disease and diabetes: An epidemiological perspective. Endocrinol Metab (Seoul) (2019) 34:226–33. doi: 10.3803/EnM.2019.34.3.226
11. Wang M, Tan Y, Shi Y, Wang X, Liao Z, Wei P. Diabetes and sarcopenic obesity: Pathogenesis, diagnosis, and treatments. Front Endocrinol (Lausanne) (2020) 11:568. doi: 10.3389/fendo.2020.00568
12. Andreadi A, Bellia A, Di Daniele N, Meloni M, Lauro R, Della-Morte D, et al. The molecular link between oxidative stress, insulin resistance, and type 2 diabetes: A target for new therapies against cardiovascular diseases. Curr Opin Pharmacol (2022) 62:85–96. doi: 10.1016/j.coph.2021.11.010
13. Gomes MJ, Martinez PF, Pagan LU, Damatto RL, Cezar MDM, Lima ARR, et al. Skeletal muscle aging: influence of oxidative stress and physical exercise. Oncotarget (2017) 8:20428–40. doi: 10.18632/oncotarget.14670
14. Bouchi R, Fukuda T, Takeuchi T, Minami I, Yoshimoto T, Ogawa Y. Sarcopenia is associated with incident albuminuria in patients with type 2 diabetes: A retrospective observational study. J Diabetes Investig (2017) 8:783–7. doi: 10.1111/jdi.12636
15. Yang Q, Zhang Y, Zeng Q, Yang C, Shi J, Zhang C, et al. Correlation between diabetic peripheral neuropathy and sarcopenia in patients with type 2 diabetes mellitus and diabetic foot disease: A cross-sectional study. Diabetes Metab Syndr Obes (2020) 13:377–86. doi: 10.2147/dmso.S237362
16. Cheng Q, Hu J, Yang P, Cao X, Deng X, Yang Q, et al. Sarcopenia is independently associated with diabetic foot disease. Sci Rep (2017) 7:8372. doi: 10.1038/s41598-017-08972-1
17. Vuckovic D. Current trends and challenges in sample preparation for global metabolomics using liquid chromatography-mass spectrometry. Anal Bioanal Chem (2012) 403:1523–48. doi: 10.1007/s00216-012-6039-y
18. American Diabetes Association of the Endocrinology Department. Diagnosis and classification of diabetes mellitus. Diabetes Care (2010) 33 Suppl 1:S62–9. doi: 10.2337/dc10-S062
19. Chen LK, Liu LK, Woo J, Assantachai P, Auyeung TW, Bahyah KS, et al. Sarcopenia in Asia: consensus report of the Asian working group for sarcopenia. J Am Med Dir Assoc (2014) 15:95–101. doi: 10.1016/j.jamda.2013.11.025
20. Fess EE. Clinical assessment recommendations. American Society of Hand Therapists Monograph (1981) Retrieved October 20, 2022, from https://asht.org/practice/clinical-assessment-recommendations.
21. He Q, Wang X, Yang C, Zhuang X, Yue Y, Jing H, et al. Metabolic and nutritional characteristics in middle-aged and elderly sarcopenia patients with type 2 diabetes. J Diabetes Res (2020) 2020:6973469. doi: 10.1155/2020/6973469
22. Wolk A, Furuheim M, Vessby B. Fatty acid composition of adipose tissue and serum lipids are valid biological markers of dairy fat intake in men. J Nutr (2001) 131:828–33. doi: 10.1093/jn/131.3.828
23. Collins KH, Paul HA, Hart DA, Reimer RA, Smith IC, Rios JL, et al. A high-fat high-sucrose diet rapidly alters muscle integrity, inflammation and gut microbiota in Male rats. Sci Rep (2016) 6:37278. doi: 10.1038/srep37278
24. Tardif N, Salles J, Guillet C, Tordjman J, Reggio S, Landrier JF, et al. Muscle ectopic fat deposition contributes to anabolic resistance in obese sarcopenic old rats through eIF2α activation. Aging Cell (2014) 13:1001–11. doi: 10.1111/acel.12263
25. Janani C, Ranjitha Kumari BD. PPAR gamma gene–a review. Diabetes Metab Syndr (2015) 9:46–50. doi: 10.1016/j.dsx.2014.09.015
26. Fink LN, Costford SR, Lee YS, Jensen TE, Bilan PJ, Oberbach A, et al. Pro-inflammatory macrophages increase in skeletal muscle of high fat-fed mice and correlate with metabolic risk markers in humans. Obes (Silver Spring) (2014) 22:747–57. doi: 10.1002/oby.20615
27. Wein S, Wolffram S, Schrezenmeir J, Gasperiková D, Klimes I, Seböková E. Medium-chain fatty acids ameliorate insulin resistance caused by high-fat diets in rats. Diabetes Metab Res Rev (2009) 25:185–94. doi: 10.1002/dmrr.925
28. Vandenbark AA, Ferro AJ, Barney CL. Inhibition of lymphocyte transformation by a naturally occurring metabolite: 5'-methylthioadenosine. Cell Immunol (1980) 49:26–33. doi: 10.1016/0008-8749(80)90052-0
29. Avila MA, García-Trevijano ER, Lu SC, Corrales FJ, Mato JM. Methylthioadenosine. Int J Biochem Cell Biol (2004) 36:2125–30. doi: 10.1016/j.biocel.2003.11.016
30. Singh M, George AK, Eyob W, Homme RP, Stansic D, Tyagi SC. High-methionine diet in skeletal muscle remodeling: epigenetic mechanism of homocysteine-mediated growth retardation. Can J Physiol Pharmacol (2021) 99:56–63. doi: 10.1139/cjpp-2020-0093
31. Van Schoor NM, Swart KM, Pluijm SM, Visser M, Simsek S, Smulders Y, et al. Cross-sectional and longitudinal association between homocysteine, vitamin B12 and physical performance in older persons. Eur J Clin Nutr (2012) 66:174–81. doi: 10.1038/ejcn.2011.151
32. Kolling J, Scherer EB, Siebert C, Hansen F, Torres FV, Scaini G, et al. Homocysteine induces energy imbalance in rat skeletal muscle: is creatine a protector? Cell Biochem Funct (2013) 31:575–84. doi: 10.1002/cbf.2938
33. Swart KM, Enneman AW, Van Wijngaarden JP, Van Dijk SC, Brouwer-Brolsma EM, Ham AC, et al. Homocysteine and the methylenetetrahydrofolate reductase 677C–>T polymorphism in relation to muscle mass and strength, physical performance and postural sway. Eur J Clin Nutr (2013) 67:743–8. doi: 10.1038/ejcn.2013.97
34. Da Cunha AA, Ferreira AG, Loureiro SO, Da Cunha MJ, Schmitz F, Netto CA, et al. Chronic hyperhomocysteinemia increases inflammatory markers in hippocampus and serum of rats. Neurochem Res (2012) 37:1660–9. doi: 10.1007/s11064-012-0769-2
35. Moore CW, Allen MD, Kimpinski K, Doherty TJ, Rice CL. Reduced skeletal muscle quantity and quality in patients with diabetic polyneuropathy assessed by magnetic resonance imaging. Muscle Nerve (2016) 53:726–32. doi: 10.1002/mus.24779
36. Teerlink T. HPLC analysis of ADMA and other methylated l-arginine analogs in biological fluids. J Chromatogr B Analyt Technol BioMed Life Sci (2007) 851:21–9. doi: 10.1016/j.jchromb.2006.07.024
37. Siasos G, Tousoulis D, Antoniades C, Stefanadi E, Stefanadis C. L-arginine, the substrate for NO synthesis: an alternative treatment for premature atherosclerosis? Int J Cardiol (2007) 116:300–8. doi: 10.1016/j.ijcard.2006.04.062
38. Chavoshan B, Sander M, Sybert TE, Hansen J, Victor RG, Thomas GD. Nitric oxide-dependent modulation of sympathetic neural control of oxygenation in exercising human skeletal muscle. J Physiol (2002) 540:377–86. doi: 10.1113/jphysiol.2001.013153
39. Hörster I, Weigt-Usinger K, Carmann C, Chobanyan-Jürgens K, Köhler C, Schara U, et al. The l-arginine/NO pathway and homoarginine are altered in duchenne muscular dystrophy and improved by glucocorticoids. Amino Acids (2015) 47:1853–63. doi: 10.1007/s00726-015-2018-x
40. Zoccali C, Maas R, Cutrupi S, Pizzini P, Finocchiaro P, Cambareri F, et al. Asymmetric dimethyl-arginine (ADMA) response to inflammation in acute infections. Nephrol Dial Transplant (2007) 22:801–6. doi: 10.1093/ndt/gfl719
41. Van Der Zwan LP, Scheffer PG, Dekker JM, Stehouwer CD, Heine RJ, Teerlink T. Systemic inflammation is linked to low arginine and high ADMA plasma levels resulting in an unfavourable NOS substrate-to-inhibitor ratio: the hoorn study. Clin Sci (Lond) (2011) 121:71–8. doi: 10.1042/cs20100595
42. Gleeson M. Dosing and efficacy of glutamine supplementation in human exercise and sport training. J Nutr (2008) 138:2045s–9s. doi: 10.1093/jn/138.10.2045S
43. Meynial-Denis D. Glutamine metabolism in advanced age. Nutr Rev (2016) 74:225–36. doi: 10.1093/nutrit/nuv052
44. Ilaiwy A, Quintana MT, Bain JR, Muehlbauer MJ, Brown DI, Stansfield WE, et al. Cessation of biomechanical stretch model of C2C12 cells models myocyte atrophy and anaplerotic changes in metabolism using non-targeted metabolomics analysis. Int J Biochem Cell Biol (2016) 79:80–92. doi: 10.1016/j.biocel.2016.08.012
45. Ridlon JM, Kang DJ, Hylemon PB. Bile salt biotransformations by human intestinal bacteria. J Lipid Res (2006) 47:241–59. doi: 10.1194/jlr.R500013-JLR200
46. Kuno T, Hirayama-Kurogi M, Ito S, Ohtsuki S. Reduction in hepatic secondary bile acids caused by short-term antibiotic-induced dysbiosis decreases mouse serum glucose and triglyceride levels. Sci Rep (2018) 8(1253). doi: 10.1038/s41598-018-19545-1
47. Fukizawa S, Yamashita M, Fujisaka S, Tobe K, Nonaka Y, Murayama N. Isoxanthohumol, a hop-derived flavonoid, alters the metabolomics profile of mouse feces. Biosci Microbiota Food Health (2020) 39:100–8. doi: 10.12938/bmfh.2019-045
48. Sellman R, Klemi PJ. Kidney toxicity of 3-methylxanthine in the rat. J Appl Toxicol (1984) 4:304–7. doi: 10.1002/jat.2550040605
49. Yu MD, Zhang HZ, Zhang Y, Yang SP, Lin M, Zhang YM, et al. Relationship between chronic kidney disease and sarcopenia. Sci Rep (2021) 11:20523. doi: 10.1038/s41598-021-99592-3
50. Bauchart-Thevret C, Cui L, Wu G, Burrin DG. Arginine-induced stimulation of protein synthesis and survival in IPEC-J2 cells is mediated by mTOR but not nitric oxide. Am J Physiol Endocrinol Metab (2010) 299:E899–909. doi: 10.1152/ajpendo.00068.2010
51. Da Boit M, Tommasi S, Elliot D, Zinellu A, Sotgia S, Sibson R, et al. Sex differences in the associations between l-arginine pathway metabolites, skeletal muscle mass and function, and their responses to resistance exercise, in old age. J Nutr Health Aging (2018) 22:534–40. doi: 10.1007/s12603-017-0964-6
52. Krüger A, Grüning NM, Wamelink MM, Kerick M, Kirpy A, Parkhomchuk D, et al. The pentose phosphate pathway is a metabolic redox sensor and regulates transcription during the antioxidant response. Antioxid Redox Signal (2011) 15:311–24. doi: 10.1089/ars.2010.3797
53. Sato E, Mori T, Mishima E, Suzuki A, Sugawara S, Kurasawa N, et al. Metabolic alterations by indoxyl sulfate in skeletal muscle induce uremic sarcopenia in chronic kidney disease. Sci Rep (2016) 6:36618. doi: 10.1038/srep36618
54. Baker A, Carrier DJ, Schaedler T, Waterham HR, Van Roermund CW, Theodoulou FL. Peroxisomal ABC transporters: functions and mechanism. Biochem Soc Trans (2015) 43:959–65. doi: 10.1042/bst20150127
55. Mège RM, Goudou D, Giaume C, Nicolet M, Rieger F. Is intercellular communication via gap junctions required for myoblast fusion? Cell Adhes Commun (1994) 2:329–43. doi: 10.3109/15419069409014208
Keywords: diabetes, sarcopenia, metabolomics, biomarkers, serum metabolites
Citation: Tan Y, Liu X, Yang Y, Li B, Yu F, Zhao W, Fu C, Yu X, Han Z and Cheng M (2023) Metabolomics analysis reveals serum biomarkers in patients with diabetic sarcopenia. Front. Endocrinol. 14:1119782. doi: 10.3389/fendo.2023.1119782
Received: 09 December 2022; Accepted: 08 March 2023;
Published: 22 March 2023.
Edited by:
Sandeep Mathur, Sawai ManSingh Medical College, IndiaReviewed by:
Giuseppe Della Pepa, University of Naples Federico II, ItalyMengran Du, Chongqing Medical University, China
Francesca Pacifici, University of Rome Tor Vergata, Italy
Aikaterini Andreadi, University of Rome Tor Vergata, Italy
Copyright © 2023 Tan, Liu, Yang, Li, Yu, Zhao, Fu, Yu, Han and Cheng. 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: Mei Cheng, am5jbTY1QGVtYWlsLnNkdS5lZHUuY24=