The initiation of adaptive immune responses relies on the precise recognition and interpretation of antigenic information. In this process, the specific binding of T cell receptors (TCRs) to peptide-major histocompatibility complex (pMHC) molecules represents one of the key molecular events in the initiation of adaptive immune responses. Accordingly, the structural features of TCR-pMHC complexes provide a fundamental basis for dissecting antigen recognition mechanisms and support rational vaccine design, therapeutic target discovery in TCR-based immunotherapy, and TCR identification and optimization. However, experimental determination of TCR-pMHC structures remains costly, time-consuming, and limited in coverage, making computational approaches essential for rapidly obtaining reliable structural information. Computational methods for predicting the structures of TCR-pMHC complexes have advanced rapidly in recent years, driven by progress in deep learning-based modeling frameworks and the increasing availability of structural and sequence resources. Despite these developments, most existing tools do not adequately distinguish the key structural and biophysical differences between MHC class I (MHC-I) and MHC class II (MHC-II) complexes during model construction. As a consequence, their predictive performance differs substantially between class I and class II complexes. In general, structural predictions for class I complexes outperform those for class II complexes. This discrepancy may be related to several fundamental differences between the two systems, including the architecture of the peptide-binding groove, the distribution of peptide lengths, and the properties of peptide flanking residues (PFRs). Compared with MHC-I molecules, MHC-II molecules usually bind longer antigenic peptides, which typically range from 13 to 25 amino acids in length. PFRs at both termini of these peptides participate in regulating the overall conformation of TCR-pMHC class II complexes and exert a pronounced effect on the geometric and physicochemical characteristics of the TCR-pMHC binding interface. Furthermore, within the TCR recognition interface, the complementarity-determining regions (CDRs) consist of segments that differ markedly in conformational behavior. They commonly include regions that are relatively rigid and structurally stable, together with highly flexible segments exhibiting substantial conformational plasticity. These rigidity-flexibility features constitute an essential structural basis enabling TCRs to recognize diverse peptide-MHC ligands and to accommodate conformational heterogeneity at the interface. However, many current modeling tools, in an effort to enforce global conformational stability or reduce structural noise, tend to over-constrain intrinsically flexible regions. Such oversimplification may lead to inappropriate rigidification of flexible CDR loops, resulting in local structural distortions, compromised interface geometry, or even complete modeling failure for specific complexes. Against this background, the review approaches the field from the perspective of computational differences between MHC-I and MHC-II complexes. We first systematically organize and summarize available resources related to TCRs and pMHCs, including structural datasets, sequence databases, prediction tools, and benchmarking studies. We then focus on five representative tools capable of predicting both class I and class II complexes—AlphaFold2, AlphaFold3, TCRmodel2, tFold-TCR, and TCR-pHLA_ModellerS. After excluding structures present in the training sets of these tools, we constructed a benchmark dataset comprising 25 class I and 10 class II TCR-pMHC complexes in the bound state and conducted a systematic evaluation using this dataset. We first employ widely used general evaluation metrics, including All-Atom Root Mean Square Deviation (All-Atom RMSD), Backbone RMSD, Template Modeling score (TM-score), and DockQ, to assess the global conformational accuracy and interface modeling quality of class I and class II complexes. For class II complexes, we propose for the first time a peptide flanking residue deviation index, including the PFRs-Deviation Index (PFRs-DI), N-PFR-Deviation Index (N-PFR-DI), and C-PFR-Deviation Index (C-PFR-DI), to quantitatively characterize conformational deviations in PFRs. In addition, we propose the CDR conformational consistency index (CCC) designed to qualitatively evaluate the ability of prediction tools to capture TCR CDR conformational flexibility. These metrics collectively assess a tool’s ability to model both overall conformation and critical functional regions, thereby addressing the limitations of existing evaluation criteria that overemphasize global structure while inadequately capturing modeling quality in key functional areas. This establishes a unified analytical framework for MHC-I and MHC-II complexes to guide data resource selection, modeling strategy formulation, and evaluation system development. The framework further advances computational modeling and provides crucial support for multi-scale analysis of TCR-pMHC recognition mechanisms and their biological functions.

The initiation of adaptive immune responses relies on the precise recognition and interpretation of antigenic information. In this process, the specific binding of T cell receptors (TCRs) to peptide-major histocompatibility complex (pMHC) molecules represents one of the key molecular events in the initiation of adaptive immune responses. Accordingly, the structural features of TCR-pMHC complexes provide a fundamental basis for dissecting antigen recognition mechanisms and support rational vaccine design, therapeutic target discovery in TCR-based immunotherapy, and TCR identification and optimization. However, experimental determination of TCR-pMHC structures remains costly, time-consuming, and limited in coverage, making computational approaches essential for rapidly obtaining reliable structural information. Computational methods for predicting the structures of TCR-pMHC complexes have advanced rapidly in recent years, driven by progress in deep learning-based modeling frameworks and the increasing availability of structural and sequence resources. Despite these developments, most existing tools do not adequately distinguish the key structural and biophysical differences between MHC class I (MHC-I) and MHC class II (MHC-II) complexes during model construction. As a consequence, their predictive performance differs substantially between class I and class II complexes. In general, structural predictions for class I complexes outperform those for class II complexes. This discrepancy may be related to several fundamental differences between the two systems, including the architecture of the peptide-binding groove, the distribution of peptide lengths, and the properties of peptide flanking residues (PFRs). Compared with MHC-I molecules, MHC-II molecules usually bind longer antigenic peptides, which typically range from 13 to 25 amino acids in length. PFRs at both termini of these peptides participate in regulating the overall conformation of TCR-pMHC class II complexes and exert a pronounced effect on the geometric and physicochemical characteristics of the TCR-pMHC binding interface. Furthermore, within the TCR recognition interface, the complementarity-determining regions (CDRs) consist of segments that differ markedly in conformational behavior. They commonly include regions that are relatively rigid and structurally stable, together with highly flexible segments exhibiting substantial conformational plasticity. These rigidity-flexibility features constitute an essential structural basis enabling TCRs to recognize diverse peptide-MHC ligands and to accommodate conformational heterogeneity at the interface. However, many current modeling tools, in an effort to enforce global conformational stability or reduce structural noise, tend to over-constrain intrinsically flexible regions. Such oversimplification may lead to inappropriate rigidification of flexible CDR loops, resulting in local structural distortions, compromised interface geometry, or even complete modeling failure for specific complexes. Against this background, the review approaches the field from the perspective of computational differences between MHC-I and MHC-II complexes. We first systematically organize and summarize available resources related to TCRs and pMHCs, including structural datasets, sequence databases, prediction tools, and benchmarking studies. We then focus on five representative tools capable of predicting both class I and class II complexes—AlphaFold2, AlphaFold3, TCRmodel2, tFold-TCR, and TCR-pHLA_ModellerS. After excluding structures present in the training sets of these tools, we constructed a benchmark dataset comprising 25 class I and 10 class II TCR-pMHC complexes in the bound state and conducted a systematic evaluation using this dataset. We first employ widely used general evaluation metrics, including All-Atom Root Mean Square Deviation (All-Atom RMSD), Backbone RMSD, Template Modeling score (TM-score), and DockQ, to assess the global conformational accuracy and interface modeling quality of class I and class II complexes. For class II complexes, we propose for the first time a peptide flanking residue deviation index, including the PFRs-Deviation Index (PFRs-DI), N-PFR-Deviation Index (N-PFR-DI), and C-PFR-Deviation Index (C-PFR-DI), to quantitatively characterize conformational deviations in PFRs. In addition, we propose the CDR conformational consistency index (CCC) designed to qualitatively evaluate the ability of prediction tools to capture TCR CDR conformational flexibility. These metrics collectively assess a tool’s ability to model both overall conformation and critical functional regions, thereby addressing the limitations of existing evaluation criteria that overemphasize global structure while inadequately capturing modeling quality in key functional areas. This establishes a unified analytical framework for MHC-I and MHC-II complexes to guide data resource selection, modeling strategy formulation, and evaluation system development. The framework further advances computational modeling and provides crucial support for multi-scale analysis of TCR-pMHC recognition mechanisms and their biological functions.

In the central nervous system (CNS), the myelin sheath, a specialized membrane structure that wraps around axons, is formed by oligodendrocytes through a highly coordinated spatiotemporal developmental program. The process begins with the directed differentiation of neural precursor cells into oligodendrocyte precursor cells (OPCs), followed by their migration, proliferation, differentiation, and maturation, ultimately leading to the formation of a multi-segmental myelin sheath structure. Recent single-cell sequencing research has revealed that this process involves the temporal regulation of over 200 key genes, with a regulatory network composed of transcription factors such as Sox10 and Olig2 playing a central role. The primary function of the myelin sheath is to accelerate nerve signal transmission and protect nerve fibers from damage. Its insulating properties not only increase nerve conduction speed by 50-100 times but also ensure the long-term functional integrity of the nervous system by maintaining axonal metabolic homeostasis and providing mechanical protection. The pathological effects of myelin sheath injury exhibit a cascade amplification pattern: acute demyelination leads to action potential conduction block, while chronic lesions may cause axonal damage and neuronal death in severe or long-term cases, ultimately resulting in irreversible neurological dysfunction with neurodegenerative characteristics. Multiple sclerosis (MS) is a neurodegenerative disease characterized by chronic inflammatory demyelination of the CNS. Clinically, the distribution of lesions in MS exhibits spatial heterogeneity, which is closely related to differences in the regenerative capacity of oligodendrocytes within the local microenvironment. Emerging evidence suggests that astrocytes form a dynamic “neural-immune-metabolic interface” and play a multidimensional regulatory role in myelin development and regeneration by forming heterogeneous populations composed of different subtypes. During embryonic development, astrocytes induce the targeted differentiation of OPCs in the ventricular region through the Wnt/β-catenin pathway. In the mature stage, they secrete platelet-derived growth factor AA (PDGF-AA) to establish a chemical gradient that guides the precise migration of OPCs along axonal bundles. Notably, astrocytes also provide crucial metabolic support by supplying energy substrates for high-energy myelin formation through the lactate shuttle mechanism. In addition, astrocytes play a dual role in myelin regulation. During the acute injury phase, reactive astrocytes establish a triple defense system within 72 h: upregulating glial fibrillary acidic protein (GFAP) to form scars that isolate lesions, activating the JAK-STAT3 regeneration pathway in oligodendrocytes via leukemia inhibitory factor (LIF), and releasing tumor necrosis factor-stimulated gene-6 (TSG-6) to inhibit excessive microglial activation. However, in chronic neurodegenerative diseases, the phenotypic transformation of astrocytes contributes to microenvironmental deterioration. The secretion of chondroitin sulfate proteoglycans (CSPGs) inhibits OPC migration via the RhoA/ROCK pathway, while the persistent release of reactive oxygen species (ROS) leads to mitochondrial dysfunction and the upregulation of complement C3-mediated synaptic pruning. This article reviews the mechanisms by which astrocytes regulate the development and regeneration of myelin sheaths in the CNS, with a focus on analyzing the multifaceted roles of astrocytes in this process. It emphasizes that astrocytes serve as central hubs in maintaining myelin homeostasis by establishing a metabolic microenvironment and signaling network, aiming to provide new therapeutic strategies for neurodegenerative diseases such as multiple sclerosis.

Objective To investigate the clinical efficacy and safety of facilitated percutaneous coronary intervention(PCI)with half-dose recombinant staphylokinase(r-SAK)in patients with ST-segment elevation myocardial infarction(STEMI)who are expected to undergo PCI within 120 minutes.Methods From October 2021 to August 2022,a total of 200 STEMI patients in eight centers were included and randomly assigned in a 1﹕1 ratio to either r-SAK group or control group.Patients received loading doses of aspirin and ticagrelor and intravenous heparin and were randomized to receive an intravenous bolus of either 5 mg r-SAK or normal saline prior to PCI.The outcomes were set as ST-segment resolution(STR)at 60-90 minutes after PCI,the proportion and transition of pathological Q waves on the 5th day after PCI,and the proportion of high-sensitivity cardiac troponin T(hs-cTnT)peaking within 12 hours of onset.The safety outcome was major bleeding events defined as Bleeding Academic Research Consortium(BARC)≥type 3 bleeding during hospitalization.Results Compared with the control group,the r-SAK group had a higher proportion of STR≥70%within 60-90 minutes after PCI(58.3%vs.40.3%,P=0.009);a lower proportion of pathological Q waves(59.1%vs.74.1%,P=0.040);a lower rate of Q wave progression(14.8%vs.43.2%,P＜0.001);a higher rate of Q wave disappearance(12.5%vs.3.7%,P=0.027);and a higher proportion of hs-cTnT peaking within 12 hours of symptom onset[31/40(77.5%)vs.17/33(51.5%),P=0.027].Regarding the safety outcome,no significant difference in BARC≥type 3 bleeding was found between the two groups during hospitalization(P＞0.05).Conclusions For STEMI patients who were expected to undergo primary PCI within 120 minutes of symptom onset,the facilitated PCI with half-dose r-SAK significantly increased the proportion of STR≥70%at 60-90 minutes after PCI,reduced the formation of pathological Q waves,and shortened the time to peak hs-cTnT,without increasing the risk of bleeding,which should be an alternative reperfusion strategy worthy of further study.

Aim To evaluate the role of the glucose transporter protein 1(GLUT1)-dependent glycolytic in the attenuation of oxygen-glucose deprivation-reoxygen-ation(OGD/R)injury in HK-2 cells by dexmedetomi-dine(Dex).Methods C57/BL6 mice were random-ly divided into three groups(n=6),namely,sham operation group(Sham group),renal ischemia reper-fusion group(I/R group)and Dex group(I/R+Dex group).Serum creatinine(Cr)and urea nitrogen(BUN)were measured,while the levels of key glyco-lytic enzymes HK2,PFKFB3 and GLUT1 were meas-ured.HK-2 cells were cultured and randomised into seven groups(n=6),which was treated with OGD/R,overexpression or interference with GLUT1,Dex and glycolysis inhibitor 2-DG.CCK-8 and LDH activi-ty were used to detect cellular damage.Glycolysis lev-els were detected by lactate and ECAR.The inflamma-tory level was reflected by qRT-PCR for IL-6 and TNF-α.qRT-PCR and Western blot were performed to de-tect the levels of GLUT1,HK2,and PFKFB3.Results Dex significantly ameliorated kidney injury and HK-2 cell injury(P＜0.05).Dex inhibited the OGD/R-induced rise in lactate and extracellular acidification rate(ECAR),as evidenced by suppression of the ex-pression of GLUT1,HK2 and PFKFB3(P＜0.05).In vitro experiments showed that GLUT1 knockdown sig-nificantly improved OGD/R-induced cellular damage.Lactate,ECAR,glycolysis-related mRNAs and pro-teins were inhibited by GLUT1 knockdown(P＜0.05).Significantly,there were no significant differ-ences in above indexes after Dex treatment based on GLUT1 knockdown.Overexpression of GLUT1 abroga-ted the protective effects of Dex,while reversing the inhibitory effects of Dex on the expression of GLUT1,HK2,and PFKFB3(P＜0.05).Conclusions Dexmedetomidine attenuates OGD/R induced injury in HK-2 cells by inhibiting GLUT1-dependent glycolysis.

Objective To develop dynamic prediction models based on multiple postnatal time points to support early diagnosis and individualized intervention for bronchopulmonary dysplasia(BPD)in preterm infants with gestational age<32 weeks.Methods Clinical data of 472 preterm infants with gestational age<32 weeks admitted to the Second Xiangya Hospital of Central South University between January 2016 and November 2020 were retrospectively analyzed.Multivariable logistic regression was applied to develop five independent prediction models at postnatal days 1,7,14,21,and 28.The performance of the models was assessed using the area under the receiver operating characteristic curve(AUC)and the Hosmer-Lemeshow test.Results Baseline characteristics such as gestational age and birth weight differed significantly between the BPD group(n=147)and the non-BPD group(n=325)(P<0.05).Predictors of BPD evolved across time points:on day 1,key predictors included gestational age,birth weight,Score for Neonatal Acute Physiology II(SNAP-II),invasive mechanical ventilation,and fraction of inspired oxygen>30%;by day 7,additional variables emerged,including fasting duration>2 days,mean feeding advancement rate<8.5 mL/(kg·d),neonatal respiratory distress syndrome,apnea of prematurity,and positive sputum culture;from day 14 onward,nutrition-and treatment-related indicators were incorporated additionally.The models demonstrated good discrimination at postnatal days 1,7,14,21,and 28,with AUCs of 0.917,0.927,0.939,0.944,and 0.968,respectively,and good calibration(Hosmer-Lemeshow P>0.05).Internal validation showed AUCs ranging from 0.899 to 0.958,indicating robust performance.Conclusions Dynamic postnatal prediction models incorporating indicators spanning perinatal factors,respiratory support,nutritional management,and therapeutic interventions demonstrate high predictive performance and facilitate dynamic risk assessment for BPD in preterm infants with gestational age<32 weeks.

This study examines the effects of zearalenone(ZEA)on the survival and function of lu-teinized granulosa cells,and studies the role of activating transcription factor 6(ATF6)in regula-ting apoptosis and functional abnormalities of luteinized granulosa cells induced by ZEA.An in vitro model of luteinized granulosa cells was utilized to examine the effects of ZEA treatment on apoptosis,hormone secretion,and the expression of relevant proteins.Furthermore,the expression of ATF6 was manipulated using siRNA to elucidate its regulatory function in the ZEA-induced damage of luteinized granulosa cells in mice.Our findings revealed that ZEA inhibited the activity of luteinized granulosa cells and reduced the secretion of estradiol(E2)and progesterone(P4)in a dose-dependent manner.The expression levels of p-IRE1,ATF6 and StAR in both low(20 pmol/L)and high(40 μmol/L)ZEA groups were significantly increased after 24 h(P＜0.05).GRP78 had no significant change at low concentration treatment(P＞0.05),but significantly increased at high concentration treatment(P＜0.05).Similarly,ATF4 and p-EIF2α had no significant change at low concentration treatment(P＞0.05),but significantly decreased at high concentration treat-ment(P＜0.05).HSD3B2 and CYP19A1 were significantly decreased in both low and high concentration treatments(P＜0.05).After 48 h of treatment,ATF6 and GRP78 were significantly increased in both low and high concentration treatments(P＜0.05).p-IRE1 was significantly de-creased at low concentration treatment(P＜0.05),but remained unchanged at high concentration treatment(P＞0.05).ATF4,p-EIF2α,HSD3B2 and CYP19A1 were significantly decreased in both low and high concentration treatments(P＜0.05).St AR was significantly increased in both low and high concentration treatments(P＜0.05).Interference with the expression of ATF6 could sig-nificantly reduce the apoptosis induced by low concentration group(P＜0.05),and enhanced the hormone secretion in both high and low concentration groups(P＜0.05).In conclusion,ZEA can cause damage to luteinized granulosa cells and activate ATF6 signaling pathway.Interference with ATF6 can alleviate apoptosis and hormone secretion disturbance induced by low concentration ZEA,but has limited effect on damage caused by high concentration ZEA.

This article aims to explore the effect and mechanism of Liujunzi Pills on the intestinal microbiota of rats with spleen Qi deficiency syndrome. The raw Rhei Radix et Rhizoma water extract(1 g·mL~(-1)) was used to prepare spleen Qi deficiency rat models. A total of 44 SD male rats were randomly divided into a control group, a model group, Liujunzi Pills groups at high(3.24 g·kg~(-1)), medium(1.62 g·kg~(-1)), low(0.81 g·kg~(-1)) doses, and Shenling Baizhu San(2.50 g·kg~(-1)) group. The drug effect was evaluated by observing the following aspects: spleen index, fecal water content, body weight, and intestinal propulsion index. Gut microbiota analysis and 16S rRNA gene sequencing were conducted on feces. Enzyme-linked immunosorbent assay(ELISA) and UV spectrophotometry were used to detect interleukin-1β(IL-1β) and adenosine triphosphate(ATP) levels in small intestine tissues. Hematoxylin-eosin staining and transmission electron microscopy were employed to observe changes in intestinal pathology and microstructure. The results show that, compared with the control group, fecal moisture content is significantly increased while spleen index, body weight, and intestinal propulsion index are significantly reduced in rats of the model group, indicating the successful establishment of the model. The above symptoms can be improved by both Shenling Baizhu San and Liujunzi Pills. Compared with the control group, in the model group, the gut microbiota abundance is changed with an unbalanced development: the abundance of beneficial bacteria within the Bacteroidetes phylum is reduced, accompanied by a significantly decreased Shannon index, and reduced signal levels of nicotinamide adenine dinucleotide phosphate(NADPH)-related enzymes relevant to mitochondria. However, Liujunzi Pills and Shenling Baizhu San can significantly improve the Bacteroidetes phylum abundance in gut microbiota, microbial diversity, and NADPH activity in the model group. Additionally, compared with the control group, the ATP level is decreased and the IL-1β level is increased in small intestinal tissues of the model group, with shorter small intestinal epithelial villi and decreased mitochondrial number. The above symptoms can be improved by Liujunzi Pills and Shenling Baizhu San. In conclusion, Liujunzi Pills can treat spleen Qi deficiency syndrome by enhancing mitochondrial function to regulate gut microbiota balance and diversity.
Animals ; Gastrointestinal Microbiome/drug effects* ; Drugs, Chinese Herbal/pharmacology* ; Male ; Rats, Sprague-Dawley ; Rats ; Qi ; Spleen/metabolism* ; Splenic Diseases/metabolism* ; Humans ; Interleukin-1beta/genetics* ; Bacteria/drug effects* ; Feces/microbiology* ; Adenosine Triphosphate/metabolism*

The study aims to examine the effects and potential mechanisms of disulfiram (DSF) on cardiac hypertrophic injury, focusing on the role of transforming growth factor-β-activated kinase 1 (TAK1)-mediated pan-apoptosis (PANoptosis). H9C2 cardiomyocytes were treated with angiotensin II (Ang II, 1 µmol/L) to establish an in vitro model of myocardial hypertrophy. DSF (40 µmol/L) was used to treat cardiomyocyte hypertrophic injury models, either along or in combination with the TAK1 inhibitor, 5z-7-oxozeaenol (5z-7, 0.1 µmol/L). We assessed cell damage using propidium iodide (PI) staining, measured cell viability with CCK8 assay, quantified inflammatory factor levels in cell culture media via ELISA, detected TAK1 and RIPK1 binding rates using immunoprecipitation, and analyzed the protein expression levels of key proteins in the TAK1-mediated PANoptosis pathway using Western blot. In addition, the surface area of cardiomyocytes was measured with Phalloidin staining. The results showed that Ang II significantly reduced the cellular viability of H9C2 cardiomyocytes and the binding rate of TAK1 and RIPK1, significantly increased the surface area of H9C2 cardiomyocytes, PI staining positive rate, levels of inflammatory factors [interleukin-1β (IL-1β), IL-18, and tumor necrosis factor α (TNF-α)] in cell culture media and p-TAK1/TAK1 ratio, and significantly up-regulated key proteins in the PANoptosis pathway [pyroptosis-related proteins NLRP3, Caspase-1 (p20), and GSDMD-N (p30), apoptosis-related proteins Caspase-3 (p17), Caspase-7 (p20), and Caspase-8 (p18), as well as necroptosis-related proteins p-MLKL, RIPK1, and RIPK3]. DSF significantly reversed the above changes induced by Ang II. Both 5z-7 and exogenous IL-1β weakened these cardioprotective effects of DSF. These results suggest that DSF may alleviate cardiac hypertrophic injury by inhibiting TAK1-mediated PANoptosis.
Animals ; MAP Kinase Kinase Kinases/physiology* ; Rats ; Myocytes, Cardiac/pathology* ; Disulfiram/pharmacology* ; Cardiomegaly ; Apoptosis/drug effects* ; Cell Line ; Angiotensin II ; Necroptosis/drug effects* ; Interleukin-1beta/metabolism* ; Receptor-Interacting Protein Serine-Threonine Kinases/metabolism* ; Lactones ; Resorcinols ; Zearalenone/administration & dosage*

ObjectiveTo establish a model that can quickly identify the aroma components in different parts of Nardostachys jatamansi， so as to provide a quality control basis for the market circulation and clinical use of N. jatamansi. MethodsHeadspace solid-phase microextraction-gas chromatography-mass spectrometry（HS-SPME-GC-MS） combined with Smart aroma database and National Institute of Standards and Technology（NIST） database were used to characterize the aroma components in different parts of N. jatamansi， and the aroma components were quantified according to relative response factor（RRF） and three internal standards， and the markers of aroma differences in different parts of N. jatamansi were identified by orthogonal partial least squares-discriminant analysis（OPLS-DA） and cluster thermal analysis based on variable importance in the projection（VIP） value >1 and P<0.01. The odor data of different parts of N. jatamansi were collected by Heracles Ⅱ Neo ultra-fast gas phase electronic nose， and the correlation between compound types of aroma components collected by the ultra-fast gas phase electronic nose and the detection results of HS-SPME-GC-MS was investigated by drawing odor fingerprints and odor response radargrams. Chromatographic peak information with distinguishing ability≥0.700 and peak area≥200 was selected as sensor data， and the rapid identification model of different parts of N. jatamansi was established by principal component analysis（PCA）， discriminant factor alysis（DFA）， soft independent modeling of class analogies（SIMCA） and statistical quality control analysis（SQCA）. ResultsThe HS-SPME-GC-MS results showed that there were 28 common components in the underground and aboveground parts of N. jatamansi， of which 22 could be quantified and 12 significantly different components were screened out. Among these 12 components， the contents of five components（ethyl isovalerate， 2-pentylfuran， benzyl alcohol， nonanal and glacial acetic acid，） in the aboveground part of N. jatamansi were significantly higher than those in the underground part（P<0.01）， the contents of β-ionone， patchouli alcohol， α-caryophyllene， linalyl butyrate， valencene， 1，8-cineole and p-cymene in the underground part of N. jatamansi were significantly higher than those in the aboveground part（P<0.01）. Heracles Ⅱ Neo electronic nose results showed that the PCA discrimination index of the underground and aboveground parts of N. jatamansi was 82， and the contribution rates of the principal component factors were 99.94% and 99.89% when 2 and 3 principal components were extracted， respectively. The contribution rate of the discriminant factor 1 of the DFA model constructed on the basis of PCA was 100%， the validation score of the SIMCA model for discrimination of the two parts was 99， and SQCA could clearly distinguish different parts of N. jatamansi. ConclusionHS-SPME-GC-MS can clarify the differential markers of underground and aboveground parts of N. jatamansi. The four analytical models provided by Heracles Ⅱ Neo electronic nose（PCA， DFA， SIMCA and SQCA） can realize the rapid identification of different parts of N. jatamansi. Combining the two results， it is speculated that terpenes and carboxylic acids may be the main factors contributing to the difference in aroma between the underground and aboveground parts of N. jatamansi.