Ischemic stroke （IS） represents a major global health challenge with complex pathological mechanisms. Although modern therapies such as intravenous thrombolysis and endovascular thrombectomy have advanced， their application remains constrained by narrow therapeutic time windows， hemorrhagic risks， and uneven distribution of medical resources. Traditional Chinese medicine （TCM） demonstrates unique value in the prevention and treatment of IS， owing to its multi-component， multi-target， and holistic regulatory characteristics. This review summarized the molecular mechanisms by which active ingredients and compound formulations of TCM exert therapeutic effects against IS through the regulation of inflammatory responses， oxidative stress， excitatory toxicity， apoptosis， and autophagy. Studies have indicated that components such as curcumin， baicalin， and astragaloside Ⅳ inhibit microglial activation and the nucleotide-binding oligomerization domain （NOD）-like receptor protein 3 （NLRP3） inflammasome to attenuate neuroinflammation， activate the nuclear factor erythroid 2-related factor 2（Nrf2）/heme oxygenase-1 （HO-1） pathway to alleviate oxidative stress， modulate glutamate receptor function to counteract excitatory toxicity， and regulate the B-cell lymphoma 2（Bcl-2）/Bcl-2-associated X protein （Bax）， cysteine aspartate-specific protease （Caspase）， and phosphatidylinositol 3 kinases （PI3K）/protein kinase B （Akt） signaling pathways to suppress neuronal apoptosis. Recent research has further revealed that TCM can modulate ferroptosis by targeting key proteins glutathione peroxidase 4 （GPX4） and acyl-coenzyme A synthetase long-chain family member 4 （ACSL4） to maintain iron homeostasis， intervene in the "microbiota-gut-brain axis" to ameliorate dysbiosis and reduce neuroinflammation， utilize exosomes for brain-targeted drug delivery， and influence neural repair processes through epigenetic regulation. Furthermore， the review discussed the integrated mechanisms of compound formulations， such as Buyang Huanwu Decoction， in improving cerebral microcirculation and promoting neurovascular remodeling via multi-component synergy. It also analyzed the strategy and advantages of integrating TCM with Western medicine for IS treatment， providing a novel theoretical foundation and research directions for future investigations and clinical translation of TCM in IS management.

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.

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.

Background Benign prostatic hyperplasia (BPH) is a common chronic urinary disease in middle-aged and older men, yet the impact of long-term exposure to atmospheric particulate matter (PM) on its pathogenesis remains unclear. Objective To investigate the association between PM exposure and the risk of incident BPH in middle-aged and older men. Methods Based on four waves of follow-up data (2011–2018) from the China Health and Retirement Longitudinal Study (CHARLS), 4766 participants were enrolled. Robust Poisson regression models were employed to assess the association between exposure to PM (PM₁, PM_2.5, and PM₁₀) and the risk of incident BPH. Relative risks (RR) and their corresponding 95% confidence intervals (95%CI) were calculated. Dose-response relationships were fitted using restricted cubic splines (RCS). Subgroup analyses were performed to explore potential effect modifications, and multiple imputation was used to handle missing data. Results Over a mean follow-up of 6 years, 914 incident BPH cases were identified among the4766 participants (cumulative incidence: 19.18%). After adjusting for confounders, each 10 μg·m⁻³ increase in PM₁, PM_2.5, and PM₁₀ concentrations was associated with a 13.1% (RR=1.131, 95%CI: 1.063, 1.203), 8.5% (RR=1.085, 95%CI: 1.050, 1.122), and 5.1% (RR=1.051, 95%CI: 1.034, 1.069) increased risk of BPH, respectively. RCS analysis showed that no nonlinear relationship was found between PM₁ and PM_2.5 and the risk of BPH (P>0.05); however, a nonlinear association was observed for PM₁₀ (P=0.03), with the risk increment slowing beyond 100 μg·m⁻³. Subgroup and sensitivity analyses confirmed the robustness of these findings. Conclusion Long-term exposure to ambient particulate matter may be associated with an increased risk of incident BPH in middle-aged and older men.

Background Benign prostatic hyperplasia (BPH) is a common chronic urinary disease in middle-aged and older men, yet the impact of long-term exposure to atmospheric particulate matter (PM) on its pathogenesis remains unclear. Objective To investigate the association between PM exposure and the risk of incident BPH in middle-aged and older men. Methods Based on four waves of follow-up data (2011–2018) from the China Health and Retirement Longitudinal Study (CHARLS), 4766 participants were enrolled. Robust Poisson regression models were employed to assess the association between exposure to PM (PM₁, PM_2.5, and PM₁₀) and the risk of incident BPH. Relative risks (RR) and their corresponding 95% confidence intervals (95%CI) were calculated. Dose-response relationships were fitted using restricted cubic splines (RCS). Subgroup analyses were performed to explore potential effect modifications, and multiple imputation was used to handle missing data. Results Over a mean follow-up of 6 years, 914 incident BPH cases were identified among the4766 participants (cumulative incidence: 19.18%). After adjusting for confounders, each 10 μg·m⁻³ increase in PM₁, PM_2.5, and PM₁₀ concentrations was associated with a 13.1% (RR=1.131, 95%CI: 1.063, 1.203), 8.5% (RR=1.085, 95%CI: 1.050, 1.122), and 5.1% (RR=1.051, 95%CI: 1.034, 1.069) increased risk of BPH, respectively. RCS analysis showed that no nonlinear relationship was found between PM₁ and PM_2.5 and the risk of BPH (P>0.05); however, a nonlinear association was observed for PM₁₀ (P=0.03), with the risk increment slowing beyond 100 μg·m⁻³. Subgroup and sensitivity analyses confirmed the robustness of these findings. Conclusion Long-term exposure to ambient particulate matter may be associated with an increased risk of incident BPH in middle-aged and older men.

ObjectiveTo investigate the role and underlying mechanism of the Wnt/β-catenin signaling pathway in the process of Xiaojinwan （XJW）-containing serum interfering with endothelial-mesenchymal transition （EMT） in human breast cancer MDA-MB-231 cells. MethodsXJW-containing serum was prepared by intragastric administration of XJW （4 g·kg^-1） to SD rats. MDA-MB-231 cells were cultured in vitro. β-Catenin overexpression and knockdown plasmids were constructed and transfected into MDA-MB-231 cells. The optimal XJW-containing serum （4 g·kg^-1） according to previous studies was used for β-catenin intervention and overexpression experiments. The groups for intervention/overexpression experiments included intervention/overexpression blank group， β-catenin intervention/overexpression group， intervention/overexpression blank + XJW group， and β-catenin intervention/overexpression + XJW group. Cell morphology was observed under a microscope. Flow cytometry was used to detect and record apoptosis rates. Western blot was performed to measure apoptosis-related proteins ［B-cell lymphoma 2 （Bcl-2）， Bcl-2-associated X protein （Bax）， and active-Caspase-3］ and EMT-related proteins （E-cadherin， N-cadherin， and vimentin）. Real-time polymerase chain reaction （Real-time PCR） was used to detect β-catenin mRNA expression levels in different groups. ResultsIntervention experiment results showed that compared with the intervention blank group， the β-catenin intervention group and the intervention blank + XJW group showed increased apoptosis rates （P<0.05， P<0.01）， elevated protein levels of Bax， active Caspase-3， and E-cadherin （P<0.05， P<0.01）， and decreased protein levels of Bcl-2 and Vimentin （P<0.05， P<0.01）. β-catenin mRNA expression was reduced in the β-catenin intervention group （P<0.05， 0.63 ± 0.15） and the intervention blank + XJW group （P<0.05， 0.79 ± 0.33）. Compared with the β-catenin intervention group， the β-catenin intervention + XJW group exhibited an increased apoptosis rate （P<0.01）， elevated Bax and active Caspase-3 protein levels （P<0.05， P<0.01）， and decreased Bcl-2， N-cadherin， and Vimentin protein levels （P<0.05）. β-catenin mRNA expression was further reduced in the β-catenin intervention + XJW group （P<0.01， 0.22 ± 0.08）. Compared with the intervention blank + XJW group， the β-catenin intervention + XJW group showed an increased apoptosis rate （P<0.01）， elevated Bax， active Caspase-3， and E-cadherin protein levels （P<0.05， P<0.01）， and decreased Bcl-2， N-cadherin， and vimentin protein levels （P<0.05， P<0.01）. β-catenin mRNA expression was also reduced （P<0.05）. Overexpression experiment results showed that compared with the overexpression blank group， the β-catenin overexpression group exhibited a reduced apoptosis rate （P<0.01）， while the overexpression blank + XJW group showed an increased apoptosis rate （P<0.05）. In the β-catenin overexpression and β-catenin overexpression + XJW groups， Bax and E-cadherin protein levels decreased （P<0.05， P<0.01）， while Bcl-2， N-cadherin， and vimentin levels increased （P<0.05）. In the overexpression blank + XJW group， Bax and E-cadherin levels increased （P<0.05， P<0.01）， while Bcl-2， N-cadherin， and vimentin levels decreased （P<0.05， P<0.01）. β-catenin mRNA expression increased in the β-catenin overexpression group （P<0.01， 7.26 ± 1.03） but decreased in the overexpression blank + XJW group （P<0.05， 0.73 ± 0.14）. Compared with the β-catenin overexpression group， the β-catenin overexpression + XJW group showed an increased apoptosis rate （P<0.01）， elevated Bax and E-cadherin levels （P<0.05）， and decreased N-cadherin and vimentin levels （P<0.05）. β-catenin mRNA expression was reduced （P<0.05， 4.65 ± 1.02）. Compared with the overexpression blank + XJW group， the β-catenin overexpression + XJW group exhibited a reduced apoptosis rate （P<0.01）， decreased Bax and E-cadherin levels （P<0.01）， and increased Bcl-2， N-cadherin， and Vimentin levels （P<0.05）. β-catenin mRNA expression increased （P<0.01）. ConclusionXJW-containing serum inhibits breast cancer cell proliferation by suppressing EMT， which is associated with the regulation of the Wnt/β-catenin pathway.

Objective To investigate the prevalence of spotted fever group rickettsiae (SFGR) infections in free-living ticks and among high-risk populations in the southern margin of the Gurbantunggut Desert, Xinjiang Uygur Autonomous Region, so as to provide insights into tick-borne diseases prevention and control in this region. Methods Free-living ticks were collected from five counties (cities) in the southern margin of the Gurbantunggut Desert, Xinjiang Uygur Autonomous Region. Tick species were identified by morphology combined with mitochondrial cytochrome c oxidase subunit Ⅰ (COI) gene analysis. Four target gene fragments of Rickettsia were amplified using a nested PCR assay, including outer membrane protein B (ompB), 17-kDa common surface antigen protein (17kDa), surface cell antigen 1 (sca1), and surface cell antigen 4 (sca4) genes, to detect the prevalence of SFGR infections in tick samples. Positive PCR amplification products were purified, sequenced, and subjected to sequence alignment, and phylogenetic trees were created using the neighbor-joining (NJ) method. Serum samples were collected from individuals at a high risk of tick exposure from six counties (cities, districts) in the southern margin of the Gurbantunggut Desert, Xinjiang Uygur Autonomous Region, and serum levels of anti-R. sibirica IgG antibodies were detected among these individuals using indirect immunofluorescence assay (IFA) and statistically analyzed. Results A total of 287 free-living ticks were collected, and all were identified as Hyalomma asiaticum. The overall prevalence of SFGR infection was 78.40% (225/287) in tick samples, and five Rickettsia species were characterized, including R. sibirica, which shared 99.81% (515/516) sequence homology with the strain detected in Dermacentor sinicus from China; R. massiliae, which shared 99.80% (509/510) sequence homology with the strain detected in Rhipicephalus sanguineus from Ethiopia; R. raoultii, which shared 100.00% (516/516) sequence homology with the strain detected in Hyalomma asiaticum from Xinjiang Uygur Autonomous Region in China; R. africae, which shared 100.00% (524/524) sequence homology with the strain detected in Amblyomma splendidum from Ethiopia; and Candidatus Rickettsia barbariae, which shared 99.80% (499/500) sequence homology with the strain detected in steppe polecats from Xinjiang Uygur Autonomous Region in China, with R. sibirica (33.80% detection rate, 97/287) and R. raoultii (33.10%, 95/287) as dominant species. A total of 1 489 serum samples were collected from individuals at a high risk of tick exposure, and the seroprevalence of antiR. sibirica IgG antibodies was 30.83% (459/1 489). There were no significant differences in the seroprevalence in terms of genders (χ² = 2.23, P > 0.05) or age groups (χ² = 6.45, P > 0.05); however, there was a region-specific seroprevalence rate of anti-R. sibirica IgG antibodies (χ² = 16.05, P < 0.05). Pairwise comparisons showed lower seroprevalence rate of anti-R. sibirica IgG antibodies among high-risk individuals in Midong District of Urumqi City (22.27%, 53/238) than in Wujiaqu City (34.70%, 93/268) and Manas County(37.55%, 98/261) (both P values < 0.05); however, they were not significantly different from those in the other regions (all P values > 0.05). Conclusion Multiple SFGR species are detected in H. asiaticum ticks sampled from the southern margin of the Gurbantunggut Desert, Xinjiang Uygur Autonomous Region, where there is a potential risk of local transmission of spotted fever rickettsiosis.

ObjectiveTo investigate the mechanism underlying the weight loss effect of electrothermal acupuncture （ETA） in diet-induced obese （DIO） model mice based on the hypoxia-inducible factor-1α （HIF-1α）/sphingosine-1-phosphate （S1P）sphingosine-1-phosphate receptor 3 （S1PR3） signaling pathway-mediated macrophage polarization. MethodsAn obesity mouse model was established by continuously feeding at a high-fat diet （HFD） for 8 weeks. Twenty successfully modeled obese mice were randomly divided into a model group and an electrothermal acupuncture （ETA） group with 10 mice per group， both of which continued to receive the HFD till week 12. Another 10 mice were assigned to the normal group and fed a standard chow diet continuously for 12 weeks. Starting from week 9 of the experiment， mice in the ETA group received electrothermal acupuncture stimulation at Zhongwan （CV12）， bilateral Tianshu （ST25）， and bilateral Zusanli （ST36） for 30 minutes per session， with the instrument current set at （40±0.5）mA and a constant temperature of （40±0.5）℃. Mice in the normal group and the model group were restrained in the same manner for 30 minutes without acupuncture. All interventions were performed once every other day for a total of 4 weeks. Ten hours after the final intervention， the body weight， waist circumference， and body length of mice in each group were measured， and Lee's index was calculated. Serum biochemical analysis was performed to assess fasting blood glucose （FBG）， total cholesterol （TC）， triglycerides （TG）， high-density lipoprotein cholesterol （HDL-C）， low-density lipoprotein cholesterol （LDL-C）， alanine aminotransferase （ALT）， aspartate aminotransferase （AST）， blood urea nitrogen （BUN）， and creatinine （Cr）. Serum leptin （LEP） and adiponectin （ADPN） levels， as well as sphingosine-1-phosphate （S1P） levels in epididymal white adipose tissue （eWAT）， were detected by ELISA. Histomorphological changes in eWAT were monitored via HE staining. Flow cytometry was used to analyze macrophage subtype-related indexes in eWAT， including the total proportion of macrophages， the proportion of M1 macrophages， the proportion of M2 macrophages， and the M1/M2 ratio. Antibody array technology was used to determine the levels of differentially expressed inflammatory factors in eWAT between the model and blank groups， as well as between the ETA and model groups. Western blot and qPCR were used to detect the protein and mRNA expression levels of HIF-1α and S1PR3 in the eWAT， respectively. ResultsCompared to the normal group， the model group showed significantly increased body weight， Lee's index， waist circumference， and FBG， TC， TG， HDL-C， LDL-C， AST， and ALT levels， upgraded LEP and S1P levels， decreased ADPN level and ADPN/LEP ratio， elevated proportion of M1 macrophages and the M1/M2 ratio， reduced proportion of M2 macrophages， and upregulated protein and mRNA expression levels of HIF-1α and S1PR3 （P＜0.05 or P＜0.01）. With the exception of HDL-C， all the abovementioned indicators in the ETA group were significantly improved compared to the model group （P＜0.05 or P＜0.01）. Analysis of differentially expressed inflammatory cytokines revealed that， compared to the normal group， the levels of interleukin-1β （IL-1β）， interferon-γ （IFN-γ）， interleukin-6 （IL-6）， and interleukin-12p70 （IL-12p70） in the eWAT of the model group were elevated； conversely， the levels of interleukin-17 （IL-17）， interleukin-5 （IL-5）， IL-6， and interleukin-23 （IL-23） in the eWAT of the ETA group decreased compared to the model group （P＜0.05 or P＜0.01）. ConclusionETA exerts a weight-reducing effect in obese model mice， and the mechanism may involve inhibition of the HIF-1α/S1P/S1PR3 signaling pathway activation， which shifts macrophage polarization from the M1 to M2 phenotype， thereby alleviating chronic inflammation.

ObjectiveTo investigate the regulatory effect of miR-21-5p in bone marrow mesenchymal stem cell-derived exosome （BMSC-Exo） on angiogenesis through the phosphatidylinositol 3-kinase （PI3K）/protein kinase B （Akt） signaling pathway. MethodsBone marrow mesenchymal stem cells （BMSCs） were isolated from rats， cultured， and divided into negative control group （NC group） and inhibitor group， which were transfected with miR-21-5p negative control and miR-21-5p inhibitor， respectively， and then exosomes were extracted and identified. Exosomes were co-cultured with human umbilical vein endothelial cells （HUVECs）. MTT assay was used to measure cell proliferation （24， 48， and 72 hours）， tube formation assay was used to assess angiogenesis ability， and Western blot was used to measure the protein expression levels of phosphorylated Akt （p-Akt）， total Akt， and vascular endothelial growth factor （VEGF）. The change in the angiogenesis ability of HUVECs was observed after treatment with the PI3K inhibitor LY294002. The independent samples t-test was used for comparison of continuous data between two groups， and a one-way analysis of variance was used for comparison between multiple groups， while the least significant difference t-test was used for further comparison between two groups. ResultsFlow cytometry showed that the third-generation BMSCs showed high expression levels of CD29 （98.18%） and CD90 （97.85%） and a low expression level of CD45 （0.22%）. Quantitative real-time reverse transcription-polymerase chain reaction showed that compared with the NC group， there was a significant reduction in the relative expression level of miR-21-5p in BMSCs after transfection with miR-21-5p inhibitor （0.38±0.05 vs 1.00±0.08， t=11.38， P<0.01）. MTT assay showed that compared with the control group， the 25， 50， 100， and 200 μg/mL exosome treatment groups had a significant increase in the proliferative activity of HUVECs （all P<0.05）， among which 100 μg/mL exosome treatment for 72 hours showed the most significant pro-proliferative effect， while the 400 μg/mL exosome treatment group showed a significant reduction in the proliferative activity of HUVECs （P<0.01）. Tube formation assay showed that compared with the control group and the NC-Exo group， the inhibitor-Exo group showed a significant reduction in angiogenesis ability （both P<0.05）， and compared with the inhibitor-Exo group， there was a further significant reduction in angiogenesis ability after addition of the PI3K inhibitor LY294002 （P<0.01）. Western blot showed that compared with the control group and the NC-Exo group， the inhibitor-Exo group had significant reduction in the protein expression levels of p-Akt and VEGF （all P<0.05）， and compared with the inhibitor-Exo group， the protein expression levels of p-Akt and VEGF further decreased significantly after the addition of LY294002 （both P<0.01）. ConclusionBMSC-Exo targets and activates the PI3K/Akt signaling pathway through miR-21-5p， upregulates the expression of VEGF， and promotes vascular endothelial cell proliferation and angiogenesis.

ObjectiveTo investigate the modifiable areal unit problem （MAUP） in the spatial pattern of Pseudostellaria heterophylla production regions and reveal the impact of statistical scales on the spatial distribution characteristics of this medicinal plant species. MethodsUsing multi-source data （literature records， field surveys， and statistical data）， we systematically analyzed the spatial patterns across three administrative levels （provincial， prefectural， and county scales）. Spatial autocorrelation （Moran's I） analysis， high-low clustering （Getis-Ord General G）， and hot/cold spot analysis （Getis-Ord Gi*） were employed. ResultsThe literature-based analysis showed that the production regions of P. heterophylla presented random distribution on the provincial scale and significant aggregation on the prefectural scale. The field survey data showed that the production regions displayed random distribution on the provincial scale but significant aggregation on both prefectural and county scales. The statistical data revealed that the production regions lacked spatial autocorrelation on the provincial scale but demonstrated significant aggregation on prefectural and county scales. ConclusionMAUP effects have substantive implications for understanding and decision-making in the arrangement of medicinal plant production regions. The county scale proves to be the most sensitive and explanatory level for analyzing the spatial pattern of P. heterophylla production regions， providing a critical foundation for habitat modeling， suitability evaluation， and ecological cultivation planning of medicinal plants.