ObjectiveTo observe the effect of Liuwei Dihuangwan on behavioral impairments in the rat model of autism spectrum disorder （ASD） and explore the mechanism of action. MethodsTwelve SD pregnant rats were intraperitoneally injected with valproic acid （VPA） （10 rats） or normal saline （2 rats）， and male offspring were selected to establish the model of ASD and the control rats. Rats were randomly assigned into model， low-dose （0.75 g·kg^-1） and high-dose （1.5 g·kg^-1） Liuwei Dihuangwan， vitamin D （positive drug， 3.7×10^-5 g·kg^-1）， and blank groups. Each group was administrated with the corresponding concentration of drugs or the same volume of normal saline by gavage for 2 weeks. After the intervention， the three-chamber social test was conducted to evaluate social interaction and social preference. The open field test was carried out to observe spontaneous behavior and anxiety state. Hematoxylin-eosin staining （HE） was used to observe the pathological changes of the prefrontal tissue. Transmission electron microscopy was employed to observe the ultrastructure of mitochondria in prefrontal neurons. Immunofluorescence was used to detect the expression of ionized calcium-binding adapter molecule-1 （Iba-1） in the prefrontal tissue. Enzyme-linked immunosorbent assay （ELISA） was adopted to measure the levels of tumor necrosis factor-α （TNF-α） and interleukin-6 （IL-6）. Western blot was employed to assess the expression differences of phosphorylated adenosine monophosphate-activated protein kinase （p-AMPK）， adenosine monophosphate-activated protein kinase （AMPK）， phosphorylated Unc-51-like autophagy-activating kinase 1 （p-ULK1）， Unc-51-like autophagy-activating kinase 1 （ULK1）， and FUN14 domain-containing protein 1 （FUNDC1）. ResultsCompared with the blank group， the model group spent less time sniffing stranger 1 and stranger 2 in the three-chamber social test （P<0.01） and showed reductions in the total distance traveled， average speed， distance traveled in the central area， and time spent in the central area in the open field test （P<0.01）. In addition， the model group showed extensive apoptosis of neurons， with shrunken nuclei and red-stained cytoplasm， and extensive necrosis of neurons in the prefrontal tissue， mitochondrial swelling， decreased matrix density， disrupted cristae， and autophagic lysosomes in neurons， increases in the rate of Iba-1 positive cells in the prefrontal area （P<0.01） and the levels of TNF-α and IL-6 （P<0.01）， and down-regulation in the expression of p-AMPK/AMPK， p-ULK1/ULK1， and FUNDC1 （P<0.01）. Compared with the model group， low-dose and high-dose Liuwei Dihuangwan and the vitamin D prolonged the time spent sniffing stranger 1 and stranger 2 in the three-chamber social test （P<0.05， P<0.01）， increased the total distance traveled， average speed， distance traveled in the central area， and time spent in the central area in the open field test （P<0.05， P<0.01）， restored the morphology of neurons in the prefrontal tissue， decreased the number of apoptotic cells， alleviated the swelling of mitochondria in neurons， increased the matrix density， mitigated the fragmentation and disorder of cristae， and increased the number of autophagosomes. Moreover， the drugs decreased the rate of Iba-1 positive cells in the prefrontal area （P<0.01）， lowered the levels of TNF-α and IL-6 （P<0.01）， and up-regulated the expression of p-AMPK/AMPK， p-ULK1/ULK1， and FUNDC1 （P<0.01）. ConclusionLiuwei Dihuangwan ameliorate autism-like behaviors and reduce neuronal apoptosis and neuroinflammatory damage in the rat model of ASD by promoting mitophagy mediated by the AMPK/ULK1/FUNDC1 pathway.

Objective:To construct a "local-systemic" dual-track prediction model integrating the resistance index (RI) score of bilateral sacroiliac joints and the systemic immune-inflammation index (SII), and to evaluate its predictive efficacy for the active stage of ankylosing spondylitis (AS).Methods:A total of 205 patients with ankylosing spondylitis (AS) from the Second Hospital of Lanzhou University between April 2022 and April 2025 were retrospectively enrolled and categorized into an active group ( n=113) and a remission group ( n=92). Hematological parameters and ultrasound data were collected. The resistance index (RI) of the synovial area in bilateral sacroiliac joints was measured by Doppler ultrasound and scored as follows: RI < 0.5: 3 points; RI 0.5~0.55: 2 points; RI > 0.55: 1 point; undetectable blood flow: 0 points. A total bilateral RI score (range 0 to 6) was calculated. The systemic immune-inflammation index (SII) was derived as (neutrophils× platelets)/lymphocytes. Normality was tested for all continuous variables; normally distributed data were compared using the t-test, while non-normally distributed data were analyzed with the Mann-Whitney U test. Categorical variables were compared using the χ2 test or analysis of variance.Variable selection was performed using Lasso regression, and a multivariate logistic regression model was developed to assess predictive performance. Results:The proportion of patients with a bilateral RI total score≥5 was significantly higher in the active group compared to the remission group (50 of 113, 44.3% vs 2 of 92, 2.2%, χ2=55.63, P<0.001). Multivariate logistic regression analysis, after adjustment for confounding variables, identified the SII [ OR(95% CI)=1.01(1.00, 1.01), P<0.001], bilateral RI total score [ OR(95% CI)=1.67(1.29, 2.26), P<0.001], erythrocyte sedimentation rate [ OR(95% CI)=1.19(1.11, 1.30), P<0.001], and mean corpuscular hemoglobin concentration [ OR(95% CI)=1.09(1.03, 1.17), P<0.001] as independent risk factors for active AS. Conversely, lymphocyte count [ OR(95% CI)=0.42(0.18, 0.92), P=0.030] and globulin [ OR(95% CI)=0.89(0.80, 0.99), P=0.040] were significantly associated with protective effects. The bilateral RI total score demonstrated the strongest predictive effect, with each 1-point increase associated with a 67% elevation in the risk of active disease. ROC curve analysis indicated that the area under the curve (AUC) for predicting whether AS is in the active disease phase was 0.94 for the combined model (SII+bilateral RI total score), compared with 0.93 for the SII-alone model and 0.92 for the bilateral RI total score-alone model, demonstrating superior predictive performance of the combined model (SII+bilateral RI total score). An online prediction tool has been developed based on the combined model. Conclusion:The dual-track prediction model, which integrates local joint hemodynamic characteristics and systemic immune-inflammatory status, facilitates a multidimensional assessment of the risk of active AS and provides an objective basis for early identification.

ObjectiveTo investigate the neuroprotective mechanism of baicalein （BAI） on Parkinson's disease （PD） model rats by regulating endoplasmic reticulum stress pathway. MethodsSeventy-two Sprague-Dawley （SD） rats were randomly divided into normal group， model group， BAI low-dose group （80 mg·kg^-1）， medium-dose group （120 mg·kg^-1）， high-dose group （160 mg·kg^-1）， and levodopa-benserazide group （51 mg·kg^-1）， with 12 rats per group. Except for the normal group， PD rat models were established by subcutaneous injection of rotenone solution （2 mg·kg^-1） into the neck back of rats in the rest of groups for consecutive 28 days. Concurrently， rats in all groups received corresponding drugs via gavage for 28 days. After treatment， behavioral changes were assessed by using the open field and pole climbing tests. Neuronal pathology and apoptosis in the substantia nigra were observed via hematoxylin-eosin （HE） staining and TdT-mediated dUTP nick-end labeling （TUNEL） assay. α-Synuclein and tyrosine hydroxylase （TH） expressions were detected by immunohistochemistry （IHC）. Inflammatory factors such as interleukin-6 （IL-6）， interleukin-1β （IL-1β）， and tumor necrosis factor-α （TNF-α） were measured by enzyme-linked immunosorbent assay （ELISA）. RNA-like endoplasmic reticulum kinase （PERK）， activating transcription factor 4 （ATF4）， C/EBP-homologous protein （CHOP）， and Bcl-2-associated X protein （Bax） expressions were analyzed by Western blot. ResultsCompared with the normal group， the model group exhibited significantly reduced locomotion distance （P<0.01） and elevated pole-climbing scores （P<0.01）， with increased neuronal apoptosis rate （P<0.01）， significantly enhanced α-Synuclein expression （P<0.01）， decreased TH expression （P<0.01）， upregulated release of inflammatory factors （P<0.05，P<0.01）， and increased protein expressions of PERK/ATF4 pathway proteins and pro-apoptotic Bax （P<0.05，P<0.01）. Compared with the model group， medium/high-dose BAI groups and levodopa-benserazide group showed obviously improved motor function （P<0.05，P<0.01）， reduced pole-climbing scores （P<0.05）， decreased neuronal apoptosis （P<0.01）， downregulated α-Synuclein expression （P<0.01）， upregulated TH expression （P<0.05，P<0.01）， suppressed release of inflammatory factors （P<0.05，P<0.01）， and decreased protein expressions of PERK/ATF4 pathway proteins and pro-apoptotic Bax （P<0.05，P<0.01）. ConclusionBAI reduces the release of neuroinflammatory factors and neuronal apoptosis to improve the neurological function of PD model rats， and its mechanism may be related to alleviating endoplasmic reticulum stress and apoptosis by regulating the PERK/ATF4 pathway.

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.

ObjectiveTo investigate the neuroprotective mechanism of baicalein （BAI） on Parkinson's disease （PD） model rats by regulating endoplasmic reticulum stress pathway. MethodsSeventy-two Sprague-Dawley （SD） rats were randomly divided into normal group， model group， BAI low-dose group （80 mg·kg^-1）， medium-dose group （120 mg·kg^-1）， high-dose group （160 mg·kg^-1）， and levodopa-benserazide group （51 mg·kg^-1）， with 12 rats per group. Except for the normal group， PD rat models were established by subcutaneous injection of rotenone solution （2 mg·kg^-1） into the neck back of rats in the rest of groups for consecutive 28 days. Concurrently， rats in all groups received corresponding drugs via gavage for 28 days. After treatment， behavioral changes were assessed by using the open field and pole climbing tests. Neuronal pathology and apoptosis in the substantia nigra were observed via hematoxylin-eosin （HE） staining and TdT-mediated dUTP nick-end labeling （TUNEL） assay. α-Synuclein and tyrosine hydroxylase （TH） expressions were detected by immunohistochemistry （IHC）. Inflammatory factors such as interleukin-6 （IL-6）， interleukin-1β （IL-1β）， and tumor necrosis factor-α （TNF-α） were measured by enzyme-linked immunosorbent assay （ELISA）. RNA-like endoplasmic reticulum kinase （PERK）， activating transcription factor 4 （ATF4）， C/EBP-homologous protein （CHOP）， and Bcl-2-associated X protein （Bax） expressions were analyzed by Western blot. ResultsCompared with the normal group， the model group exhibited significantly reduced locomotion distance （P<0.01） and elevated pole-climbing scores （P<0.01）， with increased neuronal apoptosis rate （P<0.01）， significantly enhanced α-Synuclein expression （P<0.01）， decreased TH expression （P<0.01）， upregulated release of inflammatory factors （P<0.05，P<0.01）， and increased protein expressions of PERK/ATF4 pathway proteins and pro-apoptotic Bax （P<0.05，P<0.01）. Compared with the model group， medium/high-dose BAI groups and levodopa-benserazide group showed obviously improved motor function （P<0.05，P<0.01）， reduced pole-climbing scores （P<0.05）， decreased neuronal apoptosis （P<0.01）， downregulated α-Synuclein expression （P<0.01）， upregulated TH expression （P<0.05，P<0.01）， suppressed release of inflammatory factors （P<0.05，P<0.01）， and decreased protein expressions of PERK/ATF4 pathway proteins and pro-apoptotic Bax （P<0.05，P<0.01）. ConclusionBAI reduces the release of neuroinflammatory factors and neuronal apoptosis to improve the neurological function of PD model rats， and its mechanism may be related to alleviating endoplasmic reticulum stress and apoptosis by regulating the PERK/ATF4 pathway.

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.

OBJECTIVE To observe the effects of endogenous metabolite mesaconate combined with Sodium lactate Ringer’s injection （LR） on prolonging the golden treatment time window and its resuscitation efficacy in rats with hemorrhagic shock under high-altitude conditions. METHODS Rats were divided into the shock group， LR group， and 5， 20， 50 mg/kg mesaconate+LR groups， with 20 rats in each group， to investigate the effect of additional use of mesaconate on the golden treatment time window. After establishing a model of uncontrolled hemorrhagic shock under high-altitude conditions in all groups by housing in a hypobaric hypoxia chamber combined with splenic artery transection， rats in the shock group received no resuscitation， while rats in the LR group and mesaconate+LR groups underwent low-pressure resuscitation with LR or mesaconate combined with LR. Blood pressure control， fluid infusion volume， blood loss rate and survival status were observed in each group. Rats were further divided into the normal group， shock group and mesaconate （50 mg/kg）+LR group， with 10 or 20 rats in each group， to evaluate the resuscitation effects after extending the golden treatment time window by additionally using mesaconate. Except for the normal group， the other groups underwent the same procedure to establish an uncontrolled hemorrhagic shock model under high-altitude conditions. Rats in the shock group received no resuscitation. In the mesaconate+LR group， after 3 h of low-pressure resuscitation， bleeding control was performed by ligation of the spleen artery， and the infusion volume and blood loss rate were recorded； subsequently， the rats received LR resuscitation with twice the volume of blood loss. Then， blood gas indicators of the mesaconate+LR group were measured at different time points. Survival rates， indicators related to sublingual microcirculatory perfusion， liver and kidney blood flow， indicators related to the function of vital organs， and lung and brain water content were observed in all groups. RESULTS LR infusion alone could effectively maintain mean arterial pressure （MAP） within 50-60 mmHg for approximately 1 h. The administration of mesaconate combined with LR during hypotensive resuscitation could maintain MAP within 50-60 mmHg for over 3 h， with significantly reduced fluid infusion volume and blood loss rate in 50 mg/kg mesaconate+LR group， compared to the LR group （ P ＜0.05）. In the LR group， rats maintained low pressure for up to 1 hour with a survival rate of 52.94%， and no rats survived beyond 2 h. In the 5， 20 and 50 mg/kg mesaconate+LR groups， rats maintained low pressure for up to 1 h with a survival rate exceeding 80%； in the 20 and 50 mg/kg mesaconate+LR groups， rats maintained low pressure for up to 3 h with a survival rate exceeding 70%. After complete resuscitation with mesaconate combined with LR， the 72 h survival rate of rats was 43.75%， and significant improvements in blood gas parameters were observed compared to the end of the shock phase （ P ＜0.05）. Compared to the shock group， the mesaconate+LR group showed significant recovery in sublingual microcirculatory indicators， and liver/kidney blood flow after complete resuscitation （ P ＜0.05）， with significant reductions in heart， liver and kidney function-related indicators and lung water content （ P ＜0.05）. CONCLUSIONS Mesaconate combined with LR significantly extends the golden treatment time window for hemorrhagic shock in rats under high-altitude conditions， improves blood gas parameters， sublingual microcirculatory perfusion， and liver/kidney blood flow， mitigates vital organ impairment and pulmonary edema， and increases the survival rate of shocked rats.

Objective:To construct a "local-systemic" dual-track prediction model integrating the resistance index (RI) score of bilateral sacroiliac joints and the systemic immune-inflammation index (SII), and to evaluate its predictive efficacy for the active stage of ankylosing spondylitis (AS).Methods:A total of 205 patients with ankylosing spondylitis (AS) from the Second Hospital of Lanzhou University between April 2022 and April 2025 were retrospectively enrolled and categorized into an active group ( n=113) and a remission group ( n=92). Hematological parameters and ultrasound data were collected. The resistance index (RI) of the synovial area in bilateral sacroiliac joints was measured by Doppler ultrasound and scored as follows: RI < 0.5: 3 points; RI 0.5~0.55: 2 points; RI > 0.55: 1 point; undetectable blood flow: 0 points. A total bilateral RI score (range 0 to 6) was calculated. The systemic immune-inflammation index (SII) was derived as (neutrophils× platelets)/lymphocytes. Normality was tested for all continuous variables; normally distributed data were compared using the t-test, while non-normally distributed data were analyzed with the Mann-Whitney U test. Categorical variables were compared using the χ2 test or analysis of variance.Variable selection was performed using Lasso regression, and a multivariate logistic regression model was developed to assess predictive performance. Results:The proportion of patients with a bilateral RI total score≥5 was significantly higher in the active group compared to the remission group (50 of 113, 44.3% vs 2 of 92, 2.2%, χ2=55.63, P<0.001). Multivariate logistic regression analysis, after adjustment for confounding variables, identified the SII [ OR(95% CI)=1.01(1.00, 1.01), P<0.001], bilateral RI total score [ OR(95% CI)=1.67(1.29, 2.26), P<0.001], erythrocyte sedimentation rate [ OR(95% CI)=1.19(1.11, 1.30), P<0.001], and mean corpuscular hemoglobin concentration [ OR(95% CI)=1.09(1.03, 1.17), P<0.001] as independent risk factors for active AS. Conversely, lymphocyte count [ OR(95% CI)=0.42(0.18, 0.92), P=0.030] and globulin [ OR(95% CI)=0.89(0.80, 0.99), P=0.040] were significantly associated with protective effects. The bilateral RI total score demonstrated the strongest predictive effect, with each 1-point increase associated with a 67% elevation in the risk of active disease. ROC curve analysis indicated that the area under the curve (AUC) for predicting whether AS is in the active disease phase was 0.94 for the combined model (SII+bilateral RI total score), compared with 0.93 for the SII-alone model and 0.92 for the bilateral RI total score-alone model, demonstrating superior predictive performance of the combined model (SII+bilateral RI total score). An online prediction tool has been developed based on the combined model. Conclusion:The dual-track prediction model, which integrates local joint hemodynamic characteristics and systemic immune-inflammatory status, facilitates a multidimensional assessment of the risk of active AS and provides an objective basis for early identification.

Yin-Yang theory is one of the core components of the foundational theory of traditional Chinese medicine (TCM), yet its classification as philosophical or scientific has long been controversial. To move beyond the longstanding either-philosophy-or-science debate surrounding Yin-Yang theory in TCM, this article, for the first time, proposes a three-level framework spanning macro, meso, and micro perspectives. At the macro level, Yin and Yang are a theoretical achievement of ancient Chinese natural philosophy, used to explain universal laws governing how the cosmos operates, which constitutes Yin-Yang theory in philosophy; at the meso level, Huangdi Neijing (《黄帝内经》, Inner Canon of Huangdi) translated Yin-Yang from a cosmological framework into medical discourse, making it a core methodological basis for clinical pattern differentiation and treatment, which constitutes Yin-Yang theory in medicine; at the micro level, Yin-Yang theory in TCM has increasingly converged with modern science, showing scientific features that are testable and reproducible, which constitutes Yin-Yang theory in science. The article further demonstrates that Yin-Yang theory in TCM differs markedly from Yin-Yang as a general philosophical doctrine in core application and practical orientation: in the course of its medical transformation, elements related to life phenomena were selectively incorporated, while, within a medical context, abstract propositions such as cosmological speculation were bracketed and rendered concrete, thereby achieving a practice-oriented transition from philosophy to medicine. Based on this, the present study conducts analyses from four aspects, namely philosophical roots, clinical application, modern scientific interpretation, and implications for life sciences; it substantiates the logical basis for Yin-Yang theory in TCM as the overarching framework for clinical pattern differentiation and treatment, reviews modern research progress from perspectives such as systems science and network regulation, and explores its potential value in advancing a new paradigm of state-based medicine. The study suggests that an accurate understanding of Yin-Yang theory in TCM requires moving beyond the either-philosophy-or-science binary and analyzing specific issues on a case-by-case basis. Across different perspectives, Yin-Yang theory in TCM has philosophical, medical, and scientific attributes. Only by grounding the theory in actual clinical efficacy and engaging with modern science can the innovative development of Yin-Yang theory in TCM be promoted, thereby providing valuable Eastern insights for building a life science system with original Chinese contributions.

Background: Advances in intensive care have markedly improved survival after severe brain injury, leading to a growing population of patients with prolonged disorders of consciousness (pDOC). Current management of pDOC remains largely supportive, and evidence-based neuromodulatory interventions are limited; moreover, existing guidelines provide insufficiently explicit recommendations regarding mechanisms of action and objective biomarkers of treatment response. Transcutaneous auricular vagus nerve stimulation (taVNS) has emerged as a potential noninvasive intervention; however, its modulatory effects on brain function in pDOC are not yet well characterized, and the paucity of integrative mechanistic evidence has constrained its translation into routine clinical practice. Objectives: Within a multimodal assessment framework, this study aims to systematically elucidate the neurobiological mechanisms by which taVNS modulates brain function and autonomic activity in patients with pDOC, and to evaluate its clinical potential to enhance levels of consciousness. Methods: In this randomized controlled trial, 60 patients with vegetative state/minimally conscious state will be enrolled and randomly allocated to a taVNS group, a transcutaneous nonauricular vagus nerve stimulation group (sham), or a control group (n = 20 per group) for a 4-week intervention. The primary outcome will be changes in the Coma Recovery Scale-Revised scores from baseline to weeks 1, 2, and 4 of treatment. Secondary outcomes will include functional brain activity assessed by electroencephalography and functional near-infrared spectroscopy, as well as autonomic modulation indexed by heart rate variability. Functional prognosis will be evaluated using the Glasgow Outcome Scale-Extended at the end of treatment and at a 6-month follow-up. Safety will be assessed by continuous monitoring and documentation of adverse events throughout the study period. Results and discussion: By integrating electroencephalography–functional near-infrared spectroscopy with heart rate variability, this study will characterize the effects of taVNS on functional brain networks and consciousness recovery in pDOC across complementary behavioral, electrophysiological, hemodynamic, and autonomic domains, while interrogating potential sources of clinical and neurobiological heterogeneity. The findings are expected to provide a mechanistic and evidence-based foundation for the mechanism-driven clinical implementation of taVNS and the optimization of stimulation protocols in pDOC. Clinical trial registration: International Traditional Medicine Clinical Trial Registry, ITMCTR20250021041, https://itmctr.ccebtcm.org.cn.