- 1Co-construction Collaborative Innovation Center for Chinese Medicine, Respiratory Diseases by Henan and Education Ministry of P. R. China, Henan University of Chinese Medicine, Zhengzhou, China
- 2The First Affiliated Hospital of Henan University of Chinese Medicine, Zhengzhou, China
- 3Henan Key Laboratory of Chinese Medicine for Respiratory Disease, Henan University of Chinese Medicine, Zhengzhou, China
Background: Traditional Chinese medicine (TCM) is effective for the treatment of acute exacerbation of chronic obstructive pulmonary disease (AECOPD); however, there is no objective index for the evaluation of TCM syndrome efficacy. This study aimed to screen biomarkers related to the efficacy of TCM syndrome using metabolomics.
Methods: We recruited AECOPD patients with phlegm-heat congesting lung (PH)/phlegm-damp amassing lung (PD) syndrome and treated them with Chinese herbal medicine (Qingre Huatan or Zaoshi Huatan granules) in addition to conventional medicine for 7 days. Data on clinical symptoms and sign scores, modified British Medical Research Council (mMRC), COPD assessment test (CAT), and inflammation indicators, including white blood cell (WBC) count, percentage of neutrophil count (NEU%), and C-reactive protein (CRP), were collected before and after treatment to evaluate the therapeutic effect. Serum samples were collected before and after treatment for metabolomic analysis to screen differential metabolites.
Results: A total of 69 patients with AECOPD were enrolled, including 41 and 28 patients in the PH and PD groups, respectively. The clinical symptoms and sign scores, CAT, mMRC, NEU%, and CRP levels after treatment were lower than those before treatment in both groups (p < 0.05). Serum metabolomics analysis showed that there were 13 differential metabolites in the PH group and 16 differential metabolites in the PD group before and after treatment (p < 0.05, variable importance projection (VIP) ≥ 1.00). In the PH group, lysophosphatidylcholine (LPC) (16:0), LPC (17:1), LPC (18:3), LPC (18:2), and LPC (17:0) negatively correlated with clinical symptoms and sign scores (p < 0.05); LPC (16:0), LPC (17:1), LPC (16:1), and LPC (17:0) negatively correlated with WBC (p < 0.05) and NEU% (p < 0.05); and LPC (16:0) negatively correlated with CRP levels. In the PD group, L-phenylalanine positively correlated with CRP levels (p < 0.05), and 2-methylbutyroylcarnitine positively correlated with clinical symptoms and sign (p < 0.05) and CAT scores (p < 0.05). DL-carnitine positively correlated with clinical symptoms and sign scores (p < 0.05).
Conclusion: Serum metabolites may be potential indicators to objectively evaluate the efficacy of TCM syndromes; however, further large controlled trials are required to verify these findings.
1 Introduction
Chronic obstructive pulmonary disease (COPD) is a worldwide public health challenge due to its high prevalence, disability and mortality rates (GBD, 2015 Chronic Respiratory Disease Collaborators, 2017; GBD, 2015 Disease and Injury Incidence and Prevalence Collaborators, 2016; GOLD, 2023). In China, the overall prevalence of COPD was 13.7% in adults aged 40 years or older in 2015, while the total number of patients was nearly 100 million (Wang et al., 2018). COPD is the third leading cause of death globally and ranked the seventh as a cause of disease burden as measured by disability-adjusted life years (DALYs) in 2019 (GBD Chronic Respiratory Disease Collaborators, 2020). Acute exacerbation of COPD (AECOPD) is a major event affecting the natural course of COPD. Patients with COPD have about 0.5–3.5 acute exacerbations per year, which have a serious negative impact on the quality of life, disease progression, and socioeconomic burden of the patients. According to evidence-based guidelines, early detection and timely treatment of deterioration are key to ensuring the best outcomes and reducing the burden of chronic obstructive pulmonary disease (Expert group on Acute Exacerbations of Chronic Obstructive Pulmonary Disease, 2017; GOLD, 2023).
Studies have shown that Chinese herbal medicines are safe and effective in the treatment of AECOPD, and that the treatment of AECOPD by integrating Chinese and Western medicine is superior to conventional Western medicine in terms of increasing clinical efficiency, shortening hospital stay, relieving symptoms, improving lung function, and alleviating inflammation (Meng, 2017; Bai et al., 2022). Therefore, Traditional Chinese medicine (TCM) can be used as a supplementary or alternative treatment for patients with AECOPD (Li et al., 2020a; Ma, 2022; GOLD, 2023). TCM syndromes, as the basis for the diagnosis and treatment of diseases, summarize the etiology, location, disease, and condition of the body at a certain stage in disease development, which is the key index for evaluating the clinical efficacy of TCM. Changes in the multilevel metabolic networks of the body are important (Tao et al., 2018; Tao et al., 2020). The evaluation index of the curative effect of TCM syndromes is the main component of the TCM clinical curative effect evaluation system (Guo et al., 2017). At present, the main criterion for judging the efficacy of TCM syndromes is to evaluate the degree of improvement in TCM syndromes before and after treatment. The TCM syndrome scores and syndrome efficacy evaluation scale are widely used to evaluate the intervention effect of TCM; however, due to many problems, including unclear positioning of the measured outcome level, confusion between diagnosis and evaluation attributes of TCM syndromes, and difficulty in measuring dynamic changes of TCM syndromes, the objectivity and scientificity of TCM syndromes efficacy evaluation has always been controversial (Hu, 2012).
Metabolomics technology can accurately describe the dynamic changes of in endogenous small-molecule metabolites in the body to understand the overall functional state of the body. It has phase, dynamic, and integral characteristics, which are consistent with the TCM syndrome concept (Liu et al., 2022). It has been widely used for disease diagnosis (Wang et al., 2013), differential diagnosis (Deja et al., 2014), and therapeutic interventions (Ren et al., 2016). Several studies have explored biomarkers and pathophysiological mechanisms associated with the diagnosis and classification of TCM syndromes in AECOPD (Yang et al., 2018; Zhang et al., 2020). However, there is a lack of biomarker research on the efficacy of TCM syndromes.
We hypothesized that AECOPD patients with different syndromes have different metabolic profiles before and after treatment, and that metabolomic analysis based on LC-MS can screen biomarkers related to the evaluation of TCM syndrome efficacy. The most common TCM syndromes in patients with AECOPD are phlegm-heat congesting lung (PH) and phlegm-damp amassing lung (PD) (Li et al., 2017). Our previous researches showed that Qingre Huatan granules and Zaoshi Huatan granules have good clinical efficacy in the treatment of PH and PD in AECOPD and can significantly improve clinical symptoms, improve lung function, and reduce inflammatory reactions (Li et al., 2010; Li et al., 2010; Li et al., 2010; Li et al., 2011; Li et al., 2020b). Therefore, this study analyzed the serum metabolomics of AECOPD patients with PH and PD syndrome to identify biomarkers with high specificity and sensitivity and to provide an objective basis for TCM syndrome differentiation and treatment and TCM syndrome efficacy evaluation of AECOPD.
2 Materials and methods
2.1 Polulation
This study was approved by the Ethics Committee of the First Affiliated Hospital of Henan University of Traditional Chinese Medicine according to the Declaration of Helsinki (2017HL-086-01). All the participants signed an informed consent form before participanting the study.
The subjects were AECOPD patients aged from 40 to 80 years old in the general ward of the respiratory department of the First Affiliated Hospital of Henan University of Chinese Medicine from August 2018 to September 2019. All subjects were diagnosed according to “Global Initiative for Chronic Obstructive Lung Disease (GOLD): global strategy for the diagnosis, management, and prevention of COPD (2017 REPORT)” (GOLD, 2017) and “Expert Consensus on the Diagnosis and Treatment of Acute Exacerbation of Chronic Obstructive Pulmonary Disease (AECOPD) in China (Updated in 2017)” (Expert group on Acute Exacerbations of Chronic Obstructive Pulmonary Disease, 2017), and divided into the PH group and PD group based on TCM syndrome differentiation referred to the “Diagnostic Criteria for TCM Syndromes of Chronic Obstructive Pulmonary Disease (2011 edition)” (Li et al., 2012).
The exclusion criteria were as follows: a) pregnant and lactating women; b) patients with severe cardiac insufficiency; c) patients with bronchial asthma, bronchiectasis, active tuberculosis, or diffuse panbronchiolitis; d) patients with pneumothorax, pleural effusion, or pulmonary embolism; e) patients with definite neoplastic disease; f) patients with severe liver and kidney disease; g) participating in clinical trials of other interventions.
2.2 Intervention
All the participants received conventional Western medical treatment, including controlled oxygen therapy, bronchodilators, glucocorticoids, antibiotics and so on. The specific treatment measures were individualized according to “Expert Consensus on the Diagnosis and Treatment of Acute Exacerbation of Chronic Obstructive Pulmonary Disease (AECOPD) in China (Updated in 2017)” (Expert group on Acute Exacerbations of Chronic Obstructive Pulmonary Disease, 2017). The treatment principles are shown in Supplementary Material. They were administered the Chinese herbal Qingre Huatan granules in the PH group and Zaoshi Huatan granules in the PD group. The Chinese herbal medicine granules were supplied by Jiangyin Tian Jiang Pharmaceutical Co., Ltd. Jiangsu, China. Chinese herbal medicine granules were manufactured in strict accordance with the standards of the Chinese Pharmacopoeia (2010). Each granule type was dissolved in warm water and administered orally at a dose of 10 g twice daily for 14 days. The components of the Chinese herbal medicine granules list in Table 1.
2.3 Clinical efficacy evaluation index
Clinical symptoms and sign scores, modified British Medical Research Council (mMRC), COPD assessment test (CAT), and inflammatory indicators, including white blood cell (WBC) count, percentage of neutrophil count (NEU%), and C-reactive protein (CRP) levels, were recorded before and after treatment to evaluate the therapeutic effect. The occurrence of adverse events and electrocardiogram (ECG) and laboratory examinations, such as routine tests of blood, urine, and stool and function tests of the liver and kidney, were used as safety indicators observed and recorded before and after treatment of AECOPD.
2.4 Statistical analysis
SPSS software (version 23.0, IBM, United States) was used to build the database and perform statistical analyses. Categorical variables were presented as numbers (percentages). The paired t-test was used to assess the continuous variables of a normal distribution, which were described by mean (M) ± standard deviation (SD). The continuous variables of non-normal distribution were described as median (interquartile range [IQR]) and analyzed using the Wilcoxon signed-rank test. A p-value of less than 0.05 (p < 0.05) was deemed statistically significant for results.
2.5 Metabolomics analysis
2.5.1 Biochemical and instrumental investigations
Formic acid (FA), acetonitrile (ACN), and ammonium acetate (NH4OAc) were obtained from Thermo Fisher, United States. Tryptophan-d5, valine-d8, carnitine C8:0-d3, carnitine C16:0-d3, cholic acid-d4, and stearic acid-d3 (98% purity) were purchased from Cambridge Isotope Laboratories (Tewksbury, MA, United States). Lysophosphatidylcholine LPC (19:0) was obtained from Avanti Polar Lipids (Alabaster, AL, United States). Pure water was purchased from Wahaha ((Hangzhou, China).
2.5.2 Serum samples preparation
Fasting blood of all subjects was collected in the early morning before and after treatment and allowed to clot naturally at room temperature, which was then centrifuged at 3,000 r/min at 4°C for 15 min to obtain the serum. The samples were kept at −80°C until analysis.
Frozen serum samples were thawed at 4°C. 100 μL serum was transferred, and 400 μL ACN containing seven internal standards was added to remove the proteins. After centrifugation at 13,000 r/min and 4°C for 15 min, two aliquots (200 μL for each) of the supernatant were transferred and then dried in a vacuum centrifuge. The residues were reconstituted with 40 μL 20% ACN before LC-MS analysis.
Quality control (QC) samples were prepared by mixing 10 μL of each serum sample and pretreated as the real samples. At the beginning of the sequence, eight QC samples were run to balance the instrument system and one QC sample was run every six real samples to monitor the stability of the analysis (Luo et al., 2016).
2.5.3 LC-MS analysis
A Dionex Ultimate 3,000 ultra-performance liquid chromatograph (UPLC) coupled to a Q-Exactive Mass Spectrometer (MS) (Thermo Fisher Scientific, Bremen, Germany) via an electrospray ionization source was used for metabolomic analysis. An ACQUITY BEH C18 column (2.1 mm × 100 mm, 1.7 µm) (Waters, Milford, MA, United States) was employed for chromatographic separation. Mobile phases A and B were 0.1% FA in H2O and 0.1% FA in ACN in the electrospray ionization-positive (ESI+) mode, and 5 mM ammonium acetate in water and 5 mM ammonium acetate in 95% ACN in the electrospray ionization-negative (ESI-) mode.
Gradient elution for both modes of analysis was performed as follows: maintained at 5% B for 1 min, then linearly increased to 100% B for 24 min, followed by a 4 min wash at 100% B and a 5 min equilibration at 5% B. The column temperature and flow rate were set at 30°C and 0.35 mL/min, respectively. The injection volume was 2 μL.
MS acquisition was conducted in ESI+ and ESI- modes with the spray voltage at 3.5 and 2.8 kV respectively. The capillary temperature, flow rate of the sheath gas and auxiliary gas, and S-lens RF level were set respectively at 350°C, 40 arbitrary units, 10 arbitrary units, and 55 for both modes. Sequence analysis was performed in full MS mode from 70 to 1,000 m/z at the resolution of 70,000. For identification analysis, the full MS-ddMS2 mode was used, and the parameter settings were as follows: resolution of MS and MS2, 70,000 and 17,500; Top N, 5; dynamic exclusion time, 6 s; normalized collision energy, 20, 40, 60; selection of the top 10 ions for secondary scanning; and dynamic elimination time, 10 s.
2.5.4 Data processing and differential metabolites screening
Raw data were processed using Compound Discoverer 2.1 (Thermo Fisher Scientific Co., Ltd., United States) with an untargeted workflow that included retention time alignment, compound detection, grouping, background annotation, formula prediction, and database search. The peak list was exported and re-integrated using TraceFinder (Thermo Fisher Scientific Co., Ltd., United States). The missing values were removed according to the “80% rule,” the rest of which were filled by the sequential k-nearest neighbor algorithm. The peak area of each compound in the real samples was calibrated using an internal standard, which was selected based on the performance of the internal standards [with minimum relative standard deviation (RSD)] in the QC samples.
The Wilcoxon signed-rank test was first conducted to discern differential compounds using SPSS (version 23.0, IBM, United States) with p < 0.05 defined as statistically significant. Multivariate statistical analyses, including principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA) were performed with SIMCA 14.0 (Umetrics AB, Umea, Sweden). Differential metabolites were further screened based on the variable importance projection (VIP) ≥ 1.00. Permutation tests of 200 times were performed to assess the stability of the model and risk of overfitting.
2.6 Identification of differential metabolites and analysis of metabolic pathways
For the identification of differential metabolites, candidates from the m/zCloud database in the search results of Compound Discoverer 2.1 were checked manually. Both the deviation of the precursor ions and matching of the MS/MS spectra were considered. For differential metabolites with similarity matching with carnitine or LPC, the MS/MS spectra were checked for the characterized production of 85.0291 or 184.0733 for carnitine or LPC, respectively, and the precursor ions were searched against HMDB (https://www.hmdb.ca) and LipidMaps (https://www.lipidmaps.org) with a mass tolerance of 10 ppm to determine their structures. Finally, some of the compounds were further validated by comparing their retention times, precursors, and fragment ions with those of the standards. Metabolic pathway analyses were performed via the MetaboAnalyst 5.0 (https://www.metaboanalyst.ca).
3 Results
3.1 Characteristics of study subjects
We recruited 112 patients with AECOPD, who were divided into the PH group (n = 60) and the PD group (n = 52) according to their TCM syndromes. However, 43 patients (19 in the PH group and 24 in the PD group) did not have blood samples collected after treatment for personal reasons; therefore, 69 patients were finally included in the study. Baseline characteristics of all patients are shown in Table 2.
3.2 Evaluation of clinical efficacy
In the PH group, the clinical symptoms and sign, mMRC, and CAT scores of patients after treatment were significantly lower than those before treatment (p < 0.05). The NEU% and CRP levels after treatment were significantly lower than those before treatment (p < 0.05). (Table 3).
In the PD group, the clinical symptoms and sign, mMRC, and CAT scores of the patients after treatment were lower than before treatment (p < 0.05). The NEU% and CRP levels after treatment were significantly lower than those before treatment (p < 0.05). (Table 4).
3.3 Screening of differential metabolites
PCA was performed on the QC samples to evaluate data quality. The results showed that QC samples clustered well in both ESI+ and ESI- modes (as shown in Figure 1), indicating that the instrument system was stable and the analysis results had good reproducibility.
FIGURE 1. The principal component analysis (PCA) score plots compare the QC samples with all the study samples. (A) PCA score plot in ESI + mode; (B) PCA score plot in ESI- mode.
A total of 3,095 peaks in ESI + mode and 1,168 peaks in ESI- mode were extracted using Compound Discover 2.1. After removing compounds with a RSD larger than 30%, a total of 2,438 compounds in ESI + mode and 854 compounds in ESI- mode were obtained for statistical analysis. After the Wilcoxon signed-rank test, we screened 401 and 127 differential compounds for ESI+ and ESI- modes before and after treatment in the PH group, 366 and 162 differential compounds for ESI+ and ESI- modes, respectively, before and after treatment in the PD group (p < 0.05).
The results of PCA showed a certain degree of separation in serum samples before and after treatment in the PH (Figures 2A, B) and PD groups (Figures 3A, B), although the difference was not significant. PLS-DA was used to screen for differential metabolic compounds before and after treatment. As shown in Figures 2C, E and Figures 3C, E, the serum samples before and after treatment in each group were clearly separated. The PLS-DA models were further validated using permutation tests, and the results showed that the Q2 intercepts of the permutation tests were all less than zero, indicating that the PLS-DA models were valid and reliable without overfitting (Figures 2D, F; Figures 3D, F). According to the standard of VIP ≥1.00, 131 differential compounds (ESI+) and 37 differential compounds (ESI-) were screened out in the PH group, while 95 differential compounds (ESI+) and 36 differential compounds (ESI-) in the PD group before and after treatment respectively.
FIGURE 2. Principal component analysis (PCA) score plots, partial least squares discriminant analysis (PLS-DA) score plots and permutation test of serum metabolic profiling in pre-treatment (blue squares) and post-treatment (red triangles) of the PH group. (A) PCA score plot in ESI + mode; (B) PCA score plot in ESI- mode; (C) PLS-DA score plot in ESI + mode. The model parameters were R2X = 0.221, R2Y = 0.604, Q2 = 0.323; (D) Permutation test of PLS-DA score plot ESI + mode. The parameters were R2 = 0.408, Q2 = −0.200; (E) PLS-DA score plot in ESI- mode. The model parameters were R2X = 0.353, R2Y = 0.355, Q2 = 0.107; (F) Permutation test of PLS-DA score plot ESI- mode. The parameters were R2 = 0.321, Q2 = −0.127.
FIGURE 3. Principal component analysis (PCA) score plots, partial least squares discriminant analysis (PLS-DA) score plots and permutation test of serum metabolic profiling in pre-treatment (blue squares) and post-treatment (red triangles) of the PD group. (A) PCA score plot in ESI + mode; (B) PCA score plot in ESI- mode; (C) PLS-DA score plot in ESI + mode. The model parameters were R2X = 0.393, R2Y = 0.815, Q2 = 0.653; (D) Permutation test of PLS-DA score plot ESI + mode. The parameters were R2 = 0.606, Q2 = −0.175; (E) PLS-DA score plot in ESI- mode. The model parameters were R2X = 0.456, R2Y = 0.634, Q2 = 0.517; (F) Permutation test plots of PLS-DA score plot ESI- mode. The parameters were R2 = 0.404, Q2 = −0.019.
By searching and comparing through various databases, a total of 13 different metabolites in the PH group (Table 5) and 16 different metabolites in the PD group (Table 6) before and after treatment were identified. The results of metabolites enrichment showed that in the PH group, the levels of lysophosphatidylcholine (LPC) (14:0), LPC (15:0), LPC (16:0), LPC (17:1), lysophosphatidylethanolamine (LPE) (18:2), LPC (18:3), LPC (18:2), LPC (16:1), and LPC (17:0) after treatment were upregulated compared to before treatment, while the levels of acetyl-L-carnitine, 3-hydroxyoctanoyl carnitine, hexanoylcarnitine, and 9-decenoylcarnitine after treatment were downregulated than before treatment. In the PD group, among the 16 differential metabolites before and after treatment, the levels of LPC (14:0) and LPC (15:0) were lower after treatment, whereas the levels of the other 14 differential metabolites increased after treatment (Figure 4).
FIGURE 4. Heatmap of differential metabolites. (A) PH group; (B) PD group. The heatmap was created using the balanced log2-transformed and normalized data with 41 cases in the PH group and 28 cases in the PD group before and after treatment. Each column represents one sample. Columns with a red line on top represent samples before treatment, and the columns with a green line on top represent samples after treatment. In the subsequent rows, red hues represent increased concentrations, and blue hues decreased concentrations. Color intensity increases proportionally to the magnitude of the change.
3.4 Analysis of metabolic pathways
Pathway enrichment analysis of differential metabolite showed that the metabolism involved in treatment in the PH group mainly included β-oxidation of very long-chain fatty acids, oxidation of branched-chain fatty acids, phospholipid biosynthesis, and arachidonic acid metabolism (Figure 5A). The metabolic pathways involved in treatment in the PD group were β-oxidation of very long-chain fatty acids, oxidation of branched-chain fatty acids, carnitine synthesis, purine metabolism, mitochondrial β-oxidation of long chain saturated fatty acids, phenylalanine, and tyrosine metabolism, fatty acid metabolism, arachidonic acid metabolism (Figure 5B).
FIGURE 5. Bubble map of metabolic pathway enrichment. Bubble size indicates the number of differential metabolites involved in this pathway. The bubble color represents the p-value of this metabolic pathway. (A) PH group; (B) PD group.
3.5 Correlation of differential metabolites and indicators of efficacy evaluation
In the PH group, LPC (15:0) (r = −0.403, p < 0.01), LPC (16:0) (r = −0.381, p < 0.05), LPC (17:1) (r = −0.470 p < 0.01), LPC (18:3) (r = −0.366, p < 0.05), LPC (18:2) (r = −0.311, p < 0.05), and LPC (17:0) (r = −0.405, p < 0.01) significantly negatively correlated with clinical symptoms and sign scores. Hexanoylcarnitine (r = 0.477, p < 0.01) and 9-decenoylcarnitine (r = 0.427, p < 0.01) levels negatively correlated with mMRC. LPC (16:0) (r = −0.494, p < 0.01), LPC (17:1) (r = −0.470, p < 0.01), LPC (16:1) (r = −0.495, p < 0.01), and LPC (17:0) (r = −0.484, p < 0.01) negatively correlated with WBC. LPC (16:0) (r = −0.457, p < 0.05), LPC (17:1) (r = −0.494, p < 0.01), LPC (16:1) (r = −0.461, p < 0.05), and LPC (17:0) (r = −0.426, p < 0.05) negatively correlated with NEU%. LPC (16:0) (r = −0.372, p < 0.05) negatively correlated with CRP levels. There was no significant correlation between the other differential metabolites and efficacy evaluation indicators in the phlegm heat group (Table 7).
In the PD group, correlation analysis between differential metabolites and indicators of efficacy evaluation showed that LPC (15:0) (r = −0.393, p < 0.05) and LPC (14:0) (r = −0.563, p < 0.01) negatively correlated with clinical symptoms and sign scores. 2-Methylbutyroylcarnitine (r = 0.480, p < 0.01) and DL-carnitine (r = 0.459, p < 0.05) positively correlated with clinical symptoms and sign scores. Acefylline (r = 0.382, p < 0.05) and 2-methylbutyroylcarnitine (r = 0.387, p < 0.05) positively correlated with mMRC. LPC (15:0) negatively correlated with CAT scores (r = 0.388, p < 0.05). LPC (15:0) (r = −0.444, p < 0.05) negatively correlated with CRP levels. L-phenylalanine (r = 0.463, p < 0.05) positively correlated with CRP levels. The remaining differential metabolites in the PD group were not significantly correlated with the efficacy indicators (Table 8).
4 Discussion
With the wide application of metabolomics technology to improve the objectivity and accuracy of TCM syndrome diagnosis and efficacy evaluation, research on biomarkers based on metabolomics has become a popular topic. Zhang et al. found that the expression of 15 substances including lactic acid, succinic acid, lipoic acid, neohesperid, D-erythro-sphingosine, glucose, methyl linoleate, asparagine, and urea was significantly upregulated, while threonine and taurine were significantly downregulated in patients with AECOPD phlegm-turbid syndrome (Zhang et al., 2020). Another study found that there are 50 kinds of small-molecule differential metabolites in the serum of patients with the syndrome of phlegm-heat congesting lung and syndrome of non-phlegm-heat congesting lung of AECOPD, involving metabolic disorders such as sphingolipid, glycerophospholipid, and linoleic acid metabolisms (Sun, 2019). Although several studies have shown that metabolites in patients with different TCM syndromes of AECOPD are different from those in the stable stage or in healthy people, there are still few reports on the biomarkers related to the efficacy evaluation of TCM syndromes. Therefore, on the basis of evaluating the clinical efficacy of syndrome differentiation and treatment of TCM combined with routine treatment of Western medicine, this study conducted a metabolomic analysis of the serum of AECOPD patients with the syndromes of PH and PD through metabolomics and initially explored the biomarkers related to the evaluation of TCM syndrome efficacy of AECOPD, providing an objective basis for the evaluation of syndrome efficacy.
In this study, 23 differential metabolites were screened and identified. Among them, LPC (14:0), LPC (15:0), 3-hydroxyoctanoyl carnitine, hexanoylcarnitine, acetyl-L-carnitine, and 9-decenoylcarnitine were common metabolites in both groups, and the trends of these metabolites was consistent before and after treatment in both groups. This may be due to changes in the metabolic profile caused by conventional Western medicine treatment. Additionally, acefylline is a metabolite of the bronchodilator doxofylline; therefore, it is not considered a biomarker of TCM syndrome efficacy.
In the PH group, the characteristic differential metabolites were LPCs, among them LPC (16:0) and LPC (17:0) are saturated acyl LPCs, and LPC (16:1), LPC (17:1), LPC (18:2), and LPC (18:3) are unsaturated acyl LPCs. Combined with the correlation analysis of clinical efficacy evaluation indicators, we speculated that LPC (16:0), LPC (17:0), LPC (16:1), LPC (17:1), LPC (18:2), and LPC (18:3) may be potential biomarkers of AECOPD with syndrome of PH. These LPCs are involved in glycerophospholipid metabolism and arachidonic acid metabolism, and the serum levels post-treatment of these LPCs were higher than the pre-treatment levels. This result is consistent with a previous study, that showed that the concentrations of serum LPCs, such as LPC (14:0), LPC (16:0), LPC (16:1), LPC (18:0), LPC (18:2), and LPC (18:3) in patients with AECOPD gradually increases after 10–14 days of conventional treatment (Gai et al., 2021). Studies have shown that LPC concentrations vary in different diseases. For example, the concentration of LPC increases in inflammatory diseases (Liu et al., 2020) or adrenoleukodystrophy (Huffnagel et al., 2017), but decreases in cancer (Kim et al., 2014) or infectious diseases (Ahn et al., 2017). LPCs not only exert a pro-inflammatory effect by activating a variety of signaling pathways involved in oxidative stress and inflammation, but also play an anti-inflammatory and anti-infective effects by enhancing natural immunity and adaptive immune responses (Piccirillo et al., 2019; Liu et al., 2020). The pro-inflammatory activity of saturated acyl LPC, such as LPC (16:0), can be antagonized by polyunsaturated acyl LPC, such as LPC (16:1) (Hung et al., 2012; Liu et al., 2020). The pathological change of COPD is chronic airway inflammation and AECOPD is typically associated with increased local or systemic inflammation caused by respiratory infections and air pollution (Expert group on Acute Exacerbations of Chronic Obstructive Pulmonary Disease, 2017; Christenson et al., 2022; GOLD, 2023). Abnormalities in arachidonic acid metabolism, glycerophospholipid metabolism, glycosylphosphatidylinositol-anchored biosynthesis, and linoleic acid metabolism may be related to chronic airway inflammation, vascular endothelial contraction, oxidative stress, and energy depletion in patients with AECOPD. Among the various AECOPD syndromes, the inflammatory reaction in the syndrome of PH is the most intense. Sphingomyelin metabolism and the glycerol phospholipid metabolic pathway may constitute the metabolic pathways of AECOPD PH syndrome disturbances (Sun, 2019). In this study, TCM syndrome differentiation and treatment achieved a therapeutic effect by correcting the disorders of glycerol phospholipid and arachidonic acid metabolism caused by AECOPD, which is consistent with the results of previous studies.
In the PD group, most of the differential metabolites were medium-chain acyl carnitines, such as hexanoylcarnitine, 3-hydroxyoctanoyl carnitine, acetyl-L-carnitine, decanoylcarnitine, 2-octenoylcarnitine, L-octanoylcarnitine, 2-methylbutyroyl carnitine, and 9-decenoylcarnitine, which play key roles in lipid metabolism and beta-oxidation. Correlation analysis showed that DL-carnitine was significantly and positively correlated with clinical symptoms and sign scores and 2-methylbutyroyl carnitine was positively correlated with clinical symptoms and sign scores, and CAT scores. Carnitine, a non-essential amino acid and a quaternary ammonium compound, is concentrated in skeletal and cardiac muscles, as well as in other tissues that metabolize fatty acids as an energy source. The main role of carnitine is to assist in the transportation of fatty acids to the mitochondria for β-oxidation and energy supply. During this process, the intermediate metabolite fatty acyl carnitine, which play an important role in the transport and oxidation of fatty acids, is produced. Acetyl-L-carnitine, hexanoylcarnitine, and 3-hydroxyoctanoyl carnitine are essential for fatty acid metabolism (Capitanio et al., 2020). Carnitine has two isomers (D-carnitine and L-carnitine), both of which have biological activity; however, only L-carnitine is naturally present in animals, whereas D-carnitine is toxic, because it inhibits the activity of the L-forms. Most carnitine in the body comes from dietary sources such as red meat and dairy products. Plasma carnitine concentrations correlate with the nutritional intake of carnitine in healthy adults (Lennon et al., 1986; Steiber et al., 2004). Lower levels of L-carnitine indicate more severe the heart failure, and a higher degree of inflammation, and severe arrhythmia can quickly consume stored carnitine (Vescovo et al., 2005; Dambrova et al., 2021). In addition to β oxidation, carnitine can also act as a free radical scavenger to reduce oxidative stress-induced apoptosis. Reinke, S. N. et al. through a metabolomic analysis of the urine and sputum of patients with asthma, found that carnitines showed the greatest variation associated with asthma severity (Reinke et al., 2022). In patients with COPD, the carnitine level is insufficient, and the degree of deficiency is related to the severity of COPD (Conlon et al., 2016). In the AECOPD, due to hypoxia, the heart rate and oxygen consumption of the myocardium and peripheral respiratory muscles increase. This causes the carnitine stored in the body to be consumed rapidly, ultimately leading to a decrease in carnitine content. Our study showed that after treatment of 7 days, the carnitine content was lower than before treatment. A study showed that changes in carnitine levels in patients with asthma are not affected by glucocorticoid therapy (Reinke et al., 2022). Therefore, the reduction of carnitine levels in patients with AECOPD after treatment in our study may be related to the effect of Chinese medicine components, which may be further studied and verified as the target of Chinese medicine treatment in the future.
In the PD group, we detected another potential biomarker, L-phenylalanine, an alpha-amino acid, which is highly concentrated in many high-protein foods such as meat; L-phenylalanine is one of the 20 amino acids required for protein synthesis, involved in phenylalanine and tyrosine metabolism, and is present in all organisms. In this study, the post-treatment serum phenylalanine level was lower than the pre-treatment level, which is inconsistent with the results of a previous study in which the level of plasma phenylalanine in patients with AECOPD was lower than that in the stable stage (Inoue and Ikeda, 2019). The reason for this inconsistency may be the small sample size or grouping of TCM syndromes. In addition, we analyzed the correlation between L-phenylalanine and curative effects, and inflammatory indicators. The results showed that L-phenylalanine was significantly and positively correlated with CRP levels. Based on the above findings, L-carnitine, 2-methylbutyroyl carnitine, and L-phenylalanine may be potential biomarkers related to the efficacy of AECOPD with syndrome of PD.
In this study, there are still some limitations, such as not elucidating the detailed mechanism of action of these biomarkers. In future studies, network pharmacology and multi-omics combination methods will be used to further validate these potential markers and their mechanisms of action. This can not only provide an objective index for the evaluation of the curative effect of AECOPD syndrome differentiation of TCM, but also provide guidance for improving the composition of TCM prescriptions and improving clinical efficacy.
5 Conclusion
The findings of our study showed that the differential metabolites of the PH and PD groups were significantly different, indicating that the effects of syndrome differentiation treatment on the changes in the metabolism of patients with different TCM syndromes are different. The treatment of the syndrome of phlegm-heat congesting lung by Qingre Huatan granules mainly involves glycerophospholipid metabolism, and LPC (16:0), LPC (16:1), LPC (17:0), LPC (17:1), LPC (18:2), and LPC (18:3) may be potential biomarkers related to the efficacy of syndrome of PH. The treatment of the syndrome of phlegm-damp amassing lung by Zaoshi Huatan granules mainly involves fatty acid oxidation and carnitine biosynthesis, and L-carnitine, 2-methylbutyroyl carnitine, and L-phenylalanine may be potential biomarkers related to the efficacy of syndrome of PD.
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by the Ethics Committee of the First Affiliated Hospital of Henan University of Traditional Chinese Medicine. 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
ZH: Conceptualization, Project administration, Supervision, Writing–review and editing. SY: Data curation, Formal Analysis, Investigation, Writing–original draft. DY: Data curation, Formal Analysis, Methodology, Visualization, Writing–review and editing. LX: Data curation, Formal Analysis, Writing–review and editing. LJ: Conceptualization, Project administration, Supervision, Validation, Writing–review and editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was funded by the National Natural Science Foundation of China (No. 81973781 and 81603556), the construction project of the Characteristic Backbone Discipline of Chinese Medicine in Henan Province, China (No. STG-ZYX02-202108 and STG-ZYXKY-2020007), and Special Scientific Research of Chinese Medicine in Henan Province, China (No. 2022ZY1037).
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/fphar.2024.1302950/full#supplementary-material
References
Ahn, W. G., Jung, J. S., Kwon, H. Y., and Song, D. K. (2017). Alteration of lysophosphatidylcholine-related metabolic parameters in the plasma of mice with experimental sepsis. Inflammation 40, 537–545. doi:10.1007/s10753-016-0500-6
Bai, C., Jiang, D., Wang, L., Wang, H., Li, Y., and Chen, O. (2022). Effectiveness of integrated traditional Chinese and western medicine in treating AECOPD with phlegm—dampness stagnating in the lung syndrome: a Bayesian network meta-analysis. ShangHai J. tradi. Chin. Med. 56, 10–16. doi:10.16305/j.1007-1334.2022.1906058
Capitanio, J. P., Dethloff, F., and Turck, C. W. (2020). Lipid metabolism is associated with temperament, corticosteroid, and hematological measures in infant rhesus monkeys (Macaca mulatta). Zool. Res. 41, 709–714. doi:10.24272/j.issn.2095-8137.2020.244
Christenson, S. A., Smith, B. M., Bafadhel, M., and Putcha, N. (2022). Chronic obstructive pulmonary disease. Lancet 399, 2227–2242. doi:10.1016/S0140-6736(22)00470-6
Conlon, T. M., Bartel, J., Ballweg, K., Günter, S., Prehn, C., Krumsiek, J., et al. (2016). Metabolomics screening identifies reduced L-carnitine to be associated with progressive emphysema. Clin. Sci. (Lond) 130, 273–287. doi:10.1042/CS20150438
Dambrova, M., Zuurbier, C. J., Borutaite, V., Liepinsh, E., and Makrecka-Kuka, M. (2021). Energy substrate metabolism and mitochondrial oxidative stress in cardiac ischemia/reperfusion injury. Free Radic. Biol. Med. 165, 24–37. doi:10.1016/j.freeradbiomed.2021.01.036
Deja, S., Porebska, I., Kowal, A., Zabek, A., Barg, W., Pawelczyk, K., et al. (2014). Metabolomics provide new insights on lung cancer staging and discrimination from chronic obstructive pulmonary disease. J. Pharm. Biomed. Anal. 100, 369–380. doi:10.1016/j.jpba.2014.08.020
Expert group on Acute Exacerbations of Chronic Obstructive Pulmonary Disease (2017). Expert Consensus on the diagnosis and treatment of acute exacerbation of chronic obstructive pulmonary disease (AECOPD) in China. Int. J. Respir. 37, 1041–1057. doi:10.3760/cma.j.issn.1673-436X.2017.14.001
Gai, X., Guo, C., Zhang, L., Zhang, L., Abulikemu, M., Wang, J., et al. (2021). Serum glycerophospholipid profile in acute exacerbation of chronic obstructive pulmonary disease. Front. Physiol. 12, 646010. doi:10.3389/fphys.2021.646010
GBD 2015 Disease and Injury Incidence and Prevalence Collaborators (2016). Global, regional, and national incidence, prevalence, and years lived with disability for 310 diseases and injuries, 1990-2015: a systematic analysis for the Global Burden of Disease Study 2015. Lancet 388, 1545–1602. doi:10.1016/S0140-6736(16)31678-6
GBD Chronic Respiratory Disease Collaborators (2017). Global, regional, and national deaths, prevalence, disability-adjusted life years, and years lived with disability for chronic obstructive pulmonary disease and asthma, 1990-2015: a systematic analysis for the Global Burden of Disease Study 2015. Lancet Respir. Med. 5, 691–706. doi:10.1016/S2213-2600(17)30293-X
GBD Chronic Respiratory Disease Collaborators (2020). Prevalence and attributable health burden of chronic respiratory diseases, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet Respir. Med. 8, 585–596. doi:10.1016/S2213-2600(20)30105-3
GOLD (2017). Global initiative for chronic obstructive lung disease(GOLD):Global strategy for the diagnosis, management and prevention of chronic obstructive pulmonary disease. Available at: https://goldcopd.org/wp-content/uploads/2017/02/wms-GOLD-2017-FINAL.pdf.
GOLD (2023). Global initiative for chronic obstructive lung disease (GOLD): global strategy for the diagnosis, management and prevention of chronic obstructive pulmonary disease. Available at: https://goldcopd.org/2023-gold-report-2/.
Guo, L., Han, J., and Liu, N. (2017). Discussion on the bottleneck of clinical curative effect evaluation system of TCM syndromes under modern medical architecture. China J. Tradit. Chin. Med. Pharm. 32, 3357–3360.
Hu, J. (2012). Discussion on the key scientific problems in the evaluation of syndrome therapeutic effect. World. Sci. technol./moder. Tradit. Chin. Med. Mater. Medica 14, 1242–1245. doi:10.3969/j.issn.1674-3849.2012.01.019
Huffnagel, I. C., Van De Beek, M. C., Showers, A. L., Orsini, J. J., Klouwer, F. C. C., Dijkstra, I. M. E., et al. (2017). Comparison of C26:0-carnitine and C26:0-lysophosphatidylcholine as diagnostic markers in dried blood spots from newborns and patients with adrenoleukodystrophy. Mol. Genet. Metab. 122, 209–215. doi:10.1016/j.ymgme.2017.10.012
Hung, N. D., Sok, D. E., and Kim, M. R. (2012). Prevention of 1-palmitoyl lysophosphatidylcholine-induced inflammation by polyunsaturated acyl lysophosphatidylcholine. Inflamm. Res. 61, 473–483. doi:10.1007/s00011-012-0434-x
Inoue, S., and Ikeda, H. (2019). Differences in plasma amino acid levels in patients with and without bacterial infection during the early stage of acute exacerbation of COPD. Int. J. Chron. Obstruct Pulmon Dis. 14, 575–583. doi:10.2147/COPD.S188422
Kim, S. C., Kim, M. K., Kim, Y. H., Ahn, S. A., Kim, K. H., Kim, K., et al. (2014). Differential levels of L-homocysteic acid and lysophosphatidylcholine (16:0) in sera of patients with ovarian cancer. Oncol. Lett. 8, 566–574. doi:10.3892/ol.2014.2214
Lennon, D. L., Shrago, E. R., Madden, M., Nagle, F. J., and Hanson, P. (1986). Dietary carnitine intake related to skeletal muscle and plasma carnitine concentrations in adult men and women. Am. J. Clin. Nutr. 43, 234–238. doi:10.1093/ajcn/43.2.234
Li, B., Hou, Z., Li, J., Zhao, J., Chen, P., Li, S., et al. (2010a). Huatan Qingre prescription treatment of chronic obstructive pulmonary disease with acute exacerbation and syndrome of heat-phelgm stagnating the lung the effects of cytokines. Liaoning J. Tradit. Chin. Med. 37, 403–405. doi:10.13192/j.ljtcm.2010.03.24.lib.002
Li, J., Li, B., Yu, X., Chen, P., Hou, Z., Li, S., et al. (2011). Treatment of acute embittering chronic obstructive pulmonary disease with heat-clearing and phlegm-dissolving formulas: clinical observation on 69 cases with phlegm heat in the lung. J. Tradit. Chin. Med. 52, 203–207. doi:10.13288/j.11-2166/r.2011.03.029
Li, J., Li, S., Chen, P., Li, B., Ma, L., Yu, H., et al. (2010b). Clinical evaluation of dampness and phlegm prescriptions treatment of chronic obstructive pulmonary disease with acute exacerbation and syndrome of phlegm obstructing lung. Liaoning J. Tradit. Chin. Med. 37, 1721–1725. doi:10.13192/j.ljtcm.2010.09.93.lijsh.028
Li, J., Li, S., and Wang, Z. (2012). Diagnostic criteria for TCM syndromes of chronic obstructive pulmonary disease (2011 edition). J. Chin. Med. 53, 177–178. doi:10.13288/j.11-2166/r.2012.02.025
Li, J., Yu, X., and Xie, Y. (2020b). International clinical practice guideline of Chinese medicine chronic obstructive pulmonary disease. World Chin. Med. 15, 1084–1092. doi:10.3969/j.issn.1673-7202.2020.07.026
Li, J., Zhang, H., Ruan, H., Si, Y., Sun, Z., Liu, H., et al. (2020a). Effects of Chinese herbal medicine on acute exacerbations of COPD: a randomized, placebo-controlled study. Int. J. Chron. Obstruct Pulmon Dis. 15, 2901–2912. doi:10.2147/COPD.S276082
Li, J., Zhang, H., Wang, H., Cao, F., Hou, C., Pan, Y., et al. (2017). Clinical investigation on the characteristics of TCM syndrome evolution of chronic obstructive pulmonary disease. J. Tradit. Chin. Med. 58, 772–776. doi:10.13288/j.11-2166/r.2017.09.014
Li, S., Chen, P., Li, J., Li, B., Hou, Z., Zhao, J., et al. (2010c). Dampness and phlegm prescriptions treatment of chronic obstructive pulmonary disease with acute exacerbation and syndrome of phlegm obstructing lung the effects of cytokines. J. Chin. Med. 25, 122–124. doi:10.16368/j.issn.1674-8999.2010.01.047
Liu, P., Zhu, W., Chen, C., Yan, B., Zhu, L., Chen, X., et al. (2020). The mechanisms of lysophosphatidylcholine in the development of diseases. Life Sci. 247, 117443. doi:10.1016/j.lfs.2020.117443
Liu, X., Chen, M., Li, L., Wei, H., Song, W., and Lao, H. (2022). Current situation and prospect of metabonomics in study of lung system diseases and syndromes. Shandong J. Tradit. Chin. Med. 41, 1142–1147. doi:10.16295/j.cnki.0257-358x.2022.10.021
Luo, P., Yu, H., Zhao, X., Bao, Y., Hong, C. S., Zhang, P., et al. (2016). Metabolomics study of roux-en-Y gastric bypass surgery (RYGB) to treat type 2 diabetes patients based on ultraperformance liquid chromatography-mass spectrometry. J. Proteome Res. 15, 1288–1299. doi:10.1021/acs.jproteome.6b00022
Ma, Y. (2022). Research progress on traditional Chinese medicine in the treatment of acute exacerbation chronic obstructive pulmonary disease. Guangming J. Chin. Medi. 37, 4008–4011. doi:10.3969/j.issn.1003-8914.2022.21.057
Meng, Y. (2017). Systematic review of randomized controlled trials of TCM syndrome differentiation for acute exacerbation of chronic obstructive pulmonary disease. Master Thesis. Beijing: University of Chinese Medicine.
Piccirillo, A. R., Hyzny, E. J., Beppu, L. Y., Menk, A. V., Wallace, C. T., Hawse, W. F., et al. (2019). The lysophosphatidylcholine transporter MFSD2A is essential for CD8(+) memory T cell maintenance and secondary response to infection. J. Immunol. 203, 117–126. doi:10.4049/jimmunol.1801585
Reinke, S. N., Naz, S., Chaleckis, R., Gallart-Ayala, H., Kolmert, J., Kermani, N. Z., et al. (2022). Urinary metabotype of severe asthma evidences decreased carnitine metabolism independent of oral corticosteroid treatment in the U-BIOPRED study. Eur. Respir. J. 59, 2101733. doi:10.1183/13993003.01733-2021
Ren, X., Ma, S., Wang, J., Tian, S., Fu, X., Liu, X., et al. (2016). Comparative effects of dexamethasone and bergenin on chronic bronchitis and their anti-inflammatory mechanisms based on NMR metabolomics. Mol. Biosyst. 12, 1938–1947. doi:10.1039/c6mb00041j
Steiber, A., Kerner, J., and Hoppel, C. L. (2004). Carnitine: a nutritional, biosynthetic, and functional perspective. Mol. Asp. Med. 25, 455–473. doi:10.1016/j.mam.2004.06.006
Sun, Y. (2019). Preliminary study on the characteristics of metabolic pattern of phlegm-heat obstructed lung syndrome of AECOPD based on UPLC-MS technique. Master Thesis. Beijing: University of Chinese Medicine.
Tao, J., Wang, S., Jiang, M., and Yuan, P. (2020). Progress and dilemma of clinical research on metabolomics of TCM syndromes. China J. Tradit. Chin. Med. Pharm. 35, 3799–3803.
Tao, J., Yuan, B., and Wang, S. (2018). Technical methods and application status of traditional Chinese medicine syndrome research. China J. Tradit. Chin. Med. Pharm. 33, 2982–2985.
Vescovo, G., Ravara, B., Gobbo, V., and Dalla Libera, L. (2005). Inflammation and perturbation of the l-carnitine system in heart failure. Eur. J. Heart Fail 7, 997–1002. doi:10.1016/j.ejheart.2004.11.010
Wang, C., Xu, J., Yang, L., Xu, Y., Zhang, X., Bai, C., et al. (2018). Prevalence and risk factors of chronic obstructive pulmonary disease in China (the China Pulmonary Health [CPH] study): a national cross-sectional study. Lancet 391, 1706–1717. doi:10.1016/S0140-6736(18)30841-9
Wang, L., Tang, Y., Liu, S., Mao, S., Ling, Y., Liu, D., et al. (2013). Metabonomic profiling of serum and urine by (1)H NMR-based spectroscopy discriminates patients with chronic obstructive pulmonary disease and healthy individuals. PLoS One 8, e65675. doi:10.1371/journal.pone.0065675
Yang, X., Huang, F., Wang, Q., Wu, M., Zhang, Y., Wang, X., et al. (2018). Serum metabonomic features in acute exacerbation of chronic obstructive pulmonary disease with syndrome of phlegm and heat accumulated in lung with 1-HNMR modality. J. Chin. Med. 33, 719–724. doi:10.16368/j.issn.1674-8999.2018.05.172
Keywords: biomarkers, AECOPD, efficacy evaluation, TCM syndrome, metabolomics
Citation: Hailong Z, Yimei S, Yan D, Xinguang L and Jiansheng L (2024) Exploration of biomarkers for efficacy evaluation of traditional Chinese medicine syndromes of acute exacerbation of chronic obstructive pulmonary disease based on metabolomics. Front. Pharmacol. 15:1302950. doi: 10.3389/fphar.2024.1302950
Received: 18 October 2023; Accepted: 15 January 2024;
Published: 26 January 2024.
Edited by:
Yuqin Chen, First Affiliated Hospital of Guangzhou Medical University, ChinaReviewed by:
Shiyan Dong, University of Texas MD Anderson Cancer Center, United StatesWei Gao, Tongji University School of Medicine, China
Copyright © 2024 Hailong, Yimei, Yan, Xinguang and Jiansheng. 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: Li Jiansheng, bGlfanM4QDE2My5jb20=
†These authors have contributed equally to this work and share first authorship