BACKGROUND:Studies have shown that rheumatoid arthritis and osteoporosis are positively correlated,but the causal relationship and related mechanisms have not yet been confirmed.With the cross-fertilization of computer science and life sciences,Mendelian randomization and bioinformatics analyses based on genome-wide association study(GWAS)and transcriptome sequencing data can assess the causal relationship between two diseases,explore the related mechanisms,and mine the therapeutic targets,which will be beneficial to the precision treatment of rheumatoid arthritis combined with osteoporosis.OBJECTIVE:To explore the causal relationship between rheumatoid arthritis and osteoporosis using two-sample Mendelian randomization and to mine potential co-morbid targets and potential targeted drugs through summary-data-based Mendelian randomization and bioinformatics analyses,aiming to provide theoretical basis for mechanism exploration and precision treatment in the field of rheumatoid arthritis combined with osteoporosis.METHODS:(1)Firstly,GWAS data of rheumatoid arthritis,osteoporosis,and cis-expression quantitative trait locus(cis-eQTL)in Asian and European populations were downloaded from the GWAS Catalog,IEU Open GWAS,FinnGen,and eQTLGen databases,and were used for two-sample Mendelian randomization analysis and summary-data-based Mendelian randomization analysis.(2)Transcriptome sequencing data of rheumatoid arthritis(GSE93272 and GSE15573)were downloaded from the GEO database for bioinformatics analysis.(3)Subsequently,forward and inverse Mendelian randomization analyses between rheumatoid arthritis and osteoporosis were performed,and inverse variance weighted was used as the main metric for the analyses,and the results were corroborated with MR Egger,simple mode,weighted median and weighted mode.(4)Then,the genes closely related to rheumatoid arthritis and osteoporosis were identified based on the summary-data-based Mendelian randomization analysis,and the co-disease targets of rheumatoid arthritis and osteoporosis were mined based on cross-analysis.Meanwhile,the biological functions of the co-morbid targets were verified based on bioinformatics analysis and cellular experiments.(5)In addition,a rheumatoid arthritis risk prediction nomogram was constructed based on DYRK2,and its prediction performance was verified by receiver operating characteristic curve,correction curve and decision curve.Finally,the target potential drugs were mined based on Enrichr database and molecular docking was performed.RESULTS AND CONCLUSION:(1)Forward Mendelian randomization analysis of rheumatoid arthritis and osteoporosis showed statistically significant results except for GCST90044540 and GCST90086118,and all other results indicated a significant causal relationship and positive correlation between rheumatoid arthritis and osteoporosis.(2)Inverse Mendelian randomization analysis suggested that no significant causal relationship was seen between osteoporosis and rheumatoid arthritis.(3)Summary-data-based Mendelian randomization analysis identified a total of 412 and 344 genes positively associated with rheumatoid arthritis and osteoporosis,and 421 and 347 genes negatively associated.Based on the cross-analysis,26 co-morbid genes were subsequently obtained.Among them,DYRK2 was a potential therapeutic target,and subsequent bioinformatics analysis and cellular experiments confirmed its important role in the progression of rheumatoid arthritis and osteoporosis.(4)Furthermore,the constructed nomogram has excellent predictive performance.Finally,four potential DYRK2-targeting drugs(undecanoic acid,metyrapone,JNJ-38877605,and ACA)were discovered and molecular docking also demonstrated reliable targeting ability.(5)In conclusion,based on GWAS data from Asian and European populations,we successfully demonstrated that rheumatoid arthritis and osteoporosis are causally related at the genetic level,DYRK2 is a potential therapeutic target,and four small molecules are potential target drugs.

OBJECTIVE:Based on different algorithms of machine learning,the prediction model of lumbar disc herniation has become a trend and hot spot in the development of precision medicine.However,there is limited evidence on the reporting quality and methodological quality of prediction models of lumbar disc herniation outcomes using machine learning.This article is aimed to explore the performance of machine learning algorithms in predicting the prognosis of lumbar disc herniation by comprehensively analyzing the report quality and risk of bias of previous studies that developed and validated prognosis prediction models based on machine learning through a comprehensive literature search,in order to explore the performance of machine learning algorithms in predicting the prognosis of lumbar disc herniation.METHODS:The databases of CNKI,WanFang,VIP,SinOMED,PubMed,Web of Science,Embase,and The Cochrane Library were searched by computer.Studies on the use of machine learning to develop(and/or validate)prognostic prediction models for lumbar disc herniation were collected from the inception of the database to December 31,2023.Two researchers independently screened the literature,extracted data,and assessed the risk of bias of the included studies.The reporting quality and risk of bias of the included studies were assessed by the Multivariable Transparent Reporting of Predictive Models(TRIPOD)statement and the Predictive Model Risk of Bias Assessment Tool(PROBAST).The results of the evaluation were analyzed using descriptive statistics and visual charts.RESULTS:(1)A total of 23 articles were included,and the TRIPOD compliance of each study ranged from 11％to 87％,with a median compliance of 54％.The quality of reporting of titles,detailed descriptions of treatment measures,blinding of predictors,handling of missing data,details of risk stratification,specific procedures for enrollment,model interpretation,and model performance was mostly poor,with TRIPOD adherence rates ranging from 4％to 35％.(2)Of all included studies,61％had a high risk of bias and 39％had an unclear overall risk of bias.The area under the curve,accuracy,sensitivity and specificity were used to evaluate the performance of the model.The areas under the curve of 20 models were reported,ranging from 0.561 to 0.999.Three models reported the accuracy of the model,ranging from 82.07％to 89.65％.(3)Among all included studies,the statistical analysis domain was most often assessed as having a high risk of bias,mainly due to the small number of valid samples,the selection of predictors based on univariate analysis and the lack of calibration and discrimination assessment of the model in the study.CONCLUSION:These results indicate that machine learning can achieve good predictive ability in the development and validation of prognostic models for lumbar disc herniation.The commonly used algorithms include regression algorithm,support vector machine,decision tree,random forest,artificial neural network,naive Bayes and other algorithms.Reasonable algorithms combined with clinical practice can improve the accuracy of prognosis prediction of lumbar disc herniation.However,the reporting and methodological quality of prognosis prediction models based on machine learning are poor,the prediction performance of different models varies greatly,and the generalization and extrapolation of research models are unclear.There is an urgent need to improve the design,implementation and reporting of such studies.To promote the application of machine learning in the clinical practice of lumbar disc herniation prediction models,it is necessary to comprehensively consider various predictors related to the prognosis of the disease before modeling,and strictly follow the relevant standards of PROBAST tool during modeling.

Objective To identify factors affecting and key environmental factors of the Oncomelania hupensis snail density in the Yangtze River Delta region using machine learning methods. Methods Administrative village-level O. hupensis snail survey data in the Yangtze River Delta (including Shanghai Municipality, Jiangsu Province, Zhejiang Province and Anhui Province) from 2011 to 2021 were retrieved from the Information Management System for Parasitic Disease Control of Chinese Center for Disease Control and Prevention. Environmental factor data were captured from the Google Earth Engine platform, including elevation, slope, terrain, normalized difference vegetation index (NDVI), vegetation type, soil type, total petroleum hydrocarbon (TPH), ammonium nitrogen, inorganic nitrogen, dissolved oxygen, pH of water, chemical oxygen demand (COD) and inorganic phosphorus, and climatic factor data in the study region were retrieved from the Copernicus Climate Data Store, including annual precipitation, aridity index and annual mean temperature (AMT). O. hupensis snail survey data in the Yangtze River Delta region from 2011 to 2021 were randomly divided into a training set (70%) and a test set (30%), and five machine learning models were selected for machine learning model construction and comparative analysis of the O. hupensis snail density using the software R 4.3.0, including random forest (RF), eXtreme gradient boosting (XGBoost), support vector machine (SVM), gradient boosting machine (GBM) and neural network (NN). The XGBoost model was employed to construct a predictive model for the O. hupensis snail density, and the impact of each environmental factor on O. hupensis snail distribution was quantified. The SHapley Additive exPlanations (SHAPs) values were calculated to estimate the average contribution of each variable to the model prediction, and the core environmental factors affecting the O. hupensis snail population density were screened. Results Among the five machine learning models, the XGBoost model exhibited the optimal comprehensive performance, with the coefficient of determination (R²) of 0.855, mean squared error (MSE) of 0.188, root mean squared error (RMSE) of 0.434 and mean absolute error (MAE) of 0.155, respectively. Analysis of factors affecting the O. hupensis snail density with the XGBoost model showed that among the 16 environmental factors, the top four high-impact factors ranked by SHAPs values included annual precipitation, elevation, aridity index and NDVI, with cumulative SHAPs contributions of 75%, which was higher than that of other environmental factors. If NDVI was higher than 0.6, the O. hupensis snail density increased with NDVI and peaked if NDVI was 0.8 (1.60 snails/0.1 m²). The O. hupensis snail density increased with elevation if the elevation ranged from 14 to 40 m, and slowly rose if the annual precipitation ranged from 900 to 1 300 mm, and then increased rapidly to the peak (1.52 snails/0.1 m²) if the annual precipitation ranged from 1 300 to 1 500 mm. In addition, the O. hupensis snail density increased rapidly to the maximum (1.60 snails/0.1 m²) if the aridity index ranged from 0.8 to 1.1, and decreased gradually if the aridity index exceeded 1.1. Conclusions The XGBoost model shows excellent performance in prediction of the O. hupensis snail density and identification of key environmental factors in the Yangtze River Delta region. Annual precipitation, elevation, aridity index and NDVI are key environmental factors affecting the distribution and density of O. hupensis snails in the Yangtze River Delta region.

Objective To examine the impact of extreme temperature and drought stress on the survival of Bulinus globosus, so as to provide the theoretical evidence for the genomic research of Bulinus in absence of reference genes. Methods B. globosus snail samples were collected from Kiwani Shehia in Pemba Island, Zanzibar, Tanzania, and offspring snails were obtained through laboratory breeding and reproduction. A total of 120 10-week-old B. globosus snails from the same generation were selected and randomly assigned into four groups, including the high-temperature drought (HD) group, normal temperature drought (D) group, low-temperature drought (LD) group, and the control (C) group, of 30 snails in each group. Snails in HD, D, and LD groups were placed in beakers containing dry soil at the bottom and subsequently housed in climate chambers at 35, 26 ℃, and 10 ℃, respectively, while snails in Group C were maintained in 500 mL petri dishes containing dechlorinated tap water at 26 ℃. Following 3 days of breeding, living snails in each group were collected, and soft tissues were dissected and isolated. Total RNA was extracted from snail soft tissues for library construction, followed by high-throughput sequencing on the Illumina HiSeq 4000 sequencing system. De novo transcriptome assembly was performed using the Trinity software, and the longest transcripts were selected as unigenes. Gene functional annotations of unigenes were conducted using the Diamond software against Gene Ontology (GO) knowledgebase, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database, NCBI non-redundant (NR) protein sequences database, Protein Family (Pfam) database, and UniProtKB/Swiss-Prot (Swiss-Prot) knowledgebase. GO and KEGG enrichment analyses of differentially expressed genes (DEGs) were performed using the topGO and clusterProfiler software, respectively. In addition, four relevant genes were selected for validation using a real-time quantitative PCR (qRT-PCR) assay to verify the reliability of transcriptome sequencing results. Results Following 3 days of breeding, there were 7, 20, 28, and 30 survival B. globosus snails in HD, LD, D, and C groups, with corresponding survival rates of 23.33% (7/30), 66.67% (20/30), 93.33% (28/30), and 100.00% (30/30), respectively (χ² = 52.72, P < 0.001). De novo transcriptome assembly generated 176 942 unigenes, with annotation rates of 0.98%, 13.49%, 26.46%, 12.48%, and 14.39% against GO knowledgebase, KEGG pathway database, NR protein sequences database, Pfam database, and Swiss-Prot knowledgebase, respectively. There were 33 up-regulated and 72 down-regulated genes in Group D, 483 up-regulated and 815 down-regulated genes in Group HD, and 245 up-regulated and 172 down-regulated genes in Group LD relative to in Group C. Following removal of overlapping genes across groups and unmatched genes, 11 candidate genes were identified. GO and KEGG analyses revealed 3 heat shock protein (HSP)-related DEGs in these 11 candidate genes, which were annotated as HSP12.2, HSP70, and HSP20 genes and were all significantly up-regulated in each treatment group. Three immune and nervous system-related DEGs were identified, and were all significantly down-regulated in each treatment group, which were involved in the neural cell adhesion molecule L1-like protein pathway, fibrinogen binding protein pathway, and leukocyte elastase inhibitor-like protein pathway. qRT-PCR assay quantified that the expression trends of four genes related to temperature and drought stress across different treatment groups were highly consistent with transcriptome sequencing data. Conclusion The survival rate of B. globosus significantly reduces under combined stresses of extreme temperature and drought, possibly due to an imbalance in its cellular homeostasis regulatory system.

Cancer remains a leading cause of global mortality, necessitating the development of advanced therapeutic strategies with enhanced efficacy and reduced systemic toxicity. Among promising bioactive agents, lactoferrin (LF)—a multifunctional iron-binding glycoprotein abundantly found in mammalian milk and exocrine secretions—has garnered significant interest for its potent and multifaceted anti-cancer properties. This review provides a comprehensive analysis of the current understanding of LF’s role in oncology, encompassing its structural biology, diverse mechanisms of action, and groundbreaking advancements in its application through nano-engineering. LF exerts anti-tumor effects through multiple pathways, including extracellular action, intracellular action, and immune regulation. It demonstrates a remarkable affinity for cancer cell membranes, binding to overexpressed anionic components such as glycosaminoglycans and sialic acids, as well as to specific receptors including the low-density lipoprotein receptor-related protein-1 (LRP-1). This selective binding facilitates targeted uptake. Upon internalization, LF orchestrates a direct assault by inducing cell-cycle arrest in phases such as G0/G1 or S phase through the modulation of key regulators including cyclins, CDKs, and p53. Furthermore, it promotes programmed cell death via apoptotic pathways, involving caspase activation and downregulation of anti-apoptotic proteins such as survivin. A more recently elucidated mechanism is the induction of ferroptosis, an iron-dependent form of cell death characterized by overwhelming lipid peroxidation. Beyond direct cytotoxicity, LF acts as a potent immunomodulator. It enhances natural killer (NK) cell activity, modulates T-lymphocyte populations, and crucially reprograms tumor-associated macrophages (TAMs) from a pro-tumor M2 state to an anti-tumor M1 state, thereby reversing the immunosuppressive tumor microenvironment (TME). The translation of LF’s potential has been significantly accelerated by nanotechnology. The inherent biocompatibility and natural tumor-targeting capabilities of LF make it an ideal platform for sophisticated drug-delivery systems. This review details various fabrication strategies for LF-based nanoparticles (NPs), including self-assembly, sol-in-oil emulsion, and electrostatic nanocomplexes, among others. Research demonstrates that nano-formulations not only protect LF from degradation but also enhance its bioactivity and anti-cancer potency. More importantly, LF NPs serve as versatile carriers for a wide array of therapeutic agents, including conventional chemotherapeutics, natural compounds, and imaging agents. These engineered systems enable synergistic therapy and facilitate site-specific delivery. Notably, the ability of LF to bind to receptors on the blood-brain barrier (BBB) has been leveraged to develop nano-systems for glioblastoma treatment. Other innovative designs utilize LF to modulate the TME—for instance, by alleviating tumor hypoxia to sensitize cells to radiotherapy and chemotherapy. Despite compelling pre-clinical evidence, the clinical translation of LF and its nano-formulations remains nascent. While early-phase trials have established a favorable safety profile for recombinant human LF, larger Phase III studies have yielded mixed results, underscoring the complexity of its action in humans. Key challenges include enhancing drug targeting, optimizing loading efficiency, ensuring batch-to-batch reproducibility, and achieving deep tumor penetration. Future research must focus on the rational design of next-generation LF-NPs. This entails developing standardized manufacturing protocols, engineering “smart” stimuli-responsive systems for targeted drug release in the TME, and constructing multi-targeting platforms. A concerted interdisciplinary effort is paramount to bridge the gap between bench and bedside. In conclusion, LF, particularly in its nano-engineered forms, represents a highly promising and versatile agent in the oncological arsenal, holding immense potential for precise and effective cancer therapy.

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.

Cancer remains a leading cause of global mortality, necessitating the development of advanced therapeutic strategies with enhanced efficacy and reduced systemic toxicity. Among promising bioactive agents, lactoferrin (LF)—a multifunctional iron-binding glycoprotein abundantly found in mammalian milk and exocrine secretions—has garnered significant interest for its potent and multifaceted anti-cancer properties. This review provides a comprehensive analysis of the current understanding of LF’s role in oncology, encompassing its structural biology, diverse mechanisms of action, and groundbreaking advancements in its application through nano-engineering. LF exerts anti-tumor effects through multiple pathways, including extracellular action, intracellular action, and immune regulation. It demonstrates a remarkable affinity for cancer cell membranes, binding to overexpressed anionic components such as glycosaminoglycans and sialic acids, as well as to specific receptors including the low-density lipoprotein receptor-related protein-1 (LRP-1). This selective binding facilitates targeted uptake. Upon internalization, LF orchestrates a direct assault by inducing cell-cycle arrest in phases such as G0/G1 or S phase through the modulation of key regulators including cyclins, CDKs, and p53. Furthermore, it promotes programmed cell death via apoptotic pathways, involving caspase activation and downregulation of anti-apoptotic proteins such as survivin. A more recently elucidated mechanism is the induction of ferroptosis, an iron-dependent form of cell death characterized by overwhelming lipid peroxidation. Beyond direct cytotoxicity, LF acts as a potent immunomodulator. It enhances natural killer (NK) cell activity, modulates T-lymphocyte populations, and crucially reprograms tumor-associated macrophages (TAMs) from a pro-tumor M2 state to an anti-tumor M1 state, thereby reversing the immunosuppressive tumor microenvironment (TME). The translation of LF’s potential has been significantly accelerated by nanotechnology. The inherent biocompatibility and natural tumor-targeting capabilities of LF make it an ideal platform for sophisticated drug-delivery systems. This review details various fabrication strategies for LF-based nanoparticles (NPs), including self-assembly, sol-in-oil emulsion, and electrostatic nanocomplexes, among others. Research demonstrates that nano-formulations not only protect LF from degradation but also enhance its bioactivity and anti-cancer potency. More importantly, LF NPs serve as versatile carriers for a wide array of therapeutic agents, including conventional chemotherapeutics, natural compounds, and imaging agents. These engineered systems enable synergistic therapy and facilitate site-specific delivery. Notably, the ability of LF to bind to receptors on the blood-brain barrier (BBB) has been leveraged to develop nano-systems for glioblastoma treatment. Other innovative designs utilize LF to modulate the TME—for instance, by alleviating tumor hypoxia to sensitize cells to radiotherapy and chemotherapy. Despite compelling pre-clinical evidence, the clinical translation of LF and its nano-formulations remains nascent. While early-phase trials have established a favorable safety profile for recombinant human LF, larger Phase III studies have yielded mixed results, underscoring the complexity of its action in humans. Key challenges include enhancing drug targeting, optimizing loading efficiency, ensuring batch-to-batch reproducibility, and achieving deep tumor penetration. Future research must focus on the rational design of next-generation LF-NPs. This entails developing standardized manufacturing protocols, engineering “smart” stimuli-responsive systems for targeted drug release in the TME, and constructing multi-targeting platforms. A concerted interdisciplinary effort is paramount to bridge the gap between bench and bedside. In conclusion, LF, particularly in its nano-engineered forms, represents a highly promising and versatile agent in the oncological arsenal, holding immense potential for precise and effective cancer therapy.

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.

The Type III Secretion System (T3SS) serves as a pivotal virulence apparatus for numerous Gram-negative bacterial pathogens, enabling them to infect both animal and plant hosts. Functioning as a molecular syringe, the T3SS directly translocates bacterial effector proteins from the bacterial cytoplasm into the interior of eukaryotic host cells. These effectors are central weapons that precisely manipulate a wide spectrum of host cellular physiological processes, ranging from cytoskeletal dynamics to immune signaling, to establish a favorable niche for bacterial survival and proliferation. Among the diverse arsenal of T3SS effectors, the YopJ family constitutes a critical group of virulence factors. Members of this family are characterized by a conserved catalytic triad structure—a hallmark of the CE clan of cysteine proteases that has been evolutionarily repurposed to confer acetyltransferase activity. A defining and intriguing feature of these enzymes is their stringent dependence on a host-derived eukaryotic cofactor, inositol hexakisphosphate (IP₆), for allosteric activation. This requirement acts as a sophisticated molecular safeguard, ensuring enzymatic activity only within the appropriate host environment, thereby preventing detrimental effects on the bacterium itself. While seminal studies on individual members such as Yersinia’s YopJ and Salmonella’s AvrA have provided deep mechanistic insights, a systematic and integrative understanding of the structure-function relationships across the entire family remains fragmented. Key questions persist regarding how a conserved catalytic core has diverged to recognize distinct host substrates in different kingdoms of life. To address this gap, this article provides a systematic review of the YopJ family, focusing on three interconnected aspects: their structural features, their catalytic mechanism, and their divergent immunosuppressive strategies in animal versus plant hosts. By conducting a comparative analysis of the sequences and resolved three-dimensional structures of three representative members (e.g., HopZ1a, PopP2, AvrA), we elucidate regions of significant variation embedded within the conserved core catalytic architecture. These variable regions, often involving surface loops and substrate-binding interfaces, are crucial determinants of target specificity and functional specialization. The functional divergence of this effector family is most apparent when comparing their modes of action in different hosts. In animal hosts, YopJ-family effectors primarily sabotage innate immune signaling pathways. They achieve this by acetylating key serine and threonine residues within the activation loops of critical kinases in the MAPK and NF‑κB pathways. This post-translational modification blocks the phosphorylation and subsequent activation of these kinases, leading to potent suppression of inflammatory cytokine production. Conversely, in plant hosts, the strategy broadens to dismantle the two-tiered plant immune system. YopJ homologs target a more diverse set of substrates, including immune-associated receptor-like cytoplasmic kinases (RLCKs), microtubule networks via tubulin acetylation (which disrupts cellular trafficking and signaling), and transcription factors central to defense gene regulation. This multi-target approach effectively suppresses both Pattern-Triggered Immunity (PTI) and Effector-Triggered Immunity (ETI). In conclusion, this synthesis aims to deepen the mechanistic understanding of YopJ family-mediated pathogenesis by integrating structural biology with cellular function across host kingdoms. Elucidating the precise molecular basis for substrate selection—how conserved platforms achieve target diversity—is a major frontier. Furthermore, this knowledge provides a vital theoretical foundation for developing novel anti-virulence strategies. Targeting the conserved IP₆-binding pocket or the catalytic acetyltransferase activity itself represents a promising avenue for designing broad-spectrum inhibitors that could disarm this critical family of bacterial effectors, potentially offering new therapeutic approaches against a range of pathogenic bacteria.

OBJECTIVE To assess the bleeding risk signals associated with the concomitant use of direct oral anticoagulants （DOACs） and triazole antifungals， and to provide pharmacovigilance evidence for the safety evaluation and monitoring of combined clinical use. METHODS Adverse event reports involving the concomitant use of DOACs and triazole antifungals were extracted from the US FDA Adverse Event Reporting System （FAERS） from the first quarter of 2004 to the third quarter of 2025. Nine bleeding-related preferred terms （PTs） were selected. The Ω shrinkage measure， additive model， multiplicative model， and combined risk ratio method were employed to detect drug-drug interaction signals. The strength of positive signals was further analyzed based on the Ω shrinkage measure. RESULTS A total of 790 adverse event reports involving the concomitant use of DOACs and triazole antifungals were included， among which 229 reports involved nine bleeding-related PTs. A total of 13 signals were consistently identified as posit ive by all four methods. These signals involved six drug combinations： apixaban-fluconazole， apixaban-posaconazole， rivaroxaban-itraconazole， dabigatran etexilate-fluconazole， apixaban-voriconazole， and dabigatran etexilate-itraconazole. The Ω shrinkage measure showed that the apixaban-posaconazole combination exhibited stronger signals for bleeding （ Ω ＝2.73， Ω 025 ＝2.05） and hemoptysis （ Ω ＝2.17， Ω 025 ＝0.83）； the apixaban-fluconazole combination exhibited stronger signals for hematoma （ Ω ＝2.30， Ω 025 ＝1.47） and hematuria （ Ω ＝1.71， Ω 025 ＝0.74）； the rivaroxaban-itraconazole combination exhibited stronger signals for epistaxis （ Ω ＝2.01， Ω 025 ＝0.90） and hematoma （ Ω ＝1.93， Ω 025 ＝0.42）； no positive Ω signals were observed for intracranial hemorrhage or upper gastrointestinal hemorrhage. CONCLUSION S This study suggests that the concomitant use of DOACs and triazole antifungals may increase the risk of bleeding-related events， with differences in signal strength and signal distribution across various drug combinations. In clinical practice， particular attention should be paid to the concomitant use of apixaban or rivaroxaban with strong cytochrome P450 3A4 or P-glycoprotein inhibitors such as posaconazole and itraconazole. For other DOAC-triazole antifungal combinations， close monitoring for bleeding-related manifestations and timely adjustment of anticoagulation or antifungal regimens are also warranted.