ObjectiveTo elucidate the chemical composition of the reference sample of Jinshui Liujunjian and its distribution characteristics in blood and tissues of rats. MethodsUltra performance liquid chromatography-quadrupole-time-of-flight tandem mass spectrometry（UPLC-Q-TOF-MS/MS） was used to detect the reference sample solution， plasma， and tissue samples of Jinshui Liujunjian under positive and negative ion modes， respectively. Qualitative Analysis 10.0 software and a self-constructed database were employed for primary mass spectrum matching.Compound identification was further validated by comparing retention times， secondary mass spectral fragments， reference standards， and literature data to deduce fragmentation pathways. ResultsA total of 122 compounds were identified in the reference sample of Jinshui Liujunjian， including 47 flavonoids， 5 amino acids， 13 iridoids， 16 triterpenoid saponins， etc.， of which 42 compounds were confirmed by comparison with reference substances. A total of 21 prototype components were identified in blood components； 50 prototype components were identified in different tissues， among which 13， 10， 7， 21， 11， 6， 14， and 40 prototype components were identified in the heart， liver， spleen， lung， kidney， brain， large intestine， and stomach， respectively. Among them， 7 compounds such as ferulic acid， glycyrrhizic acid， and nobiletin were exposed in the target organs of lung and kidney. ConclusionThis study elucidates the material basis of the reference samples of Jinshui Liujunjian， primarily composed of flavonoids and triterpenoid saponins， along with their in vivo distribution characteristics. These findings provide a scientific basis for establishing quality evaluation indicators and offer references for subsequent pharmacodynamic and pharmacokinetic investigations.

ObjectiveTo elucidate the chemical composition of the reference sample of Jinshui Liujunjian and its distribution characteristics in blood and tissues of rats. MethodsUltra performance liquid chromatography-quadrupole-time-of-flight tandem mass spectrometry（UPLC-Q-TOF-MS/MS） was used to detect the reference sample solution， plasma， and tissue samples of Jinshui Liujunjian under positive and negative ion modes， respectively. Qualitative Analysis 10.0 software and a self-constructed database were employed for primary mass spectrum matching.Compound identification was further validated by comparing retention times， secondary mass spectral fragments， reference standards， and literature data to deduce fragmentation pathways. ResultsA total of 122 compounds were identified in the reference sample of Jinshui Liujunjian， including 47 flavonoids， 5 amino acids， 13 iridoids， 16 triterpenoid saponins， etc.， of which 42 compounds were confirmed by comparison with reference substances. A total of 21 prototype components were identified in blood components； 50 prototype components were identified in different tissues， among which 13， 10， 7， 21， 11， 6， 14， and 40 prototype components were identified in the heart， liver， spleen， lung， kidney， brain， large intestine， and stomach， respectively. Among them， 7 compounds such as ferulic acid， glycyrrhizic acid， and nobiletin were exposed in the target organs of lung and kidney. ConclusionThis study elucidates the material basis of the reference samples of Jinshui Liujunjian， primarily composed of flavonoids and triterpenoid saponins， along with their in vivo distribution characteristics. These findings provide a scientific basis for establishing quality evaluation indicators and offer references for subsequent pharmacodynamic and pharmacokinetic investigations.

ObjectiveThe accuracy of Y-chromosome haplogroup assignment is crucial for tracing paternal lineage in male samples. With the advancement of high-throughput sequencing technologies, high-density Y-SNP genotyping from whole-genome or array-based data has become a standard method for determiningY-chromosome haplogroups. This study systematically evaluated the performance of 4 commonly used high-density SNP genotyping systems—namely, the Global Screening Array (GSA), Chinese Genotyping Array (CGA), Affymetrix array, and the 1240K capture panel—for haplogroup assignment. This work provides a reference for data comparison across different systems. MethodsWe extracted genotype data for the 4 Y-SNP panels from 30× whole-genome sequencing (WGS) data of 1 590 male samples from the 1000 Genomes Project. Additionally, GSA array genotype data from 384 relative pairs (spanning 1st- to 12th-degree relationships) from 109 Chinese Han families were collected. Haplogroup assignment was performed using Y-LineageTracker v1.3.0 software. We assessed the concordance and resolution of haplogroup assignments between the four Y-SNP panels and the WGS data. The consistency and resolution of haplogroup assignments were also evaluated for both the 1000 Genomes Project samples and the 109 family samples collected in this study. Furthermore, the impact of varying numbers of Y-SNPs on haplogroup assignment was examined. ResultsThe GSA and CGA panels demonstrated superior resolution and discrimination of haplogroup subclades compared with the other two panels. The haplogroup assignments from the GSA, CGA, and 1240K panels showed high concordance with WGS data, with consistency rates exceeding 88.70%, whereas the Affymetrix platform exhibited a significantly lower consistency rate of 61.89%. Specifically, the GSA and CGA panels consistently demonstrated superior performance compared with the other two panels in the assignment of haplogroups O-M175 and H-L901, achieving complete concordance (100%) for both haplogroups. In contrast, the Affymetrix panel erroneously assigned all individuals belonging to haplogroup O-M175 to haplogroup K2-M526. Furthermore, its accuracy for haplogroup H-L901 was exceedingly low, at merely 1.41%. This poor performance was characterized by the misassignment of 98.59% of H-L901 samples—specifically, 1.41% to J-M304 and a predominant 97.18% to F-M89. For haplogroup R-M207, all four panels exhibited uniformly high levels of consistency, with concordance values exceeding 94.00%. Notably, for haplogroup E-M96, the 1240K and Affymetrix panels outperformed the GSA and CGA panels in terms of concordance, representing the first instance in which these two panels surpassed the latter. Conversely, for haplogroups J-M304, Q-M242, and I-M170, all 4 panels showed relatively elevated misclassification rates, with the Affymetrix array demonstrating the poorest overall performance. None of the four panels showed any discordant haplogroup assignments among the familial relative pairs analyzed. A positive correlation was observed between the number of Y-SNPs (ranging from 1 000 to 10 000) and classification consistency; however, classification consistency plateaued when the number of Y-SNPs exceeded 10 000. Furthermore, a random sampling analysis conducted on the GSA and CGA panels demonstrated that the haplogroup misclassification rate exhibited negligible fluctuation across the Y-SNP range of 500 to 1 000. Conversely, a marked enhancement in classification consistency was observed as the number of markers increased from 1 000 to 5 000, ultimately reaching a plateau within the interval of 5 000 to 8 000 markers. ConclusionThese findings indicate that the GSA and CGA panels provide high resolution and concordance, delivering reliable Y-haplogroup assignment for forensic investigations.

ObjectiveThe accuracy of Y-chromosome haplogroup assignment is crucial for tracing paternal lineage in male samples. With the advancement of high-throughput sequencing technologies, high-density Y-SNP genotyping from whole-genome or array-based data has become a standard method for determiningY-chromosome haplogroups. This study systematically evaluated the performance of 4 commonly used high-density SNP genotyping systems—namely, the Global Screening Array (GSA), Chinese Genotyping Array (CGA), Affymetrix array, and the 1240K capture panel—for haplogroup assignment. This work provides a reference for data comparison across different systems. MethodsWe extracted genotype data for the 4 Y-SNP panels from 30× whole-genome sequencing (WGS) data of 1 590 male samples from the 1000 Genomes Project. Additionally, GSA array genotype data from 384 relative pairs (spanning 1st- to 12th-degree relationships) from 109 Chinese Han families were collected. Haplogroup assignment was performed using Y-LineageTracker v1.3.0 software. We assessed the concordance and resolution of haplogroup assignments between the four Y-SNP panels and the WGS data. The consistency and resolution of haplogroup assignments were also evaluated for both the 1000 Genomes Project samples and the 109 family samples collected in this study. Furthermore, the impact of varying numbers of Y-SNPs on haplogroup assignment was examined. ResultsThe GSA and CGA panels demonstrated superior resolution and discrimination of haplogroup subclades compared with the other two panels. The haplogroup assignments from the GSA, CGA, and 1240K panels showed high concordance with WGS data, with consistency rates exceeding 88.70%, whereas the Affymetrix platform exhibited a significantly lower consistency rate of 61.89%. Specifically, the GSA and CGA panels consistently demonstrated superior performance compared with the other two panels in the assignment of haplogroups O-M175 and H-L901, achieving complete concordance (100%) for both haplogroups. In contrast, the Affymetrix panel erroneously assigned all individuals belonging to haplogroup O-M175 to haplogroup K2-M526. Furthermore, its accuracy for haplogroup H-L901 was exceedingly low, at merely 1.41%. This poor performance was characterized by the misassignment of 98.59% of H-L901 samples—specifically, 1.41% to J-M304 and a predominant 97.18% to F-M89. For haplogroup R-M207, all four panels exhibited uniformly high levels of consistency, with concordance values exceeding 94.00%. Notably, for haplogroup E-M96, the 1240K and Affymetrix panels outperformed the GSA and CGA panels in terms of concordance, representing the first instance in which these two panels surpassed the latter. Conversely, for haplogroups J-M304, Q-M242, and I-M170, all 4 panels showed relatively elevated misclassification rates, with the Affymetrix array demonstrating the poorest overall performance. None of the four panels showed any discordant haplogroup assignments among the familial relative pairs analyzed. A positive correlation was observed between the number of Y-SNPs (ranging from 1 000 to 10 000) and classification consistency; however, classification consistency plateaued when the number of Y-SNPs exceeded 10 000. Furthermore, a random sampling analysis conducted on the GSA and CGA panels demonstrated that the haplogroup misclassification rate exhibited negligible fluctuation across the Y-SNP range of 500 to 1 000. Conversely, a marked enhancement in classification consistency was observed as the number of markers increased from 1 000 to 5 000, ultimately reaching a plateau within the interval of 5 000 to 8 000 markers. ConclusionThese findings indicate that the GSA and CGA panels provide high resolution and concordance, delivering reliable Y-haplogroup assignment for forensic investigations.

Background Long working hours, as a common risk factor for occupational stress, is closely related to the occurrence of depressive symptoms. Understanding how long working hours affect occupational stress and depressive symptoms will inform occupational health interventions. Objective To quantify the impact of long working hours exposure on occupational stress and depressive symptoms among Internet industry employees, translate black-box outputs into actionable insights, and demonstrate the value of interpretable machine learning for early-warning occupational-health surveillance. Methods A dataset was derived from a cross-sectional survey involving 2866 internet industry employees in China. This survey was part of the project Risk Assessment Of Long Working Hour Exposure And Its Adverse Health Effects, conducted by the National Institute for Occupational Health and Poisoning Control, Chinese Center for Disease Control and Prevention, from 2021 to 2023. Working hours, occupational stress and depressive symptoms were quantified with a set of structured questionnaires including the Core Occupational Stress Scale and the Patient Health Questionnaire. Pairwise associations were screened by Mantel tests and variance-inflation factors. Key predictors identified through feature selection were fed into six machine-learning risk-prediction models. Visual interpretation was provided by feature importance, Shapley additive explanations (SHAP) and local interpretable model-agnostic explanations (LIME), while directed causal effects and intervention impacts of prolonged working hours exposure on occupational stress and depressive symptoms were dissected with causal explanation of features techniques. Results The positive rates of occupational stress and depressive symptoms among internet employees were 12.9% and 77.8% respectively. Twelve core features for occupational stress and nine for depressive symptoms were retained after selection. After these features were supplied to six predictive algorithms and evaluated on five metrics, the Light Gradient Boosting Machine (LGBM) achieved the highest accuracy—0.89 for occupational stress and 0.79 for depressive symptoms on the hold-out test set. The feature-importance rankings converged on fatigue accumulation and life satisfaction as dominant drivers for both outcomes, whereas weekly working hours and daily overtime emerged as the principal exposure-related predictors. The SHAP summary plots revealed that longer weekly hours and daily overtime systematically elevated the probability of occupational stress. The causal feature explanation further quantified that ascending one category in weekly working hours increased the probability of occupational stress by 7.04%. Conclusion Exposure to long working hours is associated with both occupational stress and depressive symptoms among internet industry employees. Interpretable machine-learning frameworks translate these associations into transparent, defensible drivers, enabling precise identification of the pivotal factors and their interplay. This evidence base equips occupational-health practitioners with actionable insights for designing targeted prevention and intervention strategies.

AIM:To develop a SAS macro tool for calculating pharmacokinetic(PK)parameters using non-compartmental analysis(NCA)that is applica-ble to various routes and frequencies of drug ad-ministration.METHODS:By incorporating the prin-ciples of NCA PK parameter calculations,the struc-ture of the SAS macro program is designed and pro-grammed.The usage of the macro program is dem-onstrated through simulated examples.RESULTS:Based on NCA,this study programmed a SAS macro using SAS macro language to calculate both single-dose and steady-state complete PK parameters.It supports three modes of drug administration:ex-travascular,intravenous infusion,and intravenous bolus injection,as well as four algorithms for calcu-lating area under the curve(AUC).By comparing with authoritative calculation software,the results were found to be highly consistent across various calculation scenarios,with more accurate calcula-tions for the time of first measurable concentration(Tlag).CONCLUSION:The％M_CREATE_PK_PARAM-ETERS macro is a SAS macro that supports compre-hensive PK parameter calculations in NCA,provid-ing an accurate,efficient,and flexible solution for PK parameter estimation with broad application prospects.

OBJECTIVE To investigate the pathogens of pulmonary infection and their influencing factors in patients with initially treated pulmonary tuberculosis(PTB),and to analyze the diagnostic efficacy of systemic inflamma-tory response index(SIRI)for pulmonary infection.METHODS The clinical data of 459 patients with initially trea-ted PTB admitted to Changsha Central Hospital from Jan.to Dec.2023 were retrospectively analyzed,and they were divided into the combined group(n=90)and the non-combined group(n=369)according to the combina-tion of other pulmonary infections at the time of consultation.The pathogens of pulmonary infection,as well as the levels of monocytes(MN),neutrophils(NE),and lymphocytes(LYM),were analyzed,and the SIRI was calculated.The risk factors of pulmonary infection in initially treated PTB patients were explored by multivariate logistic regression analysis,and the diagnostic value of SIRI for pulmonary infection in initially treated PTB pa-tients was evaluated by receiver operating characteristic(ROC)curve.RESULTS The rate of pulmonary infection in 459 initially treated PTB patients was 19.61％(90/459),and a total of 103 pathogens were detected,of which 64 strains of gram-negative bacteria accounted for 62.14％,followed by 28 strains of gram-positive bacteria ac-counting for 27.18％,and 11 strains of fungi accounting for 10.68％,with Pseudomonas aeruginosa,Klebsiella pneumoniae and Staphylococcus aureus being predominant.Age(OR=1.908,95％CI:1.375-2.647),comor-bid diabetes(OR=2.073,95％CI:1.462-2.938),pulmonary cavities(OR=2.323,95％CI:1.588-3.398),delayed medical treatment(OR=2.024,95％CI:1.467-2.791)and elevated SIRI(OR=2.855,95％CI:1.851-4.402)were risk factors for pulmonary infection in initially treated PTB patients(P＜0.05).ROC curve analysis showed that the area under the curve(AUC)of SIRI for diagnosing pulmonary infection in initially treated PTB cases was 0.876(0.824-0.928),with a cut-off value of 4.96×109/L,a specificity of 85.91％,and a sen-sitivity of 86.67％.CONCLUSIONS The main pathogens in patients with initially treated PTB complicated with pulmonary infection are Pseudomonas aeruginosa and Klebsiella pneumoniae.Age,comorbid diabetes mellitus,pulmonary cavity,delayed medical treatment and elevated SIRI are all influencing factors for pulmonary infection in patients with initially treated PTB,and monitoring of SIRI can to some extent assist in diagnosis of whether pa-tients with initial treatment PTB are complicated with other pulmonary infections.

Objective:To explore the interrelationship between brain functional networks and features in functional magnetic resonance imaging(fMRI)of patients with Alzheimer's disease(AD),and to construct mixed-function networks(MFN),and apply them in machine learning classification models to improve the accuracy of AD classification.Methods:102 AD patients and 227 healthy subjects in the Alzheimer's Neuroimaging Initiative(ADNI)dataset were retrospectively analyzed.The partial correlation brain network of the blood oxygen level dependent(BOLD)signal was calculated and fused with low-frequency wave amplitude(ALFF),fractional low-frequency wave amplitude(fALFF)and local consistency(ReHo)features to construct MFN.Network topology parameters were extracted,and a variety of machine learning classification models were constructed based on MFN topological parameters,accuracy,precision,recall and area under the curve(AUC)were used to evaluate the predictive efficiency of the models.Results:By constructed MFN and calculated intra group to inter group ratio(IIGR),35 features could be obtained from ALFF,fALFF and ReHo feature topological parameter analysis,after rank sum test and FDR correction,there were statistical differences among 28 features(P＜0.05).The classification results show that,all the five classifiers have high classification performance on the test data set.The accuracy,precision and recall rates of random forest(RF),adaptive lifting algorithm(AdaBoost),guided aggregation algorithm(Bagging)and support vector machine(SVM)were all 99.7％,and the AUC values were up to 100％,99.5％,99.1％and 99.5％,respectively.The accuracy(98.5％),precision(98.5％),recall(98.5％),and AUC(99.1％)of the multi-layer perceptron(MLP)were slightly lower than other models,but remained excellent.It was worth noting that RF has the highest AUC value of all models at 100.0％,while Bagging has the lowest AUC value(99.1％)in the integrated approach.The results of performance comparison show that,MFN classification model can significantly improve the recognition and classification of AD disease,and greatly improve the performance of various indicators of the classifier.The results showed that,MFN classification model was superior to intelligent classification based fusion,DBN-based multitask learning,PVT-TSVM,unsupervised learning and clustering,SVM and SVM of degree 3 polynomial kernel function in key indicators such as accuracy(99.13％),AUC(99.42％),recall rate(99.46％)and specificity(99.42％)with plasma proteins,machine learning algorithms.It was further proved that MFN classification model has good generalization ability and robustness in AD disease classification.Conclusion:The AD classification model constructed based on brain mixed function network topology parameters and machine learning can improve the accuracy of AD classification.

Objective:To study the role of cadherin 1 (CDH1) gene DNA methylation in autoimmune thyroiditis (AIT) in population from different water-iodine regions.Methods:From May to June 2019, the information of AIT cases and healthy individuals in Shandong Province were collected in three types of water-iodine regions: iodine-fortification (IF) region, iodine-adequate (IA) region and iodine-excess (IE) region. A case-control study design was applied to match 176 AIT cases (case group) with age, gender, body mass index, and place of residence in a 1 ∶ 1 ratio to 176 healthy individuals (control group). Fasting urine and whole blood samples were collected to test the contents of urinary iodine, thyroid function indicators [serum free triiodothyronine (FT 3), free thyroxine (FT 4), thyroid stimulating hormone (TSH)], and serum iodine. The DNA methylation levels of the target region of the CDH1 gene and its four CpG sites in whole blood were determined using methylation sequencing technology for target regions (MethylTarget TM). Results:The DNA methylation level of the target region of CDH1 gene in the case group was 0.832 ± 0.044, and that in the control group was 0.828 ± 0.049, there was no statistically significant difference between the two groups ( t = 0.76, P = 0.448). There was no statistically significant difference in DNA methylation levels of the four CpG sites in the target region of CDH1 gene between the case group and the control group ( P > 0.05). There was no statistically significant difference in the DNA methylation level of the CDH1 gene target region between the case group and the control group in IF, IA and IE regions ( P > 0.05). The detection results of DNA methylation levels at CpG sites in the target region of CDH1 gene in different water iodine regions showed that the DNA methylation level at site 83 in case group in IF region was higher than that in the control group ( t = 2.30, P = 0.023). However, there was no statistically significant difference in the DNA methylation levels of the four CpG sites between the case group and the control group in IA and IE regions ( P > 0.05). The DNA methylation level of CDH1 gene target region in AIT patients was not significantly correlated with urinary iodine, serum iodine, and serum FT 3, FT 4, and TSH contents ( P > 0.05), but was significantly negatively correlated with age ( r =-0.19, P = 0.014). Conclusions:The DNA methylation level at CpG site 83 of CDH1 gene in AIT patients in IF region is significantly higher than that in control population, indicating that DNA methylation at this locus may be involved in the occurrence and development of AIT after iodine fortification. The DNA methylation level of CDH1 gene is negatively correlated with age.

Objective:To analyze the iodine nutrition status of children and adolescents and influencing factors in Zhejiang Province, providing scientific basis for optimizing iodine deficiency disorders (IDD) prevention and control strategies.Methods:In June 2022, a multistage stratified sampling method was used to divide 16 counties (cities, districts, abbreviated as counties) in Zhejiang Province into three categories based on their geographical locations (average distance from the coastline): coastal areas (including Dinghai District, Jiaojiang District, Sanmen County, Cixi City and Lucheng District), sub-coastal areas (including Wuxing District, Haining City, Linping District, Fuyang District and Fenghua District), and inland areas(including Suichang County, Changshan County, Shengzhou City, Jindong District, Dongyang City and Yongjia County). One county was selected from each category, and one township (street) was selected from each county. Two administrative villages (neighborhood committees) were selected from each township (street). Ten households including all children and adolescents aged 6-17 in each household were selected from each administrative village (neighborhood committee). Demographic information and personal dietary characteristics were collected via questionnaires, while household salt and a random urine sample were tested for iodine level. Trend analysis was conducted using a χ 2trend test, and a multivariate logistic stepwise regression model was used to analyze the influencing factors of urinary iodine levels. Results:A total of 755 children and adolescents aged 6-17 were selected, including 387 males (51.26%) and 368 females (48.74%), with an age of (11.24 ± 3.32) years. There were 269 children and adolescents in coastal areas (35.63%) and 409 children and adolescents in urban areas (54.17%). A total of 755 household salt samples were collected, with a median salt iodine concentration of 21.80 mg/kg. These included 263 non-iodized salt samples, 38 unqualified iodized salt samples, and 454 qualified iodized salt samples. The coverage rate of iodized salt was 65.17% (492/755), and the consumption rate of qualified iodized salt was 60.13% (454/755). The distribution of salt iodine quality among children and adolescents in different geographical locations showed statistically significant differences (χ 2 = 111.95, P < 0.001), with the proportion of non-iodized salt gradually decreasing from coastal areas to inland areas (χ 2trend = 90.17, P < 0.001). A total of 755 urine samples were collected, with a median urinary iodine concentration of 186.60 μg/L. The proportions of urinary iodine < 100, 100-199, 200-299, and ≥300 μg/L were 16.95% (128/755), 37.62% (284/755), 24.37% (184/755), and 21.06% (159/755), respectively. The χ 2trend test revealed a nonlinear positive correlation between salt iodine level and urinary iodine level (χ 2regression = 21.98, P < 0.001; χ 2partial = 6.96, P < 0.001). The frequency distribution of urinary iodine in children and adolescents from different geographical locations and between urban and rural areas showed statistically significant differences (χ 2 = 29.63, 16.56, P < 0.001). Among them, the proportion of children and adolescents with urinary iodine < 100 μg/L gradually decreasing from coastal areas to inland areas (χ 2trend = 6.15, P = 0.013). The results of multivariate logistic regression analysis revealed that sub-coastal regions, inland regions, and urban-rural regions ( OR = 1.57, 1.53, 1.64, 95% CI: 1.11-2.24, 1.03-2.27, 1.17-2.32, P < 0.05) were significantly associated with urinary iodine levels in children and adolescents aged 6-17. Conclusions:In 2022, the iodine nutrition of children and adolescents in Zhejiang Province is generally suitable, but there is a risk of iodine deficiency among coastal children and adolescents. Geographic location and urban/rural areas are influencing factors on iodine nutrition status of children and adolescents in Zhejiang Province.