OBJECTIVE To form an expert consensus addressing clinical issues regarding the use of parenteral direct thrombin inhibitors （DTIs） in special populations. METHODS Led by the Sichuan Academy of Medical Sciences & Sichuan Provincial People’s Hospital（the Affiliated Hospital of UESTC）， a multidisciplinary working group was formed comprising experts from multiple fields， including clinical pharmacy， cardiac surgery， obstetrics， pediatrics and evidence-based medicine. Through literature review and the Delphi method， clinical questions regarding the efficacy and safety of parenteral DTIs used in special populations were identified. A structured design was adopted using the “Population-Intervention-Comparison-Outcome” （PICO） framework；systematic searches were conducted in CJFD， PubMed， Embase and other databases. Relevant evidence from randomized controlled trials，cohort studies and systematic reviews were included and synthesized. Evidence quality was assessed using the Grading of Recommendations Assessment，Development and Evaluation （GRADE） approach， and recommendations were formulated through three rounds of Delphi surveys and expert consensus meetings. RESULTS &CONCLUSIONS Seven clinical questions were ultimately selected （with a consensus rate exceeding 90%）， resulting in the formulation of seven recommendations on the use of parenteral DTIs in special populations， including children， pregnant women， patients with hepatic or renal impairment， patients with mesenteric venous thrombosis， and individuals with thrombophilia. These recommendations clarify the preferred agents， dosing ranges， monitoring parameters， and safety management strategies for parenteral DTIs in these special populations. This expert consensus， which is formulated based on the best available evidence， provides evidence-based guidance for standardized and individualized use of parenteral DTIs in special populations.

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.

ObjectiveThis paper aims to observe the syndrome improvement and negative antigen conversion rate of Qingxuan Daozhi formula in the treatment of influenza in children （heat accumulation in the lung and stomach syndrome）. MethodsThrough a multi-center randomized controlled methodology design，confirmed influenza cases were collected from October 2022 to April 2023 in the pediatrics department of eight hospitals，such as Dongfang Hospital of Beijing University of Chinese Medicine. A total of 180 children with influenza and heat accumulation in the lung and stomach syndrome conforming to the standard were recruited through the clinic. The sick children meeting the inclusion criteria were randomly divided into groups by a block-randomized method. The children in the experimental group were treated with Qingxuan Daozhi formula for five days，and those in the control group were treated with Oseltamivir Phosphate Granules for five days. The primary efficacy indicator was the negative conversion rate of influenza antigen detection. Secondary efficacy indicators were the Canadian acute respiratory illness and flu scale （CARIFS） and the incidence of complications，severe cases， and critical cases. Follow-up observation was conducted on the day of enrollment，48 hours after medication，72 hours after medication， and （6+1） d after medication. ResultsOne hundred and eighty participants were randomly assigned to the experimental group （90 cases） or the control group （90 cases）. All participants were followed up during the study. Comparison of influenza antigen detection results in the primary efficacy indicators showed that the average time of negative influenza antigen conversion in the experimental group was （5.29±1.25） d，and that in the control group was （5.40±1.68） d，without a statistically significant difference. After five days of intervention，52 cases in the experimental group and 51 cases in the control group converted to negative，without a statistically significant difference. CARIFS score results in the secondary efficacy indicators showed that during 72 hours after intervention，there were statistically significant differences between the experimental group and the control group in three dimensions， including headache，muscle soreness， and the need for extra care （P<0.05）. On the （6+1） days after the intervention，the differences in both the experimental group and the control group were statistically significant in 10 dimensions， including sore throat，bad sleep，uncomfortable feeling，poor spirit and fatigue，crying more than usual，the need for extra care，symptom，function，influence on parents，and total score （P<0.05）. The comparison results within the group in the dimensional scores of symptom， function， and influence on parents，as well as the CARIFS total score showed that with the delay of follow-up time，scores of both groups decreased significantly，with a statistically significant difference （P<0.01）. Inter-group comparison results showed that the mean score of the experimental group was higher than that of the control group at the time of enrollment. With the progress of intervention，the score of the experimental group was significantly decreased compared with that of the control group. At the end of follow-up，the mean score of the experimental group was lower than that of the control group，with no statistically significant difference. In terms of the incidence of complications，severe cases， and critical cases， there were no complications，severe cases， and critical cases in the two groups，without a statistically significant difference. ConclusionThe symptom improvement effect and negative antigen conversion rate of Qingxuan Daozhi formula in the treatment of influenza in children （heat accumulation in the lung and stomach syndrome） are not inferior to Oseltamivir Phosphate granules， and children's acceptance is better. It can be more widely used in clinical treatment of influenza in children （heat accumulation in the lung and stomach syndrome）.

[Objective] To investigate the seroprevalence characteristics of hepatitis E virus (HEV) among blood donors with hepatitis B virus (HBV) infection, so as to provide data support for the monitoring, prevention and treatment of HEV. [Methods] From January to December 2022, 219 samples positive for hepatitis B surface antigen (HBsAg), 142 occult hepatitis B virus infection (OBI) samples (HBV group) and 873 samples tested negative (control group) were collected. 361 samples were further tested with viral load assay and serological testing for five serological markers (HBsAg, HBsAb, HBeAg, HBeAb and HBcAb), and the DNA load was measured using real time fluorescence quantitative PCR. Commercially available enzyme-linked immunosorbent assays (ELISA) were used for the detection of anti-HEV IgG, anti-HEV IgM and HEV antigen (Ag). The Chi-square test or Fisher's exact test was used to assess the differences in the reactivity rates of anti-HEV IgG and anti-HEV IgM among different blood donor populations and different variables. Multivariable logistic regression was used to examine potential risk factors associated with anti-HEV IgG seroprevalence. [Results] In the HBV group, HBsAg positive donors exhibited low expression of antigen. The HBV DNA load of OBI infected donors ranged from 1 to 131.43 IU/mL (median 11.24 IU/mL). The prevalence of anti-HEV IgG and IgM antibody in the HBV group were 34.63% and 1.11%, respectively. Among them, the prevalence of anti-HEV IgG and anti-HEV IgM in the HBV group was 34.63% and 0, respectively (P<0.05), while in the OBI donors, they were 41.55% and 2.82%, respectively. In the normal donors, the reactivity rates for anti-HEV IgG and anti-HEV IgM were 18.67% and 1.49%, respectively. Statistical analysis showed that there was a difference in the reactivity rate of anti-HEV IgG between the HBV-infected donors and the normal donors (34.63% vs 18.67%, P<0.05), but no difference in the reactivity rate of anti-HEV IgM (1.11% vs 1.49%, P>0.05). No HEV Ag was detected in either group of blood donors. Multivariate logistic regression analysis indicated that age was an independent risk factor for anti-HEV IgG reactivity in both groups of blood donors. [Conclusion] The reactivity rate of anti-HEV IgG among HBV-infected blood donors was significantly higher than that in the normal donors in Wuhan, with age being an independent risk factor. Therefore, for HBV-infected donors, it is essential to strengthen and prioritize the prevention and treatment of HEV to reduce the spread of HEV.

Hearing loss is the most prevalent disabling disease. Cochlear implantation（CI） serves as the primary intervention for severe to profound hearing loss. This consensus systematically explores the value of genetic diagnosis in the pre-operative assessment and efficacy prognosis for CI. Drawing upon domestic and international research and clinical experience, it proposes an evidence-based medicine three-tiered prognostic classification system（Favorable, Marginal, Poor）. The consensus focuses on common hereditary non-syndromic hearing loss（such as that caused by mutations in genes like GJB2, SLC26A4, OTOF, LOXHD1） and syndromic hereditary hearing loss（such as Jervell & Lange-Nielsen syndrome and Waardenburg syndrome）, which are closely associated with congenital hearing loss, analyzing the impact of their pathological mechanisms on CI outcomes. The consensus provides recommendations based on multiple round of expert discussion and voting. It emphasizes that genetic diagnosis can optimize patient selection, predict prognosis, guide post-operative rehabilitation, offer stratified management strategies for patients with different genotypes, and advance the application of precision medicine in the field of CI.
Humans ; Cochlear Implantation ; Prognosis ; Hearing Loss/surgery* ; Consensus ; Connexin 26 ; Mutation ; Sulfate Transporters ; Connexins/genetics*

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

ObjectiveTo establish a geographical origin identification model for Foeniculi Fructus from Haiyuan， providing a new technical reference for the protection of Haiyuan's geo-authentic medicinal materials and its designation as a national geographical indication agricultural product. MethodsSamples of Foeniculi Fructus were collected from eight producing areas， including Minqin （Gansu）， Bozhou （Anhui）， Qingdao （Shandong）， Dezhou （Shandong）， Urumqi （Xinjiang）， Nujiang （Yunnan）， Gutuo （Inner Mongolia）， and Haiyuan （Ningxia）. Gas chromatography-ion mobility spectrometry （GC-IMS） was used to detect the volatile organic compounds （VOCs） in samples from these geographic origins. VOCs were qualitatively analyzed through dual matching with the National Institute of Standards and Technology （NIST） mass spectral database and the IMS drift time database. Using the Reporter module and Gallery Plot visualization tools within the LAV analytical platform， VOC fingerprint profiles characterizing geographic origins were constructed. A non-targeted analytical strategy was adopted， and 97 VOCs detected via GC-IMS were subjected to principal component analysis （PCA） and partial least squares discriminant analysis （PLS-DA） based on their differential distribution patterns to construct an origin identification model for Foeniculi Fructus from Haiyuan region. Key discriminative markers were screened using variable importance in projection （VIP） values greater than 1. ResultsA total of 97 VOCs were identified， including alcohols， aldehydes， ketones， esters， organic acids， terpenoids， ethers， alkenes， and benzenes. The PLS-DA model， based on VOCs data obtained by GC-IMS， effectively distinguished Foeniculi Fructus in Haiyuan region from those of other origins. During cross-validation， the model achieved a prediction parameter （Q²） of 0.976 and a goodness-of-fit parameter （R²） of 0.936， with no overfitting observed in permutation testing. Twelve key flavor markers with VIP > 1 were identified as characteristic indicators of Haiyuan origin. ConclusionA stable and highly predictive origin identification model for Foeniculi Fructus from Haiyuan was successfully established using GC-IMS technology， PLS-DA， and VIP-based marker screening. This model provides a novel technical strategy for accurately distinguishing Foeniculi Fructus in Haiyuan region from other regional varieties and offers new technical support for its protection as a geo-authentic medicinal material and a nationally designated geographical indication agricultural product in China.

ObjectiveTo establish a geographical origin identification model for Foeniculi Fructus from Haiyuan， providing a new technical reference for the protection of Haiyuan's geo-authentic medicinal materials and its designation as a national geographical indication agricultural product. MethodsSamples of Foeniculi Fructus were collected from eight producing areas， including Minqin （Gansu）， Bozhou （Anhui）， Qingdao （Shandong）， Dezhou （Shandong）， Urumqi （Xinjiang）， Nujiang （Yunnan）， Gutuo （Inner Mongolia）， and Haiyuan （Ningxia）. Gas chromatography-ion mobility spectrometry （GC-IMS） was used to detect the volatile organic compounds （VOCs） in samples from these geographic origins. VOCs were qualitatively analyzed through dual matching with the National Institute of Standards and Technology （NIST） mass spectral database and the IMS drift time database. Using the Reporter module and Gallery Plot visualization tools within the LAV analytical platform， VOC fingerprint profiles characterizing geographic origins were constructed. A non-targeted analytical strategy was adopted， and 97 VOCs detected via GC-IMS were subjected to principal component analysis （PCA） and partial least squares discriminant analysis （PLS-DA） based on their differential distribution patterns to construct an origin identification model for Foeniculi Fructus from Haiyuan region. Key discriminative markers were screened using variable importance in projection （VIP） values greater than 1. ResultsA total of 97 VOCs were identified， including alcohols， aldehydes， ketones， esters， organic acids， terpenoids， ethers， alkenes， and benzenes. The PLS-DA model， based on VOCs data obtained by GC-IMS， effectively distinguished Foeniculi Fructus in Haiyuan region from those of other origins. During cross-validation， the model achieved a prediction parameter （Q²） of 0.976 and a goodness-of-fit parameter （R²） of 0.936， with no overfitting observed in permutation testing. Twelve key flavor markers with VIP > 1 were identified as characteristic indicators of Haiyuan origin. ConclusionA stable and highly predictive origin identification model for Foeniculi Fructus from Haiyuan was successfully established using GC-IMS technology， PLS-DA， and VIP-based marker screening. This model provides a novel technical strategy for accurately distinguishing Foeniculi Fructus in Haiyuan region from other regional varieties and offers new technical support for its protection as a geo-authentic medicinal material and a nationally designated geographical indication agricultural product in China.

AIM: To investigate the correlation of serum interleukin-35(IL-35), immunoglobulin 4/immunoglobulin(IgG4/IgG), thyroid stimulating immunoglobulin(TSI)levels with the activity and severity of thyroid associated ophthalmopathy(TAO).METHODS:Prospective study. A total of 148 TAO patients admitted to our hospital from January 2023 to July 2024 were selected as the observation group. They were assigned into an active group(75 cases)and an inactive group(73 cases)based on their activity level, and were assigned into a severe group(95 cases)and a mild group(53 cases)based on the severity of their condition; another 148 healthy patients who underwent physical examinations were regarded as the control group. The levels of IL-35, IgG4/IgG, and TSI in serum were compared between the two groups. The correlation between serum IL-35, IgG4/IgG, and TSI levels and TAO activity and severity of illness were analyzed. A multivariate Logistic regression was performed to analyze the influencing factors of TAO patients developing severe symptoms. ROC curve was applied to analyze the diagnostic value of serum IL-35, IgG4/IgG, and TSI levels for the development of severe TAO in patients.RESULTS: Compared with the control group, the serum IL-35 level in the observation group was significantly lower, while IgG4/IgG and TSI levels were significantly higher(all P<0.01). Compared with non-active TAO patients, active TAO patients had significantly lower serum IL-35 level and significantly higher IgG4/IgG and TSI levels(all P<0.01). Compared with the mild TAO patients, severe TAO patients had significantly lower serum IL-35 level and significantly higher disease duration, IgG4/IgG, and TSI levels(all P<0.01). The serum IL-35 level was negatively correlated with TAO activity and disease severity(r=-0.529, -0.554, both P<0.01), while serum IgG4/IgG level was positively correlated with TAO activity and disease severity(r=0.625, 0.663, both P<0.01). Serum TSI levels were positively correlated with TAO activity and disease severity(r=0.594, 0.607, both P<0.01). Multivariate Logistic regression analysis showed that serum IL-35, IgG4/IgG, and TSI levels were all factors influencing the progression of TAO patients to severe disease(all P<0.01). The areas under the curve(AUC)for diagnosing the progression of TAO patients to severe disease using serum IL-35, IgG4/IgG, and TSI levels were 0.868, 0.859, and 0.830, respectively. The combined AUC for the three markers was 0.955, significantly higher than that of each individual marker(Zthree factors combination-IL-35=2.893, Zthree factors combination-IL-35=3.510, Zthree factors combination-IL-35=4.157, P=0.004, <0.01, <0.01).CONCLUSION: Serum IL-35 level is significantly downregulated in TAO patients, while IgG4/IgG and TSI levels are significantly upregulated. The levels of IL-35, IgG4/IgG, and TSI are correlated with the activity and severity of TAO, and their combination has high diagnostic value for TAO developing into severe.