Objective To examine the impact of extreme temperature and drought stress on the survival of Bulinus globosus, so as to provide the theoretical evidence for the genomic research of Bulinus in absence of reference genes. Methods B. globosus snail samples were collected from Kiwani Shehia in Pemba Island, Zanzibar, Tanzania, and offspring snails were obtained through laboratory breeding and reproduction. A total of 120 10-week-old B. globosus snails from the same generation were selected and randomly assigned into four groups, including the high-temperature drought (HD) group, normal temperature drought (D) group, low-temperature drought (LD) group, and the control (C) group, of 30 snails in each group. Snails in HD, D, and LD groups were placed in beakers containing dry soil at the bottom and subsequently housed in climate chambers at 35, 26 ℃, and 10 ℃, respectively, while snails in Group C were maintained in 500 mL petri dishes containing dechlorinated tap water at 26 ℃. Following 3 days of breeding, living snails in each group were collected, and soft tissues were dissected and isolated. Total RNA was extracted from snail soft tissues for library construction, followed by high-throughput sequencing on the Illumina HiSeq 4000 sequencing system. De novo transcriptome assembly was performed using the Trinity software, and the longest transcripts were selected as unigenes. Gene functional annotations of unigenes were conducted using the Diamond software against Gene Ontology (GO) knowledgebase, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database, NCBI non-redundant (NR) protein sequences database, Protein Family (Pfam) database, and UniProtKB/Swiss-Prot (Swiss-Prot) knowledgebase. GO and KEGG enrichment analyses of differentially expressed genes (DEGs) were performed using the topGO and clusterProfiler software, respectively. In addition, four relevant genes were selected for validation using a real-time quantitative PCR (qRT-PCR) assay to verify the reliability of transcriptome sequencing results. Results Following 3 days of breeding, there were 7, 20, 28, and 30 survival B. globosus snails in HD, LD, D, and C groups, with corresponding survival rates of 23.33% (7/30), 66.67% (20/30), 93.33% (28/30), and 100.00% (30/30), respectively (χ² = 52.72, P < 0.001). De novo transcriptome assembly generated 176 942 unigenes, with annotation rates of 0.98%, 13.49%, 26.46%, 12.48%, and 14.39% against GO knowledgebase, KEGG pathway database, NR protein sequences database, Pfam database, and Swiss-Prot knowledgebase, respectively. There were 33 up-regulated and 72 down-regulated genes in Group D, 483 up-regulated and 815 down-regulated genes in Group HD, and 245 up-regulated and 172 down-regulated genes in Group LD relative to in Group C. Following removal of overlapping genes across groups and unmatched genes, 11 candidate genes were identified. GO and KEGG analyses revealed 3 heat shock protein (HSP)-related DEGs in these 11 candidate genes, which were annotated as HSP12.2, HSP70, and HSP20 genes and were all significantly up-regulated in each treatment group. Three immune and nervous system-related DEGs were identified, and were all significantly down-regulated in each treatment group, which were involved in the neural cell adhesion molecule L1-like protein pathway, fibrinogen binding protein pathway, and leukocyte elastase inhibitor-like protein pathway. qRT-PCR assay quantified that the expression trends of four genes related to temperature and drought stress across different treatment groups were highly consistent with transcriptome sequencing data. Conclusion The survival rate of B. globosus significantly reduces under combined stresses of extreme temperature and drought, possibly due to an imbalance in its cellular homeostasis regulatory system.

ObjectiveTo explore the mechanism by which Banxia Xiexin Tang （BXT） regulates the advanced glycation end products （AGE）/receptor for advanced glycation end products （RAGE） signaling pathway to reduce neuroinflammatory responses and ameliorate cognitive dysfunction in the rat model of vascular dementia （VD）. MethodsThe components of BXT were detected by ultra performance liquid chromatography-quadrupole -orbitrap-tandem mass spectrometry（UPLC-Q-Orbitrap-MS/MS）， and the core components and key action pathways were screened out by network pharmacology and molecular docking. Sixty SPF-grade male SD rats were randomly allocated into the sham and modeling groups by the random number table method. The VD model was replicated by the modified bilateral occlusion of the common carotid arteries （2-VO） method. The successfully modeled rats were randomly allocated into the model， low-， medium-， and high-dose （3.748 5， 7.497， 14.994 g·kg^-1） BXT （BXT-L， BXT-M， and BXT-H）， and nimodipine （NMP， 0.002 7 g·kg^-1） groups according to the random number table method. The rats in the drug intervention groups were administrated with corresponding drugs by gavage， and the sham and model groups received the same amount of normal saline for 14 consecutive days. The Morris water maze， Y-maze， and new object recognition experiments were conducted to evaluate the cognitive dysfunction of rats. Hematoxylin-eosin （HE） staining was used to evaluate the histopathological changes of the hippocampal tissue in rats. The mRNA levels of AGE， RAGE， and phosphorylated nuclear factor-kappa B p65 （p-NF-κB p65） in the hippocampal tissue of rats were determined by Real-time fluorescence quantitative polymerase chain reaction （Real-time PCR）. The expression of related proteins in the AGE/RAGE pathway in the hippocampal tissue of rats was determined by Western blot and immunohistochemistry （IHC）. The levels of neurotransmitters and inflammatory mediators in the rat serum were measured by enzyme-linked immunosorbent assay （ELISA）. ResultsThe chemical components of BXT were detected by UPLC-Q-Orbitrap-MS/MS. Network pharmacology and molecular docking identified the AGE/RAGE pathway as the key pathway. The results of the water maze， Y maze， and novel object recognition tests showed that compared with the sham group， the model group demonstrated prolonged successful latency and decreases in number of platform crossings， alternation rate， number of entries into the new arm， preference index， and discrimination index （P0.01）. Compared with the model group， the BXT-H and BXT-M groups showed shortened successful latency （P0.01） and increases in number of platform crossings （P0.05）， alternation rate （P0.01）， number of entries into the new arm （P0.05）， preference index （P0.01）， and discrimination index （P0.01）. HE results showed that compared with the sham group， the cells of model rats were loosely and disorderly arranged， and the nuclei were condensed. Compared with the model group， the pathological changes of the hippocampus in the BXT group were mitigated. Real-time PCR results showed that compared with the sham group， the model group presented up-regulated mRNA levels of AGE， RAGE， and p-NF-κB p65 in the hippocampus （P0.01）， and compared with the model group， the BXT-H and BXT-M groups showcased down-regulated mRNA levels of AGE， RAGE， and p-NF-κB p65 （P0.01）. Western blot results showed that compared with the sham group， the model group presented up-regulated expression of AGE， RAGE， p-NF-κB p65， and tumor necrosis factor-α （TNF-α） （P0.05）， and compared with the model group， the BXT-H group presented down-regulated expression of AGE， RAGE， p-NF-κB p65， and TNF-α （P0.05）. IHC results showed that compared with the sham group， the model group had increased expression of RAGE （P0.01）， and compared with the model group， the BXT-H and BXT-M groups had reduced expression of RAGE （P0.01）. ELISA results showed that compared with the sham group， the model group exhibited elevated levels of TNF-α and Interleukin-1β （IL-1β） and declined levels of acetylcholine （ACh） and dopamine （DA） in the serum （P0.01）. Compared with the model group， the BXT-L， BXT-M， and BXT-H groups showed lowered levels of TNF-α and IL-1β in the serum （P0.05） and elevated levels of ACh and DA （P0.05）. ConclusionBXT may ameliorate cognitive dysfunction in the rat model of VD by down-regulating the AGE/RAGE signaling pathway， reducing neuroinflammatory responses， and regulating neurotransmitter levels.

ObjectiveTo explore the mechanism by which Yizhi Qingxin prescription improves mitochondrial dysfunction in Alzheimer's disease （AD） through regulating mitochondrial Ca²⁺ homeostasis and kinetic balance based on the protein kinase A （PKA）/calcineurin （CaN） signaling pathway. MethodsSixty three-month-old amyloid precursor protein （APP）/presenilin 1 （PS1） double transgenic mice were randomly divided into a model group， a donepezil group（0.65 mg·kg^-1）， a low-dose Yizhi Qingxin prescription group （YQF-L，2.6 g·kg^-1）， a medium-dose Yizhi Qingxin prescription group （YQF-M，5.2 g·kg^-1）， and a high-dose Yizhi Qingxin prescription group （YQF-H，10.4 g·kg^-1）， with 12 mice in each group. Twelve C57BL/6J mice with the same genetic background served as a normal group. Each treatment group received gavage administration daily， with the model and normal groups receiving equal volume of physiological saline. Intervention continued for 12 consecutive weeks. The learning and memory abilities of the mice were assessed using the novel object recognition （NOR） and Morris water maze （MWM） tests. Hematoxylin-eosin （HE）/Nissl staining was used to observe histopathological changes in the hippocampus. Transmission electron microscopy （TEM） was used to observe mitochondrial ultrastructure. Fluo-4 acetoxymethyl ester （Fluo-4 AM） Ca²⁺ probe was used to measure intracellular Ca²⁺ concentration in brain tissue. Western blot was used to determine the protein expression of PKA， CaN， sodium/calcium/lithium exchanger （NCLX）， mitochondrial calcium uniporter （MCU）， calmodulin （CaM）， dynamin-related protein 1 （Drp1）， and phosphorylated dynamin-related protein 1 （serine 637 site） ［p-Drp1（S637）］ in the hippocampus. Real-time quantitative polymerase chain reaction （Real-time PCR） was used to measure the expression of PKA， CaN， CaM， NCLX， MCU， and Drp1 mRNAs. ResultsCompared with those in the normal group， the recognition index （RI） of the model group decreased （P0.01）， and the number of crossings through the original platform area， the duration of stay in the target quadrant， and the distance were reduced （P0.01）. The protein expression of PKA， NCLX， and p-DRP1 （ser637） significantly decreased （P0.05）， and the mRNA expression of PKA and NCLX significantly decreased （P0.05）. The escape latency （EL） was prolonged （P0.05）， and the intracellular Ca²⁺ level significantly increased （P0.01）. The protein expression of CaN， CaM， MCU， and Drp1， as well as the mRNA expression of CaN， MCU， and Drp1， significantly increased （P0.05）. After intervention with Donepezil and Yizhi Qingxin prescription， compared with that in the model group， the RI of the treatment group significantly increased （P0.05）， and the number of crossings through the platform and the duration of stay in the target quadrant significantly increased （P0.05）. The protein expression of PKA， NCLX， and p-Drp1 （ser637） and the mRNA expression of PKA and NCLX significantly increased （P0.05）. On the 4th and 5th days， the EL was shortened （P0.05）， and the intracellular Ca²⁺ level decreased （P0.05）. The protein expression of CaN， CaM， MCU， and Drp1 and the mRNA expression of CaN， MCU， and Drp1 significantly decreased （P0.05）. ConclusionYizhi Qingxin prescription regulates the PKA/CaN pathway， upregulates the expression of PKA， NCLX， and p-Drp1 （ser637） proteins， reduces the expression of CaN， CaM， MCU， and Drp1 proteins， and regulates Ca²⁺ homeostasis and mitochondrial dynamic balance， thereby enhancing the spatial learning and memory abilities of AD mice.

Insomnia is a sleep disorder characterized by difficulty in falling asleep， sleep maintenance disorder and impaired daytime function. Its pathological mechanism involves multiple factors such as nerve excitability， circadian rhythm， cell apoptosis， oxidative stress injury. As a classical tyrosine kinase signaling pathway， phosphatidylinositol 3-kinase/protein kinase B（PI3K/Akt） triggers Akt phosphorylation cascade， inducing inflammatory response， apoptosis， autophagy， oxidative damage， nerve excitability， and circadian rhythm imbalance. Traditional Chinese medicine（TCM） can improve sleep by targeting the PI3K/Akt pathway. Based on this， this paper systematically reviews the research progress on the regulation of PI3K/Akt pathway by traditional Chinese medicine（TCM） for insomnia at home and abroad. These drugs can regulate neuronal excitability by regulating the PI3K/Akt pathway， affect the circadian rhythm， alleviate inflammation， apoptosis， autophagy and oxidative stress， and thus regulate sleep-wake. Furthermore， literature review indicates that the PI3K/Akt signaling pathway may represent a specific pathway underlying phlegm-turbidity disturbing the upper Jiao-type insomnia.

ObjectiveTo explore the mechanism by which Banxia Xiexin Tang （BXT） regulates the advanced glycation end products （AGE）/receptor for advanced glycation end products （RAGE） signaling pathway to reduce neuroinflammatory responses and ameliorate cognitive dysfunction in the rat model of vascular dementia （VD）. MethodsThe components of BXT were detected by ultra performance liquid chromatography-quadrupole -orbitrap-tandem mass spectrometry（UPLC-Q-Orbitrap-MS/MS）， and the core components and key action pathways were screened out by network pharmacology and molecular docking. Sixty SPF-grade male SD rats were randomly allocated into the sham and modeling groups by the random number table method. The VD model was replicated by the modified bilateral occlusion of the common carotid arteries （2-VO） method. The successfully modeled rats were randomly allocated into the model， low-， medium-， and high-dose （3.748 5， 7.497， 14.994 g·kg^-1） BXT （BXT-L， BXT-M， and BXT-H）， and nimodipine （NMP， 0.002 7 g·kg^-1） groups according to the random number table method. The rats in the drug intervention groups were administrated with corresponding drugs by gavage， and the sham and model groups received the same amount of normal saline for 14 consecutive days. The Morris water maze， Y-maze， and new object recognition experiments were conducted to evaluate the cognitive dysfunction of rats. Hematoxylin-eosin （HE） staining was used to evaluate the histopathological changes of the hippocampal tissue in rats. The mRNA levels of AGE， RAGE， and phosphorylated nuclear factor-kappa B p65 （p-NF-κB p65） in the hippocampal tissue of rats were determined by Real-time fluorescence quantitative polymerase chain reaction （Real-time PCR）. The expression of related proteins in the AGE/RAGE pathway in the hippocampal tissue of rats was determined by Western blot and immunohistochemistry （IHC）. The levels of neurotransmitters and inflammatory mediators in the rat serum were measured by enzyme-linked immunosorbent assay （ELISA）. ResultsThe chemical components of BXT were detected by UPLC-Q-Orbitrap-MS/MS. Network pharmacology and molecular docking identified the AGE/RAGE pathway as the key pathway. The results of the water maze， Y maze， and novel object recognition tests showed that compared with the sham group， the model group demonstrated prolonged successful latency and decreases in number of platform crossings， alternation rate， number of entries into the new arm， preference index， and discrimination index （P<0.01）. Compared with the model group， the BXT-H and BXT-M groups showed shortened successful latency （P<0.01） and increases in number of platform crossings （P<0.05）， alternation rate （P<0.01）， number of entries into the new arm （P<0.05）， preference index （P<0.01）， and discrimination index （P<0.01）. HE results showed that compared with the sham group， the cells of model rats were loosely and disorderly arranged， and the nuclei were condensed. Compared with the model group， the pathological changes of the hippocampus in the BXT group were mitigated. Real-time PCR results showed that compared with the sham group， the model group presented up-regulated mRNA levels of AGE， RAGE， and p-NF-κB p65 in the hippocampus （P<0.01）， and compared with the model group， the BXT-H and BXT-M groups showcased down-regulated mRNA levels of AGE， RAGE， and p-NF-κB p65 （P<0.01）. Western blot results showed that compared with the sham group， the model group presented up-regulated expression of AGE， RAGE， p-NF-κB p65， and tumor necrosis factor-α （TNF-α） （P<0.05）， and compared with the model group， the BXT-H group presented down-regulated expression of AGE， RAGE， p-NF-κB p65， and TNF-α （P<0.05）. IHC results showed that compared with the sham group， the model group had increased expression of RAGE （P<0.01）， and compared with the model group， the BXT-H and BXT-M groups had reduced expression of RAGE （P<0.01）. ELISA results showed that compared with the sham group， the model group exhibited elevated levels of TNF-α and Interleukin-1β （IL-1β） and declined levels of acetylcholine （ACh） and dopamine （DA） in the serum （P<0.01）. Compared with the model group， the BXT-L， BXT-M， and BXT-H groups showed lowered levels of TNF-α and IL-1β in the serum （P<0.05） and elevated levels of ACh and DA （P<0.05）. ConclusionBXT may ameliorate cognitive dysfunction in the rat model of VD by down-regulating the AGE/RAGE signaling pathway， reducing neuroinflammatory responses， and regulating neurotransmitter levels.

ObjectiveTo explore the mechanism by which Yizhi Qingxin prescription improves mitochondrial dysfunction in Alzheimer's disease （AD） through regulating mitochondrial Ca²⁺ homeostasis and kinetic balance based on the protein kinase A （PKA）/calcineurin （CaN） signaling pathway. MethodsSixty three-month-old amyloid precursor protein （APP）/presenilin 1 （PS1） double transgenic mice were randomly divided into a model group， a donepezil group（0.65 mg·kg^-1）， a low-dose Yizhi Qingxin prescription group （YQF-L，2.6 g·kg^-1）， a medium-dose Yizhi Qingxin prescription group （YQF-M，5.2 g·kg^-1）， and a high-dose Yizhi Qingxin prescription group （YQF-H，10.4 g·kg^-1）， with 12 mice in each group. Twelve C57BL/6J mice with the same genetic background served as a normal group. Each treatment group received gavage administration daily， with the model and normal groups receiving equal volume of physiological saline. Intervention continued for 12 consecutive weeks. The learning and memory abilities of the mice were assessed using the novel object recognition （NOR） and Morris water maze （MWM） tests. Hematoxylin-eosin （HE）/Nissl staining was used to observe histopathological changes in the hippocampus. Transmission electron microscopy （TEM） was used to observe mitochondrial ultrastructure. Fluo-4 acetoxymethyl ester （Fluo-4 AM） Ca²⁺ probe was used to measure intracellular Ca²⁺ concentration in brain tissue. Western blot was used to determine the protein expression of PKA， CaN， sodium/calcium/lithium exchanger （NCLX）， mitochondrial calcium uniporter （MCU）， calmodulin （CaM）， dynamin-related protein 1 （Drp1）， and phosphorylated dynamin-related protein 1 （serine 637 site） ［p-Drp1（S637）］ in the hippocampus. Real-time quantitative polymerase chain reaction （Real-time PCR） was used to measure the expression of PKA， CaN， CaM， NCLX， MCU， and Drp1 mRNAs. ResultsCompared with those in the normal group， the recognition index （RI） of the model group decreased （P0.01）， and the number of crossings through the original platform area， the duration of stay in the target quadrant， and the distance were reduced （P0.01）. The protein expression of PKA， NCLX， and p-DRP1 （ser637） significantly decreased （P0.05）， and the mRNA expression of PKA and NCLX significantly decreased （P0.05）. The escape latency （EL） was prolonged （P0.05）， and the intracellular Ca²⁺ level significantly increased （P0.01）. The protein expression of CaN， CaM， MCU， and Drp1， as well as the mRNA expression of CaN， MCU， and Drp1， significantly increased （P0.05）. After intervention with Donepezil and Yizhi Qingxin prescription， compared with that in the model group， the RI of the treatment group significantly increased （P0.05）， and the number of crossings through the platform and the duration of stay in the target quadrant significantly increased （P0.05）. The protein expression of PKA， NCLX， and p-Drp1 （ser637） and the mRNA expression of PKA and NCLX significantly increased （P0.05）. On the 4th and 5th days， the EL was shortened （P0.05）， and the intracellular Ca²⁺ level decreased （P0.05）. The protein expression of CaN， CaM， MCU， and Drp1 and the mRNA expression of CaN， MCU， and Drp1 significantly decreased （P0.05）. ConclusionYizhi Qingxin prescription regulates the PKA/CaN pathway， upregulates the expression of PKA， NCLX， and p-Drp1 （ser637） proteins， reduces the expression of CaN， CaM， MCU， and Drp1 proteins， and regulates Ca²⁺ homeostasis and mitochondrial dynamic balance， thereby enhancing the spatial learning and memory abilities of AD mice.

Insomnia is a sleep disorder characterized by difficulty in falling asleep， sleep maintenance disorder and impaired daytime function. Its pathological mechanism involves multiple factors such as nerve excitability， circadian rhythm， cell apoptosis， oxidative stress injury. As a classical tyrosine kinase signaling pathway， phosphatidylinositol 3-kinase/protein kinase B（PI3K/Akt） triggers Akt phosphorylation cascade， inducing inflammatory response， apoptosis， autophagy， oxidative damage， nerve excitability， and circadian rhythm imbalance. Traditional Chinese medicine（TCM） can improve sleep by targeting the PI3K/Akt pathway. Based on this， this paper systematically reviews the research progress on the regulation of PI3K/Akt pathway by traditional Chinese medicine（TCM） for insomnia at home and abroad. These drugs can regulate neuronal excitability by regulating the PI3K/Akt pathway， affect the circadian rhythm， alleviate inflammation， apoptosis， autophagy and oxidative stress， and thus regulate sleep-wake. Furthermore， literature review indicates that the PI3K/Akt signaling pathway may represent a specific pathway underlying phlegm-turbidity disturbing the upper Jiao-type insomnia.

Objective: To investigate the current status of benefit finding among young and middle-aged patients with type 2 diabetes mellitus (T2DM) and analyze its influencing factors, so as to provide a reference for improving the level of benefit finding in this population. Methods: From November 2022 to May 2023, young and middle-aged patients with T2DM aged 18-59 years hospitalized in the endocrinology departments of 2 tertiary hospitals in Hengyang City, Hunan Province were selected as survey subjects by a convenience sampling method. Basic demographic information was collected using a general questionnaire survey. Benefit finding, resourcefulness, and stigma were evaluated using the Benefit Finding Scale, the Chinese Version of the Resourcefulness Scale, and the Type 2 Diabetes Stigma Assessment Scale, respectively. A multiple linear regression model was used to analyze the influencing factors of benefit finding among young and middle-aged patients with T2DM. Results: A total of 305 young and middle-aged patients with T2DM were investigated, including 222 males (72.79%) and 83 females (27.21%). There were 231 cases aged 45-59 years, accounting for 75.74%. The scores for benefit finding, resourcefulness, and stigma were (42.86±6.06), (75.12±11.30), and (41.20±10.10), respectively. Multiple linear regression analysis showed that young and middle-aged patients with T2DM who were male (β′=0.088), aged 18-<45 years (β′=0.083), absence of diabetes complications (β′=0.124), and had higher resourcefulness scores (β′=0.679) had higher levels of benefit finding, while patients with higher stigma scores (β′=-0.097) had lower levels of benefit finding. Conclusion The level of benefit finding among young and middle-aged patients with T2DM was moderate, and was related to gender, age, diabetes complications, resourcefulness, and stigma.

Objective: To investigate the status of heat illness, knowledge awareness and adaptive behaviors of heat wave knowledge among food delivery riders, so as to provide a basis for optimizing heat wave response measures for food delivery riders. Methods: In November 2022, food delivery riders from a large food delivery platform in Ningbo City, Zhejiang Province were selected as survey subjects using a cluster sampling method. A self-designed electronic questionnaire was used to select demographic information, work status, lifestyle behaviors and disease history, heat illness status, knowledge awareness and adaptive behaviors of heatwave. Results: A total of 911 questionnaires were distributed, and 830 valid questionnaires were recovered, resulting in a valid response rate of 91.11%. Among the respondents, 796 (95.90%) were male, and 818 (98.55%) worked full-time. The mean age was (27.75±8.00) years. A total of 470 respondents (56.63%) had a work tenure of less than 1 year. The primary working hours were 8-<12 hours, with 504 people accounting for 60.72%. There were 108 cases of heatstroke, with an occurrence rate of 13.01%. And 286 people reported heat-related symptoms, with an occurrence rate of 34.46%. The overall awareness rate of heat wave knowledge was 73.22%, while the awareness rate of heat warning signal classification was relatively low at 9.04%. The heat wave cognition score was (5.86±1.31) points. There were statistically significant differences in heat wave cognition scores among food delivery riders of different ages, educational levels, family annual income, work tenures, and work durations (all P<0.05). Regarding positive adaptive behaviors, the number of riders paying attention to weather forecasts and actively learning about preventive measures was higher (734 people each, accounting for 88.43%). Regarding negative adaptive behaviors, the number of riders who often drank ice-cold beverages was higher (509 people, accounting for 61.33%). The heat wave adaptive behavior score was (6.88±1.77) points. There were statistically significant differences in adaptive behavior scores among riders with different educational levels, family annual income, work tenures, and smoking frequency (all P<0.05). Conclusions The occurrence rates of heatstroke and heat-related symptoms among food delivery riders are relatively high. The knowledge awareness and adaptive behaviors regarding heat wave are at a moderate level. It is suggested to strengthen health education, reinforce risk cognition of heat wave, and promote positive adaptive behaviors among food delivery riders.

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.