ObjectiveTo explore the molecular mechanisms by which serum containing Gegen Qinlian Tang （GQT） inhibits glycolysis and enhances chemotherapy sensitivity in 5-fluorouracil （5-FU）-resistant colorectal cancer （CRC） cells based on the cyclin-dependent kinase 16 （CDK16）/MYC proto-oncogene （MYC） pathway. MethodsHCT-116/5-FU cells were treated with different concentrations （5%， 10%， 20%， 30%） of GQT-containing serum. Cell viability and 5-FU sensitivity were assessed using the cell counting kit-8 （CCK-8） assay， and the experimental concentrations of 5-FU and GQT for subsequent experiments were determined. Cell proliferation and apoptosis under individual 5-FU， GQT， and combined 5-FU + GQT treatments were evaluated using 5-ethynyl-2′-deoxyuridine （EDU） staining and annexin V-FITC/PI double staining， respectively. Glucose consumption， adenosine triphosphate （ATP） production， and lactate levels were measured by colorimetric assays. Expression levels of glycolysis-related proteins， CDK16， MYC， and phosphorylated MYC were detected by Western blot. Co-immunoprecipitation （CoIP） was used to examine the protein interaction between CDK16 and MYC， and cycloheximide （CHX） treatment was applied to assess the effect of CDK16 overexpression on MYC protein stability. ResultsCCK-8 assays showed that 2.5 mg·L^-1 5-FU significantly inhibited HCT-116 cell viability in a dose-dependent manner. In HCT-116/5-FU cells， significant inhibition was observed only at 5 mg·L^-1 5-FU （P<0.05）， which was used for model establishment. Compared with 5-FU alone， addition of 5% GQT-containing serum significantly suppressed HCT-116/5-FU cell viability （P<0.05）， with stronger inhibition at higher serum concentrations. Thus， 5% GQT-containing serum was used in subsequent experiments. Compared with the control group， 5-FU， GQT， and 5-FU + GQT treatments all significantly reduced cell proliferation （P<0.05） and increased apoptosis （P<0.01）. The 5-FU + GQT combination showed superior inhibition of proliferation compared with 5-FU or GQT alone （P<0.01）， accompanied by more pronounced reductions in glucose consumption， ATP production， and lactate generation （P<0.01）. Additionally， compared with control， 5-FU， and GQT groups， the 5-FU + GQT group exhibited stronger suppression of MYC and its phosphorylated forms （P<0.01） and greater inhibition of glycolytic enzymes， including hexokinase 2 （HK2）， 3-phosphoinositide-dependent protein kinase 1 （PDK1）， lactate dehydrogenase A （LDHA）， and pyruvate kinase M2 （PKM2） （P<0.01）. CDK16， MYC， and MYC phosphorylation expression levels were significantly downregulated in the 5-FU + GQT group compared with the 5-FU group （all P<0.01）. MYC protein stability decreased in a time-dependent manner in the 5-FU + GQT group （P<0.05）， which was rescued by CDK16 overexpression （P<0.05）. ConclusionGQT significantly enhances the sensitivity of HCT-116/5-FU cells to 5-FU， potentially by inhibiting CDK16 and thereby reducing MYC-mediated glycolysis.

ObjectiveTo explore the possible mechanisms by which ligustilide （LIG） exerts neuroprotective effects on ischemic stroke （IS） by inhibiting the release of neutrophil extracellular traps （NETs）， promoting blood-brain barrier repair， and alleviating post-ischemic neuroinflammation， thereby providing a new direction for IS treatment. MethodsA middle cerebral artery occlusion （MCAO） model was established in rats. The rats were divided into the sham operation （Sham） group， model （Model） group， low- and high-dose LIG groups （20， 40 mg·kg^-1）， and the NET inhibitor CI-amidine group （CI-amidine， 10 mg·kg^-1）. Drug treatments were administered for 3 days. Neurological injury after ischemia was evaluated by 2，3，5-triphenyltetrazolium chloride （TTC） staining， neurological deficit scoring， and brain index measurement. Flow cytometry and Western blot were used to analyze changes in neutrophil expression. Immunofluorescence was used to observe the fluorescence intensity of the NET marker citrullinated histone H3 （H3Cit）. Western blot was performed to detect the expression of blood-brain barrier tight junction-related proteins and inflammatory factors， including interleukin-18 （IL-18） and interleukin-1β （IL-1β）. ResultsCompared with the Sham group， the Model group exhibited significant brain tissue injury （P<0.05）， significantly increased neutrophil numbers and NET expression （P<0.05）， significantly impaired blood-brain barrier permeability （P<0.05）， and significantly increased expression of inflammatory factors （P<0.05）. Compared with the Model group， both low- and high-dose LIG significantly alleviated brain tissue injury in rats （P<0.01）， inhibited neutrophil numbers and NET expression （P<0.01）， reduced blood-brain barrier damage （P<0.01）， and suppressed the expression of inflammatory factors IL-18 and IL-1β （P<0.01）， thereby ultimately exerting a neuroprotective effect. ConclusionThe neuroprotective effect of LIG in rats with cerebral ischemia-reperfusion injury may be related to inhibition of neutrophils and the NETs induced by them.

ObjectiveThis paper aims to find the identified and validated clinical biomarker data building upon a clinical study of early-phase phase Ⅱ and investigate the correlation analysis of Huanglian Jiedu Wan on syndrome improvement and clinical biomarkers in the treatment of "excess heat-toxicity" based on a machine learning model. Additionally， the effective prediction of clinical biomarker values for the main symptoms of the "excess heat-toxicity" syndrome was assessed. MethodsA total of 229 patients meeting the inclusion criteria for "excess heat-toxicity" syndrome were randomly divided into the Huanglian Jiedu Wan group and the placebo group. Syndrome score transition matrices were constructed for the Huanglian Jiedu Wan group and the placebo group based on three main symptoms of "excess heat-toxicity" syndrome， such as oral ulcers， sore throat， and gum swelling and pain. Data from the patients with these three syndromes were also integrated for an overall analysis. The corresponding syndrome score transition matrices were further constructed to visualize symptom change trends of the patients in the two groups via heatmaps. Based on the identified and validated clinical biomarkers related to inflammation， oxidative stress， and energy metabolism in the early phase， Spearman correlation analysis was employed to analyze and evaluate the associations between clinical biomarkers and syndrome improvement. Key clinical biomarkers reflecting the effect of Huanglian Jiedu Wan were screened through the comparison of differences between groups. An extreme gradient boosting （XGBoost） algorithm was used to develop a prediction model for main symptom classification， with classification performance evaluated through 10-fold cross-validation. Feature importance analysis was applied to identify variables with the greatest contribution to the prediction result. ResultsThe syndrome transition matrix results indicated that the Huanglian Jiedu Wan group showed a superior effect to the placebo group in improving oral ulcers， sore throat， and overall symptoms， with significant effects observed especially in sore throat and overall symptom analyses （P<0.01）. Spearman correlation analysis revealed that several clinical biomarkers positively correlated with "excess heat-toxicity" syndrome and its main symptom improvement， were also called "heat-related biomarkers"， including succinic acid， α-ketoglutaric acid， glycine， lactic acid， adenosine monophosphate （AMP）， tumor necrosis factor-α （TNF-α）， interferon-γ （IFN-γ）， interleukin-1β （IL-1β）， interleukin-4 （IL-4）， interleukin-6 （IL-6）， interleukin-8 （IL-8）， interleukin-10 （IL-10）， and so on. Conversely， clinical biomarkers negatively correlated with symptom severity， were also called "heat-clearing related biomarkers" after administration of Huanglian Jiedu Wan， including malic acid， fumaric acid， cis-aconitic acid， adrenocorticotropic hormone （ACTH）， IL-1β， IL-4， IL-8， succinic acid， and citric acid. The XGBoost classification model using all 52 biomarkers as variables achieved an average test accuracy of 0.754 and an average F1 score of 0.777. Feature importance analysis identified the scores of glutamic acid in saliva and IL-6 were the highest in all the variables， with importance scores of 0.081 and 0.080， respectively. After screening out 14 key variables and optimizing the parameters， model performance improved to an average accuracy of 0.758 and an F1 score of 0.798. Feature importance analysis further determined that the glutamic acid in saliva and IL-6 showed obvious changes after screening the variables， confirming the good syndrome prediction ability of the model constructed by these key clinical biomarkers. ConclusionThis study systematically elucidates the correlation between syndrome improvement and clinical biomarkers of Huanglian Jiedu Wan in the treatment of "excess heat-toxicity" syndrome. An XGBoost classification model based on key clinical biomarkers is successfully established， achieving effective prediction of the symptoms related to the "excess heat-toxicity" syndrome such as oral ulcers and sore throat and providing a new insight for objective identification of traditional Chinese medicine syndromes.

ObjectiveGastrodia elata has evolved ecological types with shortened rhizome internodes and diversified flower and fruit coloration in response to different altitudes. Studying the genetic mechanisms of different ecotype germplasm is significant for guiding variety breeding in different cultivation areas. MethodsThe bHLH gene family was identified based on the whole-genome datasets of G. elata f. elata and G. elata f. glauca. Subsequently， the gene family members were subject to analysis， including gene structure， chromosomal localization， cis-acting elements， gene synteny， and phylogeny. Combined with transcriptome data and quantitative Real-time PCR， the expression patterns of bHLH genes in the stems of the different G. elata ecotype germplasm were analyzed. Finally， correlation analysis was conducted between gene expression patterns and color to obtain the key bHLH genes regulating the color formation of stem. ResultsA total of 63 bHLH genes were identified in both G elata f. elata and G. elata f. glauca， unevenly distributed across 17 chromosomes and clustered into 16 subfamilies， with significant expansion in some family members. Obvious inversions of bHLH genes on the same chromosome and interchromosomal translocations were detected in the two ecotype germplasm. Among these genes， 12 bHLH genes （such as bHLH62-3 and bHLH74） were associated with the bright yellow color of G elata f. elata stem， while 9 bHLH genes （such as PIL13， UNE12， and bHLH130） were correlated with the red color of G. elata f. glauca stem. Compared to G. elata f. glauca， the bHLH48 expression level was significantly higher in flowers and scale leaves of G elata f. elata， and the bHLH62-3 expression level was significantly higher in all organs of G elata f. elata. ConclusionsFunctional pathway divergence of the bHLH family members has occurred across different chromosomes in G elata f. elata and G. elata f. glauca. Through synergism or antagonism with other genes， 21 bHLH genes participate in the coloration metabolic pathway regulation of stems， flowers， and fruits. Specifically， bHLH62-3 is involved in regulating stem color differentiation in the anthocyanin biosynthesis pathway of G. elata， thus relevant to the color formation of stem. Additionally， GebHLH48 positively regulates flowering-related pathways to promote the early-flowering phenotype of G. elata f. elata. These findings have laid the foundation for analyzing the genetic regulatory mechanisms underlying the color formation of the G. elata stem.

Chronic obstructive pulmonary disease（COPD） is a common chronic respiratory disorder frequently accompanied by diaphragmatic dysfunction during its course， which significantly increases respiratory burden and impairs quality of life. As the primary inspiratory muscle， the diaphragm is prone to fatigue， atrophy， inflammation， and fibrosis during the long-term progression of COPD. Its pathological mechanisms involve multiple pathways such as inflammatory responses， oxidative stress， mitochondrial dysfunction， apoptosis， ion channel abnormalities， epigenetic regulation， autophagy disorder， and protein metabolism imbalance. In recent years， traditional Chinese medicine（TCM） has demonstrated multi-targeted and systemic regulatory advantages in improving diaphragmatic function in COPD. However， related studies remain fragmented， and integrated mechanistic understanding is lacking. This paper focuses on the mechanism-target-TCM intervention framework， systematically summarizing the molecular mechanisms of diaphragmatic dysfunction， while incorporating the TCM theory of Zongqi（ancestral Qi）. It highlights the therapeutic effects of Chinese herbal formulas， single herbs， and active components in modulating inflammation， oxidative stress， mitochondrial function， ion channels， epigenetic processes， autophagy， and protein homeostasis. Additionally， the review outlines existing challenges， including insufficient study volume， unbalanced selection of herbal prescriptions， limited mechanistic depth， inconsistent disease models and experimental designs， lack of standardized diaphragmatic function assessment， and weak clinical validation. Future research should strengthen the integration of TCM and modern medicine， identify additional therapeutic targets， deepen mechanistic research， and establish unified and standardized experimental systems to advance the theoretical foundation and clinical application of TCM in the prevention and treatment of COPD-related diaphragmatic dysfunction.

Background Lung cancer, one of the most common malignant tumors worldwide, has long ranked first in cancer incidence and mortality, posing a severe challenge to public health systems. Objective To analyze the trends in incidence, mortality, and disability-adjusted life years (DALYs) of lung cancer among residents in Jing'an District, Shanghai, from 2002 to 2021, explore the impacts of population aging, population growth, and age-specific prevalence on disease burden, and provide a scientific basis for optimizing regional lung cancer prevention and control strategies. Methods Based on the cancer registration and cause-of-death surveillance data of registered residents in Jing'an District, Shanghai, from 2002 to 2021, Joinpoint regression models were used to analyze the annual change trends (APC) and average annual change trends (AAPC) of lung cancer incidence, mortality, DALY rate, and their age-standardized rates. Decomposition analysis was applied to quantify the contribution of population aging, population growth, and age-specific prevalence to changes in the number of new cases, deaths, and DALYs. Results From 2002 to 2021, the crude incidence rate of lung cancer in Jing'an District increased from 68.00 per 100000 to 144.36 per 100000 (AAPC=4.86%, P<0.001), and the age-standardized incidence rate rose from 43.09 per 100000 to 61.46 per 100000 (AAPC=2.89%, P<0.001). The growth rate of incidence was significantly higher in females (AAPC=6.73%) than in males (AAPC=3.62%). The crude mortality rate increased from 58.08 per 100000 to 66.11 per 100000 (AAPC=1.10%, P<0.001), while the age-standardized mortality rate decreased (AAPC=−1.54%, P<0.001). The crude DALY rate showed an upward trend (AAPC=0.94%, P=0.002), whereas the age-standardized DALY rate declined (AAPC=−1.63%, P<0.001). Significant differences were observed across genders and age groups: the mortality and DALY rates increased in males aged 60-64 years, while the fastest incidence growth occurred in females aged 45-64 years. The results of decomposition analysis indicated that the increase in new cases was primarily attributed to age-specific incidence (38.15%) and aging (37.31%); the growth in deaths and DALYs was driven by aging, while age-specific mortality and population growth showed offsetting effects. From 2002 to 2021, the probability of premature death from lung cancer among males (2.04%-2.52%) in Jing'an District showed no statistically significant trend (AAPC=−0.28%, P=0.440), while it presented a decreasing trend in the total population (1.37%-1.89%) and females (0.72%-1.26%) (AAPC=−0.80% and −2.23%, respectively; P=0.029 and <0.001, respectively). Conclusion From 2002 to 2021, the risk of lung cancer incidence in Jing’an District, Shanghai, continue to increase, with a disease burden exhibiting significant heterogeneity by gender and age: key features include a rapid increase in new cases and a trend toward younger onset among females, alongside rising mortality risk among males aged 60–64 years. Population aging and increasing age-specific incidence rates are primary drivers of the growing number of cases, while declining age-specific mortality rates mitigate the rise in mortality burden.

ObjectiveTo explore the molecular mechanisms by which serum containing Gegen Qinlian Tang （GQT） inhibits glycolysis and enhances chemotherapy sensitivity in 5-fluorouracil （5-FU）-resistant colorectal cancer （CRC） cells based on the cyclin-dependent kinase 16 （CDK16）/MYC proto-oncogene （MYC） pathway. MethodsHCT-116/5-FU cells were treated with different concentrations （5%， 10%， 20%， 30%） of GQT-containing serum. Cell viability and 5-FU sensitivity were assessed using the cell counting kit-8 （CCK-8） assay， and the experimental concentrations of 5-FU and GQT for subsequent experiments were determined. Cell proliferation and apoptosis under individual 5-FU， GQT， and combined 5-FU + GQT treatments were evaluated using 5-ethynyl-2′-deoxyuridine （EDU） staining and annexin V-FITC/PI double staining， respectively. Glucose consumption， adenosine triphosphate （ATP） production， and lactate levels were measured by colorimetric assays. Expression levels of glycolysis-related proteins， CDK16， MYC， and phosphorylated MYC were detected by Western blot. Co-immunoprecipitation （CoIP） was used to examine the protein interaction between CDK16 and MYC， and cycloheximide （CHX） treatment was applied to assess the effect of CDK16 overexpression on MYC protein stability. ResultsCCK-8 assays showed that 2.5 mg·L^-1 5-FU significantly inhibited HCT-116 cell viability in a dose-dependent manner. In HCT-116/5-FU cells， significant inhibition was observed only at 5 mg·L^-1 5-FU （P<0.05）， which was used for model establishment. Compared with 5-FU alone， addition of 5% GQT-containing serum significantly suppressed HCT-116/5-FU cell viability （P<0.05）， with stronger inhibition at higher serum concentrations. Thus， 5% GQT-containing serum was used in subsequent experiments. Compared with the control group， 5-FU， GQT， and 5-FU + GQT treatments all significantly reduced cell proliferation （P<0.05） and increased apoptosis （P<0.01）. The 5-FU + GQT combination showed superior inhibition of proliferation compared with 5-FU or GQT alone （P<0.01）， accompanied by more pronounced reductions in glucose consumption， ATP production， and lactate generation （P<0.01）. Additionally， compared with control， 5-FU， and GQT groups， the 5-FU + GQT group exhibited stronger suppression of MYC and its phosphorylated forms （P<0.01） and greater inhibition of glycolytic enzymes， including hexokinase 2 （HK2）， 3-phosphoinositide-dependent protein kinase 1 （PDK1）， lactate dehydrogenase A （LDHA）， and pyruvate kinase M2 （PKM2） （P<0.01）. CDK16， MYC， and MYC phosphorylation expression levels were significantly downregulated in the 5-FU + GQT group compared with the 5-FU group （all P<0.01）. MYC protein stability decreased in a time-dependent manner in the 5-FU + GQT group （P<0.05）， which was rescued by CDK16 overexpression （P<0.05）. ConclusionGQT significantly enhances the sensitivity of HCT-116/5-FU cells to 5-FU， potentially by inhibiting CDK16 and thereby reducing MYC-mediated glycolysis.

ObjectiveTo explore the possible mechanisms by which ligustilide （LIG） exerts neuroprotective effects on ischemic stroke （IS） by inhibiting the release of neutrophil extracellular traps （NETs）， promoting blood-brain barrier repair， and alleviating post-ischemic neuroinflammation， thereby providing a new direction for IS treatment. MethodsA middle cerebral artery occlusion （MCAO） model was established in rats. The rats were divided into the sham operation （Sham） group， model （Model） group， low- and high-dose LIG groups （20， 40 mg·kg^-1）， and the NET inhibitor CI-amidine group （CI-amidine， 10 mg·kg^-1）. Drug treatments were administered for 3 days. Neurological injury after ischemia was evaluated by 2，3，5-triphenyltetrazolium chloride （TTC） staining， neurological deficit scoring， and brain index measurement. Flow cytometry and Western blot were used to analyze changes in neutrophil expression. Immunofluorescence was used to observe the fluorescence intensity of the NET marker citrullinated histone H3 （H3Cit）. Western blot was performed to detect the expression of blood-brain barrier tight junction-related proteins and inflammatory factors， including interleukin-18 （IL-18） and interleukin-1β （IL-1β）. ResultsCompared with the Sham group， the Model group exhibited significant brain tissue injury （P<0.05）， significantly increased neutrophil numbers and NET expression （P<0.05）， significantly impaired blood-brain barrier permeability （P<0.05）， and significantly increased expression of inflammatory factors （P<0.05）. Compared with the Model group， both low- and high-dose LIG significantly alleviated brain tissue injury in rats （P<0.01）， inhibited neutrophil numbers and NET expression （P<0.01）， reduced blood-brain barrier damage （P<0.01）， and suppressed the expression of inflammatory factors IL-18 and IL-1β （P<0.01）， thereby ultimately exerting a neuroprotective effect. ConclusionThe neuroprotective effect of LIG in rats with cerebral ischemia-reperfusion injury may be related to inhibition of neutrophils and the NETs induced by them.

ObjectiveThis paper aims to find the identified and validated clinical biomarker data building upon a clinical study of early-phase phase Ⅱ and investigate the correlation analysis of Huanglian Jiedu Wan on syndrome improvement and clinical biomarkers in the treatment of "excess heat-toxicity" based on a machine learning model. Additionally， the effective prediction of clinical biomarker values for the main symptoms of the "excess heat-toxicity" syndrome was assessed. MethodsA total of 229 patients meeting the inclusion criteria for "excess heat-toxicity" syndrome were randomly divided into the Huanglian Jiedu Wan group and the placebo group. Syndrome score transition matrices were constructed for the Huanglian Jiedu Wan group and the placebo group based on three main symptoms of "excess heat-toxicity" syndrome， such as oral ulcers， sore throat， and gum swelling and pain. Data from the patients with these three syndromes were also integrated for an overall analysis. The corresponding syndrome score transition matrices were further constructed to visualize symptom change trends of the patients in the two groups via heatmaps. Based on the identified and validated clinical biomarkers related to inflammation， oxidative stress， and energy metabolism in the early phase， Spearman correlation analysis was employed to analyze and evaluate the associations between clinical biomarkers and syndrome improvement. Key clinical biomarkers reflecting the effect of Huanglian Jiedu Wan were screened through the comparison of differences between groups. An extreme gradient boosting （XGBoost） algorithm was used to develop a prediction model for main symptom classification， with classification performance evaluated through 10-fold cross-validation. Feature importance analysis was applied to identify variables with the greatest contribution to the prediction result. ResultsThe syndrome transition matrix results indicated that the Huanglian Jiedu Wan group showed a superior effect to the placebo group in improving oral ulcers， sore throat， and overall symptoms， with significant effects observed especially in sore throat and overall symptom analyses （P<0.01）. Spearman correlation analysis revealed that several clinical biomarkers positively correlated with "excess heat-toxicity" syndrome and its main symptom improvement， were also called "heat-related biomarkers"， including succinic acid， α-ketoglutaric acid， glycine， lactic acid， adenosine monophosphate （AMP）， tumor necrosis factor-α （TNF-α）， interferon-γ （IFN-γ）， interleukin-1β （IL-1β）， interleukin-4 （IL-4）， interleukin-6 （IL-6）， interleukin-8 （IL-8）， interleukin-10 （IL-10）， and so on. Conversely， clinical biomarkers negatively correlated with symptom severity， were also called "heat-clearing related biomarkers" after administration of Huanglian Jiedu Wan， including malic acid， fumaric acid， cis-aconitic acid， adrenocorticotropic hormone （ACTH）， IL-1β， IL-4， IL-8， succinic acid， and citric acid. The XGBoost classification model using all 52 biomarkers as variables achieved an average test accuracy of 0.754 and an average F1 score of 0.777. Feature importance analysis identified the scores of glutamic acid in saliva and IL-6 were the highest in all the variables， with importance scores of 0.081 and 0.080， respectively. After screening out 14 key variables and optimizing the parameters， model performance improved to an average accuracy of 0.758 and an F1 score of 0.798. Feature importance analysis further determined that the glutamic acid in saliva and IL-6 showed obvious changes after screening the variables， confirming the good syndrome prediction ability of the model constructed by these key clinical biomarkers. ConclusionThis study systematically elucidates the correlation between syndrome improvement and clinical biomarkers of Huanglian Jiedu Wan in the treatment of "excess heat-toxicity" syndrome. An XGBoost classification model based on key clinical biomarkers is successfully established， achieving effective prediction of the symptoms related to the "excess heat-toxicity" syndrome such as oral ulcers and sore throat and providing a new insight for objective identification of traditional Chinese medicine syndromes.

ObjectiveGastrodia elata has evolved ecological types with shortened rhizome internodes and diversified flower and fruit coloration in response to different altitudes. Studying the genetic mechanisms of different ecotype germplasm is significant for guiding variety breeding in different cultivation areas. MethodsThe bHLH gene family was identified based on the whole-genome datasets of G. elata f. elata and G. elata f. glauca. Subsequently， the gene family members were subject to analysis， including gene structure， chromosomal localization， cis-acting elements， gene synteny， and phylogeny. Combined with transcriptome data and quantitative Real-time PCR， the expression patterns of bHLH genes in the stems of the different G. elata ecotype germplasm were analyzed. Finally， correlation analysis was conducted between gene expression patterns and color to obtain the key bHLH genes regulating the color formation of stem. ResultsA total of 63 bHLH genes were identified in both G elata f. elata and G. elata f. glauca， unevenly distributed across 17 chromosomes and clustered into 16 subfamilies， with significant expansion in some family members. Obvious inversions of bHLH genes on the same chromosome and interchromosomal translocations were detected in the two ecotype germplasm. Among these genes， 12 bHLH genes （such as bHLH62-3 and bHLH74） were associated with the bright yellow color of G elata f. elata stem， while 9 bHLH genes （such as PIL13， UNE12， and bHLH130） were correlated with the red color of G. elata f. glauca stem. Compared to G. elata f. glauca， the bHLH48 expression level was significantly higher in flowers and scale leaves of G elata f. elata， and the bHLH62-3 expression level was significantly higher in all organs of G elata f. elata. ConclusionsFunctional pathway divergence of the bHLH family members has occurred across different chromosomes in G elata f. elata and G. elata f. glauca. Through synergism or antagonism with other genes， 21 bHLH genes participate in the coloration metabolic pathway regulation of stems， flowers， and fruits. Specifically， bHLH62-3 is involved in regulating stem color differentiation in the anthocyanin biosynthesis pathway of G. elata， thus relevant to the color formation of stem. Additionally， GebHLH48 positively regulates flowering-related pathways to promote the early-flowering phenotype of G. elata f. elata. These findings have laid the foundation for analyzing the genetic regulatory mechanisms underlying the color formation of the G. elata stem.