Objective: To evaluate the current situation of thermal environment in primary and secondary school classrooms during winter, and to analyze students thermal comfort needs, so as to provide a basis for improving classroom thermal environment. Methods: From December 16 to 26, 2024, a stratified cluster random sampling method was used to select 90 classrooms from 15 primary and secondary schools in centralized/air conditioned heating areas(Liaoning Province, Tianjin City, Shanghai City) and naturally ventilated areas(Anhui Province and Jiangxi Province)for on site environmental measurement. A questionnaire survey was conducted among 743 students. The differences between groups using the χ 2 test were compared. Based on actual measurement data, a predicted mean vote prepared percentage of dissatisfied (PMV-PPD) model for centralized/air conditioned classrooms and an adaptive model for naturally ventilated classrooms were established, and the thermal neutral temperature and comfort interval were calculated. Results: The average outdoor temperature during on site measurement was 4.00(0.20，7.00)℃. In classrooms with centralized or air conditioned heating systems, the measured average temperature was (19.33±2.59)℃, with a thermal comfort range of 20.35-25.35 ℃ and a thermal neutral temperature of 22.85 ℃. And 13.92% of students reported feeling cold, while 80.80% felt comfortable. In classrooms with natural ventilation, the measured average temperature was (12.26±1.83)℃, with a thermal neutral temperature of 19.67 ℃ and a thermal comfort range of 16.17-23.17 ℃. About 48.33% of students reported feeling cold, and 49.81 % felt comfortable.The results of univariate analysis showed that there were statistically significant differences in shoe thickness, temperature sensation, relative humidity sensation and wind speed sensation between centralized/air conditioned heating areas ( χ 2= 7.01 , 31.47, 13.57, 13.80,all P <0.05). There were also statistically significant differences in school stage for primary and secondary school students, body mass index, classroom location for seat, temperature sensation, relative humidity sensation and wind speed sensation between naturally ventilated areas ( χ 2=42.13, 11.13, 11.04, 60.39, 29.27, 38.46,all P <0.05). Conclusions There are differences in thermal environment and students subjective thermal comfort in primary and secondary schools under different ventilation modes in winter. The temperature standards for heated classrooms should be revised, and differentiated environmental regulation strategies should be adopted based on different ventilation methods to improve students health and comfort levels.

ObjectiveThis paper aims to investigate and evaluate the staged efficacy and safety of the representative empirical prescription of the “Zhu Ke Jiao” theory， Qijia Rougan prescription， combined with entecavir in the treatment of hepatic fibrosis in chronic hepatitis B. MethodsA multicenter randomized controlled clinical study was conducted， and 101 patients diagnosed with chronic hepatitis B-related hepatic fibrosis （CHB-HF） who met the diagnosis and inclusion criteria were randomly assigned to an observation group （Qijia Rougan prescription + entecavir） and a control group （entecavir）. The treatment duration was 24 weeks. Liver stiffness measurement （LSM）， fibrosis-4 index （FIB-4）， portal vein diameter， hepatitis B serology， biochemical indicators， hepatic fibrosis markers in serum ［hyaluronic acid （HA）， laminin （LN）， procollagen Ⅲ peptide （PⅢP）， and type Ⅳ collagen （Ⅳ-C）］， and traditional Chinese medicine syndrome scores were used as efficacy evaluation indicators. Efficacy assessments and explorations of different staged subgroups of Qijia Rougan prescription were conducted according to LSM values based on the Metavir pathological staging standard. ResultsA total of 98 cases were included for statistical analysis， with 49 cases in the observation group and 49 in the control group. The general data of the patients in both groups were comparable. Compared with the same group before treatment， the observation group showed a significant reduction in LSM and FIB-4 （P<0.01）， as well as notable improvements in LN， Ⅳ-C， and various TCM syndrome scores （P<0.05， P<0.01）. When compared to the control group after treatment， the observation group demonstrated significant improvements in LSM， FIB-4， and various TCM syndrome score indicators （P<0.05， P<0.01）， indicating that the observation group performed better than the control group. Subgroup analysis of the regression of hepatic fibrosis stages showed that compared to the same group before treatment， the observation group had better improvement in regression of stages F2 and F3 （P<0.05）. When compared to the control group after treatment， the observation group exhibited superior improvement in regression of stage F3 （P<0.05）. No adverse events occurred in either group during the treatment period. ConclusionCompared with entecavir alone， the combination of Qijia Rougan prescription and entecavir significantly improves the degree of hepatic fibrosis and clinical TCM symptoms in patients. The optimal intervention period is primarily during stage F3， which is a potential “interception” point of the “Zhu Ke Jiao” theory.

Metabolic dysfunction-associated steatotic liver disease (MASLD) has become the most prevalent chronic liver disease worldwide, affecting approximately 32%-38% of the adult population and posing a growing public health burden. MASLD represents a continuous disease spectrum ranging from simple steatosis to metabolic dysfunction-associated steatohepatitis (MASH), progressive hepatic fibrosis, cirrhosis, and ultimately hepatocellular carcinoma (HCC). The pathological core of MASLD lies in disruption of hepatic lipid metabolic homeostasis, characterized by an imbalance among de novo lipogenesis, fatty acid β-oxidation, and very-low-density lipoprotein (VLDL)-mediated lipid export. This metabolic disequilibrium subsequently drives inflammatory injury and fibrotic progression. Among the multiple regulatory pathways involved, thyroid hormone (TH) signaling has emerged as a central regulator of hepatic metabolic homeostasis. The liver is a major peripheral target organ of TH action, where TH predominantly exerts its metabolic effects through thyroid hormone receptor β (TRβ). Large-scale epidemiological studies and meta-analyses have demonstrated that hypothyroidism is significantly associated with increased MASLD prevalence, more severe histological injury, and advanced hepatic fibrosis, suggesting that dysregulation of TH signaling may participate throughout the entire MASLD disease spectrum. At the molecular level, TH regulates hepatic lipid metabolism by coordinating suppression of lipogenesis, enhancement of mitochondrial fatty acid oxidation, and promotion of VLDL assembly and secretion through integrated genomic actions of the T3-TRβ axis and non-genomic signaling pathways. Across different stages of MASLD, TH signaling exerts stage-dependent protective effects. In the steatosis stage, TH improves metabolic flexibility by modulating insulin sensitivity, glucose metabolism, and lipid droplet clearance, thereby alleviating early lipotoxic stress. During progression to MASH, TH attenuates inflammatory amplification by improving mitochondrial homeostasis, suppressing activation of the NOD-like receptor family pyrin domain containing 3 (NLRP3) inflammasome, and modulating the gut-liver axis microenvironment. In advanced stages, TH signaling influences hepatic stellate cell activation and extracellular matrix deposition, partly through interaction with the transforming growth factor-β (TGF-β)/SMAD pathway, while alterations in intrahepatic TH availability, mediated by dynamic changes in iodothyronine deiodinase 1 (DIO1), contribute to fibrosis progression and hepatocellular dedifferentiation. In hepatocellular carcinoma, coordinated downregulation of TRβ and DIO1 establishes a tumor-associated hypothyroid state that promotes metabolic reprogramming and tumor progression. The clinical relevance of TH signaling in MASLD has been underscored by the recent approval of Resmetirom, a liver-targeted TRβ‑selective agonist, for the treatment of non-cirrhotic MASH with moderate-to-severe fibrosis (F2-F3). This approval represents a landmark transition from mechanistic understanding to metabolism-centered precision therapy in MASLD. Clinical trials have demonstrated that Resmetirom not only improves key histological endpoints, including MASH resolution and fibrosis regression, but also favorably modulates atherogenic lipid profiles, highlighting the therapeutic potential of selectively targeting hepatic TH pathways. This review systematically summarizes the multidimensional regulatory roles of TH across the MASLD disease spectrum and discusses emerging diagnostic and therapeutic implications of TH-based interventions, aiming to inform future mechanistic research and optimize clinical management strategies.

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.

OBJECTIVE To compare the efficacy and safety of vedolizumab （VDZ） versus infliximab （IFX） in biologic-naive（Bio-naive） patients with moderate-to-severe active ulcerative colitis （UC）， and to analyze the factors influencing efficacy. METHODS Clinical data were retrospectively collected from Bio-naive patients with moderate-to-severe active UC who received treatment in the Department of Gastroenterology at Shanxi Provincial People’s Hospital from June 2023 to June 2024. Based on the type of biologic agent administered， the patients were divided into the IFX group （41 cases） and the VDZ group （30 cases）. Patients in the two groups received IFX （5 mg/kg） or VDZ （300 mg） for induction and maintenance of remission therapy. The two groups were compared regarding modified Mayo score， serological indicators （hemoglobin， albumin， platelet count， erythrocyte sedimentation rate， C-reactive protein）， combined medication， efficacy-related indexes （clinical response rate/remission rate， and endoscopic response rate/remission rate）， and the occurrence of adverse drug reactions （ADR）. Based on Logistic regression model， univariate and multivariate analyses were conducted to identify potential factors influencing clinical remission at week 14 and endoscopic remission at week 38. RESULTS There were no statistically significant differences in clinical response rate/remission rate， or endoscopic response rate/remission rate between the two groups at weeks 14 and 38 （ P ＞0.05）. However， at week 14 of treatment， the proportion of patients using concomitant corticosteroids in VDZ group was 26.67%， significantly higher than the 7.50% in IFX group （ P ＜0.05）. There was no statistical significance in the overall incidence of ADR between the two groups （ P ＞0.05）； all ADRs in the IFX group were grade 3 and led to treatment discontinuation （6 cases）， whereas ADR in the VDZ group was grade 2 and did not interrupt therapy （1 case）. Univariate and multivariate regression analyses revealed that disease type （relapsing） was significantly associated with clinical remission at week 14 of treatment， and a history of smoking was significantly associated with endoscopic remission at week 38 of treatment （the odds ratios were 0.08 for both， with 95% confidence intervals of 0.01-0.77 and 0.01-0.91 respectively， P ＜0.05）. CONCLUSIONS For Bio-naive patients with moderate-to-severe active UC， VDZ and IFX demonstrate comparable efficacy in inducing and maintaining clinical remission and promoting mucosal healing， as well as overall safety. Although IFX can achieve faster control of inflammation in the early stage of the disease， it causes more severe ADR. Disease type （relapsing） and smoking history are identified as independent negative predictors for short-term clinical remission and long-term endoscopic remission， respectively.

Licochalcone E is a chalcone isolated from Glycyrrhiza uralensis and G. inflata Batal. This study explored the effect of licochalcone E on improving hepatic steatosis in obese mice and evaluated the role of licochalcone E in regulating lipid accumulation in hepatocytes. In vitro, oleic acid–induced hepatocytes were treated with licochalcone E to investigate its effect on lipid metabolic pathways. In animal experiments, male C57BL/6 mice were fed with a high-fat diet (HFD) and treated with licochalcone E by intraperitoneal injection for 12 weeks to assess its effects on biochemical indexes and hepatic steatosis. Furthermore, mice were fed a methionine/choline-deficient (MCD) diet and administered licochalcone E, followed by evaluation of liver fibrosis. Licochalcone E effectively reduced body weight, epididymal and inguinal fat weight, and adipocyte size in HFD-induced obese mice. Licochalcone E treatment of obese mice also reduced hepatic lipid accumulation and improved hepatocyte steatosis. Licochalcone E regulated the expression of lipogenesis- and lipolysis-related genes in the livers of obese mice and increased AMPK phosphorylation and Sirt1 expression in the liver. Licochalcone E also attenuated hepatic inflammation and oxidative stress in obese mice. Furthermore, treatment of MCD-induced mice with licochalcone E reduced the number of lipid vacuoles and the extent of fibrosis and inhibited liver inflammation. In FL83B hepatocytes, licochalcone E could regulate lipogenesis and lipolysis, and increase the phosphorylation of AMPK and ACC. These findings provide new insights into the role of licochalcone E in regulating lipid metabolism and preventing hepatic steatosis.

ObjectiveTo investigate the association of the traditional Chinese medicine （TCM） syndrome types of primary liver cancer （PLC） with clinical features and multimodal quantitative radiological features on computed tomography （CT） and magnetic resonance imaging （MRI）， and to provide a reference for the objectification of TCM syndrome differentiation and precise diagnosis and treatment. MethodsA retrospective analysis was performed for the clinical data of 312 patients who were diagnosed with PLC in The First Affiliated Hospital of Hunan University of Chinese Medicine from March 2020 to June 2025， and according to the TCM syndrome type， they were divided into stagnation of liver Qi group with 40 patients， stagnation of liver Qi and spleen deficiency group with 109 patients， Qi stagnation and blood stasis group with 62 patients， dampness-heat toxin amassment group with 81 patients， and liver-kidney Yin deficiency group with 20 patients. Clinical features and multimodal quantitative radiological features were compared between the patients with different TCM syndrome types. A one-way analysis of variance was used for comparison of normally distributed continuous data between multiple groups， and the least significant difference t-test was used for further comparison between two groups； the Kruskal-Wallis H test was used for comparison of non-normally distributed continuous data between multiple groups， and the Dunn’s multiple test was used for further comparison between two groups； the chi-square test was used for comparison of categorical data between groups， and the Bonferroni method was used for further comparison between two groups. ResultsThere were significant differences between the patients with different TCM syndrome types in China liver cancer staging （CNLC）， Child-Pugh class， alanine aminotransferase， aspartate aminotransferase， albumin， direct bilirubin， total bilirubin， prothrombin time， neutrophil， and albumin-bilirubin score （all P<0.05）. In the stagnation of liver Qi group， the patients with Child-Pugh class A accounted for 75.00%； among the patients with CNLC stage I PLC， the patients with stagnation of liver Qi accounted for 60.00%， and those with Qi stagnation and blood stasis syndrome accounted for 59.68%， while among the patients with CNLC stage IV PLC， the distribution proportion of dampness-heat toxin amassment （27.16%） and liver-kidney Yin deficiency （30.00%） was significantly higher than that of stagnation of liver Qi （2.50%） （all P<0.05）. Radiological examination showed that there were significant differences between the patients with different TCM syndrome types in the number of tumors， ascites， venous tumor thrombus， maximum tumor diameter， intrahepatic metastasis， and lymph node metastasis in the hepatic hilar and retroperitoneal regions （all P<0.05）. Compared with the patients with stagnation of liver Qi， the patients with liver depression and spleen deficiency or liver-kidney Yin deficiency were more likely to develop intrahepatic metastasis； the patients with liver depression and spleen deficiency， dampness-heat toxin amassment， or liver-kidney Yin deficiency were more likely to develop lymph node metastasis in the hepatic hilar and retroperitoneal regions； the patients with liver-kidney Yin deficiency were more likely to experience multiple tumors； the patients with liver depression and spleen deficiency or dampness-heat toxin amassment were more likely to develop ascites （all P<0.05）. Compared with the patients with Qi stagnation and blood stasis syndrome， the patients with liver depression and spleen deficiency had a significantly longer maximum tumor diameter and a significantly higher proportion of patients with venous tumor thrombus （both P<0.05）. Furthermore， among the 184 patients with MRI diffusion-weighted imaging sequences， the patients with dampness-heat toxin amassment or Qi stagnation and blood stasis syndrome had significantly higher ADC values and relative ADC values than those with stagnation of liver Qi （all P<0.05）. ConclusionThere are significant differences in CT/MRI radiological features and clinical features between PLC patients with different TCM syndrome types， among whom the patients with liver depression and spleen deficiency， dampness-heat toxin amassment， and liver-kidney Yin deficiency tend to exhibit progressive radiological features， and those with dampness-heat toxin amassment or Qi stagnation and blood stasis syndrome tend to have higher ADC values. These findings provide an objective basis for TCM syndrome differentiation in PLC.

Objective This study investigated the diversity of rodent populations and the infection status of Bartonella at three ports along the China-Mongolia border in 2022. Methods Rodents at three Belt and Road ports along the China-Mongolia border, namely Ganqimaodu, Erenhot, and Zhuengadabuqi, were morphologically identified. The Bartonella citrate synthase(gltA)gene was amplified by nested polymerase chain reaction(PCR), and PCR-positive products were sequenced. The resulting sequences were then analyzed for genetic characteristics. Phylogenetic analysis was performed using MEGA 11.0 software using Neighbor-Joining and Maximum-Likelihood method. Results A total of 94 rodents were captured, representing seven species from five families and six genera: Mus musculus, Meriones unguiculatus, Meriones meridianus, Spermophilus dauricus, Allactaga sibirica, Dipus sagitta, and Phodopus roborovskii. Among them, M. unguiculatus was the most abundant species, with 58 rodents, accounting for 61.70% of the total. Overall,11 positive Bartonella pathogen sequences were obtained from the three ports, with a positive detection rate of 11.70%(11/94). The infected rodent species included M. unguiculatus, M. meridianus, and A. sibirica, and two species of Bartonella were detected. Conclusions Rodents at Ganqimaodu, Erenhot, and Zhuengadabuqi ports along the China-Mongolia border were naturally infected with Bartonella. Rodent monitoring and pathogen prevention and control in this area should be strengthened.

Objective This study reports on the obligatory parasitism of water mites Arrenurus sp. on Anopheles messeae at the Manzhouli Port, Inner Mongolia, China. Methods Duing July 2024, a survey on the mosquito diversity was conducted at the Manzhouli Port. Captured mosquitoes and their ectoparasites were identified to species level. Results A total of 1840 adult mosquitoes were collected, representing species from three genera: Culex(Cx. modestus, Cx. pipiens pallens), Aedes(Ae. dorsalis, Ae. flavidorsalis, Ae. flavescens), and Anopheles (An. messeae). Among all the mosqutioes specimens,3 out of 150 captured An. messeae were found to carry ectoparasitic mites, with number of 2,4,27 mites separately. Morphological and molecular identification reached the same result as water mites(Hydrachnidiae, Hydracrina). COI gene sequence showed 94% similarity with the closest species Arrenurus truncatellus. Conlusions Literature review suggests water mites are host-specific parasitism of mosquito species and herein with the first record of Arrenurus sp. parasiting on An. Messeae in the most high-latitude region globally.