ObjectiveSpinal cord injury (SCI) directly impairs the regulatory function of the autonomic nervous system, induces intestinal dysfunction, and significantly reduces patients’ quality of life. Preclinical studies have shown that electroacupuncture (EA) therapy can regulate the brain-gut axis and is used to treat central nervous system diseases such as major depressive disorder, Alzheimer’s disease and Parkinson’s disease. Recent research has established that fecal microbiota transplantation (FMT) from EA-treated SCI rats restored intestinal motility and colonic morphology. However, it remains unclear whether the regulation of gut microbiota by EA therapy directly contributes to neural repair after SCI. This study aims to explore whether gut microbiota mediates the neuroprotective effect of EA in the treatment of SCI and its possible mechanism. MethodsThe study employed RNA transcriptome analysis of spinal cord tissue to characterize gene expression profiles and to identify key signaling pathways following EA treatment for SCI. Hematoxylin-Eosin (HE) staining and Nissl staining were used to observe the morphological changes in spinal cord tissue. Western blot (WB) and enzyme-linked immunosorbent assay (ELISA) were applied to detect the effects of EA on the expression of proteins related to nucleotide-binding domain leucine-rich repeat and pyrin domain-containing receptor 3 (NLRP3) -dependent pyroptosis. Using 16S rDNA sequencing, the study observed alterations in gut microbiota diversity and community composition in SCI rats. Prior to establishing SCI models, rats were pretreated with an antibiotic cocktail to induce gut dysbiosis, and the effects on intestinal function and spinal cord neural repair were evaluated. FMT was performed to investigate the regulatory effects of post-EA FMT on motor function, general status, liver and spleen indices, and NLRP3-mediated pyroptosis in SCI rats. ResultsEA improved motor function and reduced regulated neuronal cell death in SCI rats. Transcriptomic analysis demonstrated the activation of immune- and inflammation-related pathways post-SCI, including NOD-like receptors, nuclear factor-kappa B(NF-κB), and Toll-like receptor (TLR) pathways. EA primarily influenced intestinal inflammation and autoimmune functions. 16S rDNA sequencing illustrated that EA did not alter the diversity of gut microbiota. However, EA altered the gut microbiota composition in SCI rats, increasing Lactobacillus and Akkermansia genera while rebalancing the Firmicutes/Bacteroidetes ratio. Furthermore, depletion of gut microbiota by antibiotics disrupted the intestinal barrier, reduced the expression of intestinal barrier proteins Zonula Occludens-1 (ZO-1) and Occludin, elevated serum lipopolysaccharide-binding protein (LBP) levels, exacerbated spinal cord tissue damage, and hindered motor function recovery in SCI rats. FMT from donors treated with EA reduced LBP levels in the intestine, blood, and spinal cord of rats, inhibited the TLR4 myeloid differentiation primary response protein 88 (MyD88)-NF‑κB pathway and NLRP3-dependent pyroptosis, and improved motor function. On the other hand, FMT treatment resulted in decreased body weight and food intake, whereas FMT using EA-treated donors effectively alleviated these alterations. ConclusionEA effectively alleviated neuroinflammatory responses in rats with SCI, primarily through regulating the gut microbiota and suppressing the NLRP3-dependent pyroptosis signaling pathway.

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.

ObjectiveSpinal cord injury (SCI) directly impairs the regulatory function of the autonomic nervous system, induces intestinal dysfunction, and significantly reduces patients’ quality of life. Preclinical studies have shown that electroacupuncture (EA) therapy can regulate the brain-gut axis and is used to treat central nervous system diseases such as major depressive disorder, Alzheimer’s disease and Parkinson’s disease. Recent research has established that fecal microbiota transplantation (FMT) from EA-treated SCI rats restored intestinal motility and colonic morphology. However, it remains unclear whether the regulation of gut microbiota by EA therapy directly contributes to neural repair after SCI. This study aims to explore whether gut microbiota mediates the neuroprotective effect of EA in the treatment of SCI and its possible mechanism. MethodsThe study employed RNA transcriptome analysis of spinal cord tissue to characterize gene expression profiles and to identify key signaling pathways following EA treatment for SCI. Hematoxylin-Eosin (HE) staining and Nissl staining were used to observe the morphological changes in spinal cord tissue. Western blot (WB) and enzyme-linked immunosorbent assay (ELISA) were applied to detect the effects of EA on the expression of proteins related to nucleotide-binding domain leucine-rich repeat and pyrin domain-containing receptor 3 (NLRP3) -dependent pyroptosis. Using 16S rDNA sequencing, the study observed alterations in gut microbiota diversity and community composition in SCI rats. Prior to establishing SCI models, rats were pretreated with an antibiotic cocktail to induce gut dysbiosis, and the effects on intestinal function and spinal cord neural repair were evaluated. FMT was performed to investigate the regulatory effects of post-EA FMT on motor function, general status, liver and spleen indices, and NLRP3-mediated pyroptosis in SCI rats. ResultsEA improved motor function and reduced regulated neuronal cell death in SCI rats. Transcriptomic analysis demonstrated the activation of immune- and inflammation-related pathways post-SCI, including NOD-like receptors, nuclear factor-kappa B(NF-κB), and Toll-like receptor (TLR) pathways. EA primarily influenced intestinal inflammation and autoimmune functions. 16S rDNA sequencing illustrated that EA did not alter the diversity of gut microbiota. However, EA altered the gut microbiota composition in SCI rats, increasing Lactobacillus and Akkermansia genera while rebalancing the Firmicutes/Bacteroidetes ratio. Furthermore, depletion of gut microbiota by antibiotics disrupted the intestinal barrier, reduced the expression of intestinal barrier proteins Zonula Occludens-1 (ZO-1) and Occludin, elevated serum lipopolysaccharide-binding protein (LBP) levels, exacerbated spinal cord tissue damage, and hindered motor function recovery in SCI rats. FMT from donors treated with EA reduced LBP levels in the intestine, blood, and spinal cord of rats, inhibited the TLR4 myeloid differentiation primary response protein 88 (MyD88)-NF‑κB pathway and NLRP3-dependent pyroptosis, and improved motor function. On the other hand, FMT treatment resulted in decreased body weight and food intake, whereas FMT using EA-treated donors effectively alleviated these alterations. ConclusionEA effectively alleviated neuroinflammatory responses in rats with SCI, primarily through regulating the gut microbiota and suppressing the NLRP3-dependent pyroptosis signaling pathway.

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 innate immune system can be boosted in response to subsequent triggers by pre-exposure to microbes or microbial products, known as “trained immunity”. Compared to classical immune memory, innate trained immunity has several different features. Firstly, the molecules involved in trained immunity differ from those involved in classical immune memory. Innate trained immunity mainly involves innate immune cells (e.g., myeloid immune cells, natural killer cells, innate lymphoid cells) and their effector molecules (e.g., pattern recognition receptor (PRR), various cytokines), as well as some kinds of non-immune cells (e.g., microglial cells). Secondly, the increased responsiveness to secondary stimuli during innate trained immunity is not specific to a particular pathogen, but influences epigenetic reprogramming in the cell through signaling pathways, leading to the sustained changes in genes transcriptional process, which ultimately affects cellular physiology without permanent genetic changes (e.g., mutations or recombination). Finally, innate trained immunity relies on an altered functional state of innate immune cells that could persist for weeks to months after initial stimulus removal. An appropriate inducer could induce trained immunity in innate lymphocytes, such as exogenous stimulants (including vaccines) and endogenous stimulants, which was firstly discovered in bone marrow derived immune cells. However, mature bone marrow derived immune cells are short-lived cells, that may not be able to transmit memory phenotypes to their offspring and provide long-term protection. Therefore, trained immunity is more likely to be relied on long-lived cells, such as epithelial stem cells, mesenchymal stromal cells and non-immune cells such as fibroblasts. Epigenetic reprogramming is one of the key molecular mechanisms that induces trained immunity, including DNA modifications, non-coding RNAs, histone modifications and chromatin remodeling. In addition to epigenetic reprogramming, different cellular metabolic pathways are involved in the regulation of innate trained immunity, including aerobic glycolysis, glutamine catabolism, cholesterol metabolism and fatty acid synthesis, through a series of intracellular cascade responses triggered by the recognition of PRR specific ligands. In the view of evolutionary, trained immunity is beneficial in enhancing protection against secondary infections with an induction in the evolutionary protective process against infections. Therefore, innate trained immunity plays an important role in therapy against diseases such as tumors and infections, which has signature therapeutic effects in these diseases. In organ transplantation, trained immunity has been associated with acute rejection, which prolongs the survival of allografts. However, trained immunity is not always protective but pathological in some cases, and dysregulated trained immunity contributes to the development of inflammatory and autoimmune diseases. Trained immunity provides a novel form of immune memory, but when inappropriately activated, may lead to an attack on tissues, causing autoinflammation. In autoimmune diseases such as rheumatoid arthritis and atherosclerosis, trained immunity may lead to enhance inflammation and tissue lesion in diseased regions. In Alzheimer’s disease and Parkinson’s disease, trained immunity may lead to over-activation of microglial cells, triggering neuroinflammation even nerve injury. This paper summarizes the basis and mechanisms of innate trained immunity, including the different cell types involved, the impacts on diseases and the effects as a therapeutic strategy to provide novel ideas for different diseases.

Alzheimer’s disease (AD) is a chronic neurodegenerative disorder that severely affects the health of the elderly, marked by its incurability, high prevalence, and extended latency period. The current approach to AD prevention and treatment emphasizes early detection and intervention, particularly during the pre-AD stage of mild cognitive impairment (MCI), which provides an optimal “window of opportunity” for intervention. Clinical detection methods for MCI, such as cerebrospinal fluid monitoring, genetic testing, and imaging diagnostics, are invasive and costly, limiting their broad clinical application. Speech, as a vital cognitive output, offers a new perspective and tool for computer-assisted analysis and screening of cognitive decline. This is because elderly individuals with cognitive decline exhibit distinct characteristics in semantic and audio information, such as reduced lexical richness, decreased speech coherence and conciseness, and declines in speech rate, voice rhythm, and hesitation rates. The objective presence of these semantic and audio characteristics lays the groundwork for computer-based screening of cognitive decline. Speech information is primarily sourced from databases or collected through tasks involving spontaneous speech, semantic fluency, and reading, followed by analysis using computer models. Spontaneous language tasks include dialogues/interviews, event descriptions, narrative recall, and picture descriptions. Semantic fluency tasks assess controlled retrieval of vocabulary items, requiring participants to extract information at the word level during lexical search. Reading tasks involve participants reading a passage aloud. Summarizing past research, the speech characteristics of the elderly can be divided into two major categories: semantic information and audio information. Semantic information focuses on the meaning of speech across different tasks, highlighting differences in vocabulary and text content in cognitive impairment. Overall, discourse pragmatic disorders in AD can be studied along three dimensions: cohesion, coherence, and conciseness. Cohesion mainly examines the use of vocabulary by participants, with a reduction in the use of nouns, pronouns, verbs, and adjectives in AD patients. Coherence assesses the ability of participants to maintain topics, with a decrease in the number of subordinate clauses in AD patients. Conciseness evaluates the information density of participants, with AD patients producing shorter texts with less information compared to normal elderly individuals. Audio information focuses on acoustic features that are difficult for the human ear to detect. There is a significant degradation in temporal parameters in the later stages of cognitive impairment; AD patients require more time to read the same paragraph, have longer vocalization times, and produce more pauses or silent parts in their spontaneous speech signals compared to normal individuals. Researchers have extracted audio and speech features, developing independent systems for each set of features, achieving an accuracy rate of 82% for both, which increases to 86% when both types of features are combined, demonstrating the advantage of integrating audio and speech information. Currently, deep learning and machine learning are the main methods used for information analysis. The overall diagnostic accuracy rate for AD exceeds 80%, and the diagnostic accuracy rate for MCI also exceeds 80%, indicating significant potential. Deep learning techniques require substantial data support, necessitating future expansion of database scale and continuous algorithm upgrades to transition from laboratory research to practical product implementation.

The current study aimed to clarify the roles of apolipoprotein A5 (ApoA5) and milk fat globule-epidermal growth factor 8 (Mfge8) in regulating myocardial lipid deposition and the regulatory relationship between them. The serum levels of ApoA5 and Mfge8 in obese and healthy people were compared, and the obesity mouse model induced by the high-fat diet (HFD) was established. In addition, primary cardiomyocytes were purified and identified from the hearts of suckling mice. The 0.8 mmol/L sodium palmitate treatment was used to establish the lipid deposition cardiomyocyte model in vitro. ApoA5-overexpressing adenovirus was used to observe its effects on cardiac function and lipids. The expressions of the fatty acid uptake-related molecules and Mfge8 on transcription or translation levels were detected. Co-immunoprecipitation was used to verify the interaction between ApoA5 and Mfge8 proteins. Immunofluorescence was used to observe the co-localization of Mfge8 protein with ApoA5 or lysosome-associated membrane protein 2 (LAMP2). Recombinant rMfge8 was added to cardiomyocytes to investigate the regulatory mechanism of ApoA5 on Mfge8. The results showed that participants in the simple obesity group had a significant decrease in serum ApoA5 levels (P < 0.05) and a significant increase in Mfge8 levels (P < 0.05) in comparison with the healthy control group. The adenovirus treatment successfully overexpressed ApoA5 in HFD-fed obese mice and palmitic acid-induced lipid deposition cardiomyocytes, respectively. ApoA5 reduced the weight of HFD-fed obese mice (P < 0.05), shortened left ventricular isovolumic relaxation time (IVRT), increased left ventricular ejection fraction (LVEF), and significantly reduced plasma levels of triglycerides (TG) and cholesterol (CHOL) (P < 0.05). In myocardial tissue and cardiomyocytes, the overexpression of ApoA5 significantly reduced the deposition of TG (P < 0.05), transcription of fatty acid translocase (FAT/CD36) (P < 0.05), fatty acid-binding protein (FABP) (P < 0.05), and fatty acid transport protein (FATP) (P < 0.05), and protein expression of Mfge8 (P < 0.05), while the transcription levels of Mfge8 were not significantly altered (P > 0.05). In vitro, the Mfge8 protein was captured using ApoA5 as bait protein, indicating a direct interaction between them. Overexpression of ApoA5 led to an increase in co-localization of Mfge8 with ApoA5 or LAMP2 in cardiomyocytes under lipid deposition status. On this basis, exogenous added recombinant rMfge8 counteracted the improvement of lipid deposition in cardiomyocytes by ApoA5. The above results indicate that the overexpression of ApoA5 can reduce fatty acid uptake in myocardial cells under lipid deposition status by regulating the content and cellular localization of Mfge8 protein, thereby significantly reducing myocardial lipid deposition and improving cardiac diastolic and systolic function.
Animals ; Humans ; Mice ; Myocytes, Cardiac/metabolism* ; Obesity/physiopathology* ; Male ; Apolipoprotein A-V/blood* ; Lipid Metabolism/physiology* ; Milk Proteins/blood* ; Myocardium/metabolism* ; Diet, High-Fat ; Antigens, Surface/physiology* ; Mice, Inbred C57BL ; Cells, Cultured ; Female

AIM To establish the quantitative models for gallic acid,mononucleoside,loganin,resveratrol,and rhein in Yishen Xiezhuo Mixture.METHODS HPLC was adopted in the content determination of various effective components,after which the near-infrared spectroscopy(NIRS)data were collected in 128 batches of samples and pretreatment was conducted,competitive adaptive reweighting sampling(CARS)algorithm was used for screening wavelength,partial least square method(PLS)regression analysis was performed.RESULTS There were no significant differences between the predicted values obtained by PLS models and measured values obtained by HPLC for various effective components(P＞0.05).CONCLUSION The quantitative models established by NIRS combined with chemometrics display good predictive performance,which can be used for the rapid determination of effective components in Yishen Xiezhuo Mixture,and provide a reference for the rapid monitoring of other traditional Chinese medicine preparations in production processes.

This study was analyzed the spatial-temporal distribution and aggregation characteristics of Yersinia pestispositive host animals and vector pathogens in marmot natural foci in the Qilian-Altun mountains,Gansu Province,to provide a scientific basis for precise plague prevention and control.Y.pestissurveillance data for marmot natural foci in Qilian-Altun Mountains of Gansu Province from 2014 to 2023 were obtained from the Disease Control and Prevention Center of Gansu Province.Origin 2024 software was used for data visualization and presentation.Global and local spatial autocorrelation analyses and trend analyses were conducted in ArcGIS 10.8 software,with townships as the spatial scale.Cumulatively,440 strains of Y.pestis were isolated from the natural marmot foci in the Qilian-Altun mountainsof Gansu Province from 2014 to 2023.Most strains was isolated from marmots(345 strains,78.41%),and the remainder were isolated from vectors.Temporal distribution analysis indicated that the highest number of detected bacteria was reported in July and August(both 121 strains,27.50%).Regional distribution analysis revealed that Aksai County reported the highest number of detected bacteria(255 strains,57.95%).Global spatial autocorrelation analysis showed a spatially clustered distribution of the number of bacteria detected annually in the townships containing natural foci,except in2014,2016,and 2021-2023.The strongest spatial clustering was observed in 2020(Moran's I=0.521 2,Z=14.397 0,P<0.001).Local spatial autocorrelation analysis indicated a"high-high"aggregation area in the natural foci every year from 2014 to 2023,primarily in Hongliuwan Town of Aksai County and Dangchengwan Town of Subei County.The distribution of the"low-low"aggregation area was essentially consistent with the low activity area of the Yersinia pestisepidemic.The trend in annual total bacterial count gradually increased from east to west,and peaked in the western part of the epidemic focus.Clear spatial aggregation characteristics of the number of Y.pestis were detected in the marmot natural foci in the Qilian-Altun mountains at the townshiplevel as a whole in Gansu Province from 2014 to 2023.The aggregation area was mainly in the western section of Qilian Mountain to the Altun mountain section of the epidemic source area.Monitoring and prevention and control efforts should be focused in this key area,with prevention and control measures tailored to the local conditions,and classified guidance to decrease the risk of plague occurrence and spread.

This study was aimed at studying the drug resistance,serotypes,and molecular characteristics of Vibrio parahaemo-lyticus(VP)in Suzhou,to provide basic data for the prevention and control of VP-related diseases.Drug susceptibility testing of 177 VP strains isolated from Suzhou City in 2023 was performed with the microbroth dilution method.Virulence genes,serotypes,and multi-locus sequence typing(MLST)were analyzed on the basis of whole genome sequencing results.The drug resistance rate of 177 VP strains was highest against cefazolin(100.00%),followed by ampicillin(77.97%),and polymyxine E(63.84%),and the multiple drug resistance rate was 53.67%.In clinical isolates,O10∶K4(37.41%)was the most abundant serotype,and was followed by O3∶K6(28.78%),and ST3 was the dominant ST type.The main virulence genes of clinical isolates were tlh+,tdh+,and trh-(79.86%),whereas the virulence genes in food isolates were all tlh+,tdh-,and trh-.Strains of the same serotype clustered together in the SNP phylogenetic tree.The environmental isolates showed no obvious dominant serotype or ST type.In Suzhou,VP has a high proportion of multi-drug resistance,the clinical isolates have prevalent serotypes and ST types,and most isolates carried virulence genes;there-fore,monitoring should be strengthened.