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.

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: CRISPR-Cas protects bacteria from exogenous DNA invasion and is associated with bacterial biofilm formation and pathogenicity. METHODS: We analyzed the type I-F CRISPR-Cas system of Legionella pneumophila WX48, including Cas1, Cas2-Cas3, Csy1, Csy2, Csy3, and Cas6f, along with downstream CRISPR arrays. We explored the effects of the CRISPR-Cas system on the in vitro growth, biofilm-forming ability, and pathogenicity of L. pneumophila through constructing gene deletion mutants. RESULTS: The type I-F CRISPR-Cas system did not affect the in vitro growth of wild-type or mutant strains. The biofilm formation and intracellular proliferation of the mutant strains were weaker than those of the wild type owing to the regulation of type IV pili and Dot/Icm type IV secretion systems. In particular, Cas6f deletion strongly inhibited these processes. CONCLUSION The type I-F CRISPR-Cas system may reduce biofilm formation and intracellular proliferation in L. pneumophila.
Legionella pneumophila/pathogenicity* ; CRISPR-Cas Systems ; Biofilms/growth & development* ; Phenotype ; Bacterial Proteins/metabolism* ; Gene Deletion

Objective To develop a predictive model for postoperative mortality risk in patients with acute aortic dissection(AAD)using the Extreme Gradient Boosting(XGBoost)algorithm combined with Shapley Additive Explanation(SHAP),and to establish a prediction website to serve as a diagnostic and therapeutic support platform for clinicians and patients.Methods A retrospective cohort study design was adopted.Data from 782 AAD patients who underwent surgical treatment at the Fourth Hospital of Hebei Medical University from January 2013 to December 2023 were collected,including basic information and initial serum biomarker test results.Patients were randomly divided into training and test sets at a 7:3 ratio.An external validation set consisting of 313 AAD patients admitted to the Second Hospital of Hebei Medical University from January 2020 to December 2023 was also established for further model validation.Variables were screened using LASSO regression,and an XGBoost machine learning model was constructed and interpreted using SHAP.The predictive performance of the model was evaluated using receiver operating characteristic(ROC)curve analysis.Using the Shiny package,the XGBoost model was deployed to shinyapps.io to create a prediction website for postoperative mortality risk in AAD patients.One patient was selected by simple random sampling from the test set and the external validation set respectively for the prediction example on the Shiny webpage.Results The XGBoost model demonstrated high predictive performance for postoperative mortality in AAD patients,with area under the ROC curve(AUC)values of 0.928(95%CI 0.901-0.956)in the training set,0.919(95%CI 0.891-0.949)in the test set,and 0.941(95%CI 0.915-0.967)in the external validation set.SHAP values indicated the following order of variable importance in the model(from highest to lowest):"lactate dehydrogenase""blood chlorine""multiple organ injury""carbon dioxide combining power""prothrombin time""α-hydroxybutyric acid""creatine kinase isoenzyme""Stanford classification""combined use of bedside blood purification""gender""acute kidney injury""gastrointestinal bleeding""brain injury"and"shock".A risk prediction website for adverse postoperative outcomes in AAD patients was developed using XGBoost-SHAP method(https://dun-dunxiaolu.shinyapps.io/document/)and validated with examples.One randomly selected patient from each of the test and external validation sets was applied:the predicted mortality risk value for patient 1(who died postoperatively)was 0.9539,and that for patient 2(who survived postoperatively)was 0.0206.Conclusions The XGBoost-SHAP model demonstrates high accuracy in predicting postoperative mortality risk for AAD patients.The online prediction tool established based on this model enhances the identification efficiency of high-risk postoperative mortality patients.

Objective To evaluate the efficacy and safety of transcatheter aortic valve implantation(TAVI)for the treatment of primary aortic valve regurgitation(NAVR)and to compare the difference in the choice of prosthetic valve size and the difference in complications with aortic stenosis(AS).Methods According to the definition of Valve Academic Research Consortium(VARC-3),143 patients with NAVR/AS treated with TAVI and patients with NAVR treated with surgical aortic valve replacement(SAVR)at Tongji Hospital,Tongji Medical College,Huazhong University of Science and Technology,China,from March 2019 to September 2024 were selected,and clinical data on baseline,perioperative,and primary endpoint events were were retrospectively collected and compared.Results Forty-three patients with NAVR were treated with TAVI,with a device success rate of 86.0％and a surgical success rate of 95.3％.Subgroup comparisons:(1)NAVR-TAVI group than NAVR-SAVR group:patients in the TAVI group had a significantly shorter operative time than those in the SAVR group(P＜0.001);complete left bundle branch block was more likely to occur after TAVI(P=0.042),and complete right bundle branch block was more likely to occur after SAVR(P=0.044).SAVR postoperatively The incidence of congestive heart failure was higher(P=0.013),and the mortality rate was significantly higher in the SAVR group than in the TAVI group(P=0.019).(2)NAVR-TAVI group than AS-TAVI group:the differences in access selection,THV size[28(22,34)mm vs.24(22,32)mm,P=0.044]and proportion of THV overdiameter[14％(7％,20％)vs.7％(3％,11％),P＜0.001]were statistically significant.patients in AS and NAVR groups had 1 case of permanent pacing after TAVI treatment.In the AS and NAVR groups,there was 1 case of permanent pacemaker implantation after TAVI.2 patients in the AS group were converted to surgical treatment,and 6 patients died.Conclusions The use of"off-label"(transfemoral)and"on-label"(transapical)TAVI devices(both from domestic sources)is safer than SAVR for the treatment of NAVR,especially in elderly and high-risk patients.Compared with patients with AS treated with TAVI,larger diameter annulas are usually selected for NAVR,with higher rates of valve migration,but overall safety and efficacy are comparable to AS.

Objective To evaluate the efficacy and safety of transcatheter aortic valve implantation(TAVI)for the treatment of primary aortic valve regurgitation(NAVR)and to compare the difference in the choice of prosthetic valve size and the difference in complications with aortic stenosis(AS).Methods According to the definition of Valve Academic Research Consortium(VARC-3),143 patients with NAVR/AS treated with TAVI and patients with NAVR treated with surgical aortic valve replacement(SAVR)at Tongji Hospital,Tongji Medical College,Huazhong University of Science and Technology,China,from March 2019 to September 2024 were selected,and clinical data on baseline,perioperative,and primary endpoint events were were retrospectively collected and compared.Results Forty-three patients with NAVR were treated with TAVI,with a device success rate of 86.0％and a surgical success rate of 95.3％.Subgroup comparisons:(1)NAVR-TAVI group than NAVR-SAVR group:patients in the TAVI group had a significantly shorter operative time than those in the SAVR group(P＜0.001);complete left bundle branch block was more likely to occur after TAVI(P=0.042),and complete right bundle branch block was more likely to occur after SAVR(P=0.044).SAVR postoperatively The incidence of congestive heart failure was higher(P=0.013),and the mortality rate was significantly higher in the SAVR group than in the TAVI group(P=0.019).(2)NAVR-TAVI group than AS-TAVI group:the differences in access selection,THV size[28(22,34)mm vs.24(22,32)mm,P=0.044]and proportion of THV overdiameter[14％(7％,20％)vs.7％(3％,11％),P＜0.001]were statistically significant.patients in AS and NAVR groups had 1 case of permanent pacing after TAVI treatment.In the AS and NAVR groups,there was 1 case of permanent pacemaker implantation after TAVI.2 patients in the AS group were converted to surgical treatment,and 6 patients died.Conclusions The use of"off-label"(transfemoral)and"on-label"(transapical)TAVI devices(both from domestic sources)is safer than SAVR for the treatment of NAVR,especially in elderly and high-risk patients.Compared with patients with AS treated with TAVI,larger diameter annulas are usually selected for NAVR,with higher rates of valve migration,but overall safety and efficacy are comparable to AS.

Objective To find potential biomarkers and analyzing metabolic pathways of the treatment by honey-processed Hedysari Radix,the cecal contents of rats with spleen-Qi deficiency were used as samples for analysis.Methods Sixty male SD rats were randomly divided into blank,model,experimental and control groups.The rats in other groups except the control group were carried out by using the three-factor compound modeling method of bitter-cold diarrhea,excessive exertion and hunger and satiety disorders.Experimental group was given 12.60 g·kg-1 honey-processed Hedysari Radix;control group was given 0.63 g·kg-1 lactobacillus bifidum triplex tabletsa;control and model groups received with equal volume of distilled water for a total of 15 days.Measure body weight,anal temperature,immune organ index of rats.Ultra-pressure liquid chromatography-quadrupole-exactive-mass spectrometry technology was used to measure the levels of endogenous metabolites in cecum contents.Orthogonal partial least squares discriminant analysis and database"Kyoto Encyclopedia of Genes and Genomes"were used to identify potential differential metabolites and possible metabolic pathways.Results After the intervention,the average body weight of the experimental,control,model and blank groups was(216.87±7.85),(210.96±9.03),(159.47±5.18)and(293.51±22.98)g;anal temperature was(36.14±0.48),(35.40±0.64),(34.50±0.78)and(36.61±0.34)℃;the thymus indexes were(1.19±0.20),(1.24±0.25),(0.47±0.15)and(1.31±0.21)mg·g-1;the spleen indexes were(1.95±0.33),(2.18±0.28),(1.61±0.27)and(2.29±0.24)mg·g-1.Compared with the model group,the above indexes of the experimental group and the control group were significantly increased(all P＜0.01).A total of 14 potential biomarkers of Honey-processed Hedysari Radix in treating spleen-Qi deficiency syndrome were screened out in this study,which mainly involved amino acid metabolism such as tryptophan and glutamate,riboflavin metabolism and adenosine 5'-monophosphate-activated protein kinase metabolism.Conclusion Honey-processed Hedysari Radix can further protect the intestinal mucosal barrier and reduce the intestinal inflammatory response by improving the metabolic level of cecum contents in rats with spleen-Qi deficiency in cecum contents,thus exerting the effect of strengthening the spleen and tonifying the Qi.

Endo-beta-N-acetylglucosaminidase (ENGase) is widely distributed in various organisms. The first reported ENGase activity was detected in Diplococcus pneumoniae in 1971. The protein (Endo D) was purified and its peptide sequence was determined in 1974. Three ENGases (Endo F1-F3) were discovered in Flavobacterium meningosepticum from 1982 to 1993. After that, the activity was detected from different species of bacteria, yeast, fungal, plant, mice, human, etc. Multiple ENGases were detected in some species, such as Arabidopsis thaliana and Trichoderma atroviride. The first preliminary crystallographic analysis of ENGase was conducted in 1994. But to date, only a few ENGases structures have been obtained, and the structure of human ENGase is still missing. The currently identified ENGases were distributed in the GH18 or GH85 families in Carbohydrate-Active enZyme (CAZy) database. GH18 ENGase only has hydrolytic activity, but GH85 ENGase has both hydrolytic and transglycosylation activity. Although ENGases of the two families have similar (β/α)8-TIM barrel structures, the active sites are slightly different. ENGase is an effective tool for glycan detection andglycan editing. Biochemically, ENGase can specifically hydrolyze β‑1,4 glycosidic bond between the twoN-acetylglucosamines (GlcNAc) on core pentasaccharide presented on glycopeptides and/or glycoproteins. Different ENGases may have different substrate specificity. The hydrolysis products are oligosaccharide chains and a GlcNAc or glycopeptides or glycoproteins with a GlcNAc. Conditionally, it can use the two products to produce a new glycopeptides or glycoprotein. Although ENGase is a common presentation in cell, its biological function remains unclear. Accumulated evidences demonstrated that ENGase is a none essential gene for living and a key regulator for differentiation. No ENGase gene was detected in the genomes of Saccharomyces cerevisiae and three other yeast species. Its expression was extremely low in lung. As glycoproteins are not produced by prokaryotic cells, a role for nutrition and/or microbial-host interaction was predicted for bacterium produced enzymes. In the embryonic lethality phenotype of the Ngly1-deficient mice can be partially rescued by Engase knockout, suggesting down regulation of Engase might be a solution for stress induced adaptation. Potential impacts of ENGase regulation on health and disease were presented. Rabeprazole, a drug used for stomach pain as a proton inhibitor, was identified as an inhibitor for ENGase. ENGases have been applied in vitro to produce antibodies with a designated glycan. The two step reactions were achieved by a pair of ENGase dominated for hydrolysis of substrate glycoprotein and synthesis of new glycoprotein with a free glycan of designed structure, respectively. In addition, ENGase was also been used in cell surface glycan editing. New application scenarios and new detection methods for glycobiological engineering are quickly opened up by the two functions of ENGase, especially in antibody remodeling and antibody drug conjugates. The discovery, distribution, structure property, enzymatic characteristics and recent researches in topical model organisms of ENGase were reviewed in this paper. Possible biological functions and mechanisms of ENGase, including differentiation, digestion of glycoproteins for nutrition and stress responding were hypothesised. In addition, the role of ENGase in glycan editing and synthetic biology was discussed. We hope this paper may provide insights for ENGase research and lay a solid foundation for applied and translational glycomics.

Objective:To investigate the pathogenic role and possible mechanism of NADPH oxidase 4 (Nox4) in type 1 diabetic keratopathy mouse models.Methods:Forty Nox4 knockout ( Nox4-/-) heterozygous male mice were selected and 120 age- and sex-matched wild-type C57BL/6 ( Nox4+ /+ ) mice were selected as controls. Nox4-/- and Nox4+ /+ mice were randomized into diabetic group (DM group) and non-DM group by random number method.Type 1 DM model was established in DM groups by intraperitoneal injection of streptozotocin.The DM and non-DM groups of Nox4+ /+ mice were randomized into regular feed group and Nox4 inhibitor GKT137831 (GKT) supplementary feed group by random number method.At 16 weeks after modeling, tear secretion of mice in different groups was measured by the phenol red thread test.Corneal epithelial integrity was evaluated by fluorescent staining.Changes in corneal never fiber density were observed by the in vivo laser scanning confocal microscopy.Reactive oxygen species (ROS) products in corneal epithelium were assayed by CellROX staining.The expressions of E-Cadherin and nuclear factor-κB (NF-κB) proteins were detected by immunofluorescence staining.Central corneal nerve fiber density was examined by flatmount staining with TUBB3 antibody.The use and care of laboratory animals complied with ARVO statement.The study protocol was approved by Laboratory Animal Care Committee of Xi'an Jiaotong University (No.XJTULAC201301). Results:In Nox4+ /+ mice, the tear secretion was (2.40±1.18)mm/minute in DM group, which was significantly less than (5.30±1.02)mm/minute in non-DM group ( P<0.01).The tear secretion was (4.19±0.63)mm/minute in DM group of Nox4-/- mice, which was significantly more than that in DM group of Nox4+ /+ mice ( P<0.05).Significant difference was found between (2.23±0.83)mm/minute of regular feed group and (4.02±0.71)mm/minute of GKT supplementary feed group ( P<0.01).In Nox4+ /+ mice, the DM group showed significantly increased corneal staining score, reduced corneal nerve fiber density, increased fluorescence intensity of ROS in corneal epithelium, weakened fluorescence intensity of E-Cadherin protein expression, and enhanced fluorescence of NF-κB protein expression compared with non-DM group.In Nox4-/- mice and mice fed with GKT supplementary feed, the increased fluorescence of ROS and decreased fluorescence of E-Cadherin protein expression were seen in the corneal epithelium of the DM groups compared with non-DM groups.In Nox4-/- mice and mice fed with GKT supplementary feed, NF-κB protein fluorescence was weak in corneal epithelial cells in DM groups, which was similar to that in non-DM groups.Immunofluorescence staining of corneal flatmount showed that the density of TUBB3-stained nerve fibers in DM group of Nox4+ /+ mice was significantly lower than that in non-DM group of Nox4+ /+ mice, and there was no significant reduction of nerve fibers in the corneal stromal layer in DM group of Nox4-/- mice or mice fed with GKT supplementary feed. Conclusions:Nox4 is involved in the pathogenic process of diabetic keratopathy, and its mechanism may be related to oxidative stress-induced aggregation of ROS products and activation of NF-κB-mediated inflammatory responses.

Objective To investigate the risk factors for bronchopulmonary dysplasia(BPD)in twin preterm infants with a gestational age of<34 weeks,and to provide a basis for early identification of BPD in twin preterm infants in clinical practice.Methods A retrospective analysis was performed for the twin preterm infants with a gestational age of<34 weeks who were admitted to 22 hospitals nationwide from January 2018 to December 2020.According to their conditions,they were divided into group A(both twins had BPD),group B(only one twin had BPD),and group C(neither twin had BPD).The risk factors for BPD in twin preterm infants were analyzed.Further analysis was conducted on group B to investigate the postnatal risk factors for BPD within twins.Results A total of 904 pairs of twins with a gestational age of<34 weeks were included in this study.The multivariate logistic regression analysis showed that compared with group C,birth weight discordance of>25%between the twins was an independent risk factor for BPD in one of the twins(OR=3.370,95%CI:1.500-7.568,P<0.05),and high gestational age at birth was a protective factor against BPD(P<0.05).The conditional logistic regression analysis of group B showed that small-for-gestational-age(SGA)birth was an independent risk factor for BPD in individual twins(OR=5.017,95%CI:1.040-24.190,P<0.05).Conclusions The development of BPD in twin preterm infants is associated with gestational age,birth weight discordance between the twins,and SGA birth.