Potential mechanisms and drug prediction of Rheumatoid Arthritis and primary Sjögren's Syndrome: A public databases-based study.

Rheumatoid arthritis (RA) and primary Sjögren's syndrome (pSS) are the most common systemic autoimmune diseases, and they are increasingly being recognized as occurring in the same patient population. These two diseases share several clinical features and laboratory parameters, but the exact mechanism of their co-pathogenesis remains unclear. The intention of this study was to investigate the common molecular mechanisms involved in RA and pSS using integrated bioinformatic analysis. RNA-seq data for RA and pSS were picked up from the Gene Expression Omnibus (GEO) database. Co-expression genes linked with RA and pSS were recognized using weighted gene co-expression network analysis (WGCNA) and differentially expressed gene (DEG) analysis. Then, we screened two public disease-gene interaction databases (GeneCards and Comparative Toxicogenomics Database) for common targets associated with RA and pSS. The DGIdb database was used to predict therapeutic drugs for RA and pSS. The Human microRNA Disease Database (HMDD) was used to screen out the common microRNAs associated with RA and pSS. Finally, a common miRNA-gene network was created using Cytoscape. Four hub genes (CXCL10, GZMA, ITGA4, and PSMB9) were obtained from the intersection of common genes from WGCNA, differential gene analysis and public databases. Twenty-four drugs corresponding to hub gene targets were predicted in the DGIdb database. Among the 24 drugs, five drugs had already been reported for the treatment of RA and pSS. Other drugs, such as bortezomib, carfilzomib, oprozomib, cyclosporine and zidovudine, may be ideal drugs for the future treatment of RA patients with pSS. According to the miRNA-gene network, hsa-mir-21 may play a significant role in the mechanisms shared by RA and pSS. In conclusion, we identified commom targets as potential biomarkers in RA and pSS from publicly available databases and predicted potential drugs based on the targets. A new understanding of the molecular mechanisms associated with RA and pSS is provided according to the miRNA-gene network.

Integrating scRNA-seq and bulk RNA-seq to characterize infiltrating cells in the colorectal cancer tumor microenvironment and construct molecular risk models.

Colorectal cancer (CRC) is a malignancy that is both highly lethal and heterogeneous. Although the correlation between intra-tumoral genetic and functional heterogeneity and cancer clinical prognosis is well-established, the underlying mechanism in CRC remains inadequately understood. Utilizing scRNA-seq data from GEO database, we re-isolated distinct subsets of cells, constructed a CRC tumor-related cell differentiation trajectory, and conducted cell-cell communication analysis to investigate potential interactions across cell clusters. A prognostic model was built by integrating scRNA-seq results with TCGA bulk RNA-seq data through univariate, LASSO, and multivariate Cox regression analyses. Eleven distinct cell types were identified, with Epithelial cells, Fibroblasts, and Mast cells exhibiting significant differences between CRC and healthy controls. T cells were observed to engage in extensive interactions with other cell types. Utilizing the 741 signature genes, prognostic risk score model was constructed. Patients with high-risk scores exhibited a significant correlation with unfavorable survival outcomes, high-stage tumors, metastasis, and low responsiveness to chemotherapy. The model demonstrated a strong predictive performance across five validation cohorts. Our investigation involved an analysis of the cellular composition and interactions of infiltrates within the microenvironment, and we developed a prognostic model. This model provides valuable insights into the prognosis and therapeutic evaluation of CRC.

Distinct molecular subtypes of systemic sclerosis and gene signature with diagnostic capability.

Background: As Systemic Sclerosis (SSc) is a connective tissue ailment that impacts various bodily systems. The study aims to clarify the molecular subtypes of SSc, with the ultimate objective of establishing a diagnostic model that can inform clinical treatment decisions. Methods: Five microarray datasets of SSc were retrieved from the GEO database. To eliminate batch effects, the combat algorithm was applied. Immune cell infiltration was evaluated using the xCell algorithm. The ConsensusClusterPlus algorithm was utilized to identify SSc subtypes. Limma was used to determine differential expression genes (DEGs). GSEA was used to determine pathway enrichment. A support vector machine (SVM), Random Forest(RF), Boruta and LASSO algorithm have been used to select the feature gene. Diagnostic models were developed using SVM, RF, and Logistic Regression (LR). A ROC curve was used to evaluate the performance of the model. The compound-gene relationship was obtained from the Comparative Toxicogenomics Database (CTD). Results: The identification of three immune subtypes in SSc samples was based on the expression profiles of immune cells. The utilization of 19 key intersectional DEGs among subtypes facilitated the classification of SSc patients into three robust subtypes (gene_ClusterA-C). Gene_ClusterA exhibited significant enrichment of B cells, while gene_ClusterC showed significant enrichment of monocytes. Moderate activation of various immune cells was observed in gene_ClusterB. We identified 8 feature genes. The SVM model demonstrating superior diagnostic performance. Furthermore, correlation analysis revealed a robust association between the feature genes and immune cells. Eight pertinent compounds, namely methotrexate, resveratrol, paclitaxel, trichloroethylene, formaldehyde, silicon dioxide, benzene, and tetrachloroethylene, were identified from the CTD. Conclusion: The present study has effectively devised an innovative molecular subtyping methodology for patients with SSc and a diagnostic model based on machine learning to aid in clinical treatment. The study has identified potential molecular targets for therapy, thereby offering novel perspectives for the treatment and investigation of SSc.

A compendium of mitochondrial molecular characteristics provides novel perspectives on the treatment of rheumatoid arthritis patients.

Rheumatoid arthritis (RA) is an autoimmune disease that exhibits a high degree of heterogeneity, marked by unpredictable disease flares and significant variations in the response to available treatments. The lack of optimal stratification for RA patients may be a contributing factor to the poor efficacy of current treatment options. The objective of this study is to elucidate the molecular characteristics of RA through the utilization of mitochondrial genes and subsequently construct and authenticate a diagnostic framework for RA. Mitochondrial proteins were obtained from the MitoCarta database, and the R package limma was employed to filter for differentially expressed mitochondrial genes (MDEGs). Metascape was utilized to perform enrichment analysis, followed by an unsupervised clustering algorithm using the ConsensuClusterPlus package to identify distinct subtypes based on MDEGs. The immune microenvironment, biological pathways, and drug response were further explored in these subtypes. Finally, a multi-biomarker-based diagnostic model was constructed using machine learning algorithms. Utilizing 88 MDEGs present in transcript profiles, it was possible to classify RA patients into three distinct subtypes, each characterized by unique molecular and cellular signatures. Subtype A exhibited a marked activation of inflammatory cells and pathways, while subtype C was characterized by the presence of specific innate lymphocytes. Inflammatory and immune cells in subtype B displayed a more modest level of activation (Wilcoxon test P < 0.05). Notably, subtype C demonstrated a stronger correlation with a superior response to biologics such as infliximab, anti-TNF, rituximab, and methotrexate/abatacept (P = 0.001) using the fisher test. Furthermore, the mitochondrial diagnosis SVM model demonstrated a high degree of discriminatory ability in distinguishing RA in both training (AUC = 100%) and validation sets (AUC = 80.1%). This study presents a pioneering analysis of mitochondrial modifications in RA, offering a novel framework for patient stratification and potentially enhancing therapeutic decision-making.

Mendelian randomization analysis suggests no associations of herpes simplex virus infections with systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) characterized by immune dysfunction is possibly more vulnerable to herpes simplex virus (HSV) infection. The infection has been intensively considered a common onset and exacerbation of SLE. This study is aimed at elucidating the causal association between SLE and HSV. A bidirectional two-sample Mendelian Randomization (TSMR) analysis was systematically conducted to explore the causal effect of SLE and HSV on each other. The causality was estimated by inverse variance weighted (IVW), MR-Egger and weighted median methods based on the summary-level genome-wide association studies (GWAS) data from a publicly available database. Genetically proxied HSV infection exhibited no causal association with SLE in the forward MR analysis using IVW method (odds ratio [OR] = 0.987; 95% confidence interval [CI]: 0.891-1.093; p = 0.798), nor did HSV-1 IgG (OR = 1.241; 95% CI: 0.874-1.762; p = 0.227) and HSV-2 IgG (OR = 0.934; 95% CI: 0.821-1.062; p = 0.297). Similar null results with HSV infection (OR = 1.021; 95% CI: 0.986-1.057; p = 0.245), HSV-1 IgG (OR = 1.003; 95% CI: 0.982-1.024; p = 0.788) and HSV-2 IgG (OR = 1.034; 95% CI: 0.991-1.080; p = 0.121) were observed in the reverse MR where SLE served as the exposure. Our study demonstrated no causal association between the genetically predicted HSV and SLE.

[Construction of a Prognostic Model of Multiple Myeloma Based on Metabolism-Related Genes].

OBJECTIVE: To screen the prognostic biomarkers of metabolic genes in patients with multiple myeloma (MM), and construct a prognostic model of metabolic genes. METHODS: The histological database related to MM patients was searched. Data from MM patients and healthy controls with complete clinical information were selected for analysis.The second generation sequencing data and clinical information of bone marrow tissue of MM patients and healthy controls were collected from human protein atlas (HPA) and multiple myeloma research foundation (MMRF) databases. The gene set of metabolism-related pathways was extracted from Molecular Signatures Database (MSigDB) by Perl language. The biomarkers related to MM metabolism were screened by difference analysis, univariate Cox risk regression analysis and LASSO regression analysis, and the risk prognostic model and Nomogram were constructed. Risk curve and survival curve were used to verify the grouping effect of the model. Gene set enrichment analysis (GSEA) was used to study the difference of biological pathway enrichment between high risk group and low risk group. Multivariate Cox risk regression analysis was used to verify the independent prognostic ability of risk score. RESULTS: A total of 8 mRNAs which were significantly related to the survival and prognosis of MM patients were obtained (P<0.01). As molecular markers, MM patients could be divided into high-risk group and low-risk group. Survival curve and risk curve showed that the overall survival time of patients in the low-risk group was significantly better than that in the high risk group (P<0.001). GSEA results showed that signal pathways related to basic metabolism, cell differentiation and cell cycle were significantly enriched in the high-risk group, while ribosome and N polysaccharide biosynthesis signaling pathway were more enriched in the low-risk group. Multivariate Cox regression analysis showed that the risk score composed of the eight metabolism-related genes could be used as an independent risk factor for the prognosis of MM patients, and receiver operating characteristic curve (ROC) showed that the molecular signatures of metabolism-related genes had the best predictive effect. CONCLUSION: Metabolism-related pathways play an important role in the pathogenesis and prognosis of patients with MM. The clinical significance of the risk assessment model for patients with MM constructed based on eight metabolism-related core genes needs to be confirmed by further clinical studies.

A Compendium of Mucosal Molecular Characteristics Provides Novel Perspectives on the Treatment of Ulcerative Colitis.

BACKGROUND AND AIMS: Ulcerative colitis [UC] is a complex heterogeneous disease. This study aims to reveal the underlying molecular features of UC using genome-scale transcriptomes of patients with UC, and to develop and validate a novel stratification scheme. METHODS: A normalised compendium was created using colon tissue samples (455 patients with UC and 147 healthy controls [HCs]), covering genes from 10 microarray datasets. Upregulated differentially expressed genes [DEGs] were subjected to functional network analysis, wherein samples were grouped using unsupervised clustering. Additionally, the robustness of subclustering was further assessed by two RNA sequencing datasets [100 patients with UC and 16 HCs]. Finally, the Xgboost classifier was applied to the independent datasets to evaluate the efficacy of different biologics in patients with UC. RESULTS: Based on 267 upregulated DEGs of the transcript profiles, UC patients were classified into three subtypes [subtypes A-C] with distinct molecular and cellular signatures. Epithelial activation-related pathways were significantly enriched in subtype A [named epithelial proliferation], whereas subtype C was characterised as the immune activation subtype with prominent immune cells and proinflammatory signatures. Subtype B [named mixed] was modestly activated in all the signalling pathways. Notably, subtype A showed a stronger association with the superior response of biologics such as golimumab, infliximab, vedolizumab, and ustekinumab compared with subtype C. CONCLUSIONS: We conducted a deep stratification of mucosal tissue using the most comprehensive microarray and RNA sequencing data, providing critical insights into pathophysiological features of UC, which could serve as a template for stratified treatment approaches.

Characteristics of the Gut Microbiome and Its Relationship With Peripheral CD4+ T Cell Subpopulations and Cytokines in Rheumatoid Arthritis.

This study investigated the association between intestinal microbiota abundance and diversity and cluster of differentiation (CD)4+ T cell subpopulations, cytokine levels, and disease activity in rheumatoid arthritis RA. A total of 108 rheumatoid arthritis (RA) patients and 99 healthy control (HC) subjects were recruited. PICRUSt2 was used for functional metagenomic predictions. Absolute counts of peripheral CD4+ T cell subpopulations and cytokine levels were detected by flow cytometry and with a cytokine bead array, respectively. Correlations were analyzed with the Spearman rank correlation test. The results showed that the diversity of intestinal microbiota was decreased in RA patients compared to HCs. At the phylum level, the abundance of Firmicutes, Fusobacteriota, and Bacteroidota was decreased while that of Actinobacteria and Proteobacteria was increased and at the genus level, the abundance of Faecalibacterium, Blautia, and Escherichia-Shigella was increased while that of Bacteroides and Coprococcus was decreased in RA patients compared to HC subjects. The linear discriminant analysis effect size indicated that Bifidobacterium was the most significant genus in RA. The most highly enriched Kyoto Encyclopedia of Genes and Genomes pathway in RA patients was amino acid metabolism. The relative abundance of Megamonas, Monoglobus, and Prevotella was positively correlated with CD4+ T cell counts and cytokine levels; and the relative numbers of regulatory T cells (Tregs) and T helper (Th17)/Treg ratio were negatively correlated with disease activity in RA. These results suggest that dysbiosis of certain bacterial lineages and alterations in gut microbiota metabolism lead to changes in the host immune profile that contribute to RA pathogenesis.

Development and validation of a GRGPI model for predicting the prognostic and treatment outcomes in head and neck squamous cell carcinoma.

Background: Head and neck squamous cell carcinoma (HNSCC) is among the most lethal and most prevalent malignant tumors. Glycolysis affects tumor growth, invasion, chemotherapy resistance, and the tumor microenvironment. Therefore, we aimed at identifying a glycolysis-related prognostic model for HNSCC and to analyze its relationship with tumor immune cell infiltrations. Methods: The mRNA and clinical data were obtained from The Cancer Genome Atlas (TCGA), while glycolysis-related genes were obtained from the Molecular Signature Database (MSigDB). Bioinformatics analysis included Univariate cox and least absolute shrinkage and selection operator (LASSO) analyses to select optimal prognosis-related genes for constructing glycolysis-related gene prognostic index(GRGPI), as well as a nomogram for overall survival (OS) evaluation. GRGPI was validated using the Gene Expression Omnibus (GEO) database. A predictive nomogram was established based on the stepwise multivariate regression model. The immune status of GRGPI-defined subgroups was analyzed, and high and low immune groups were characterized. Prognostic effects of immune checkpoint inhibitor (ICI) treatment and chemotherapy were investigated by Tumor Immune Dysfunction and Exclusion (TIDE) scores and half inhibitory concentration (IC50) value. Reverse transcription-quantitative PCR (RT-qPCR) was utilized to validate the model by analyzing the mRNA expression levels of the prognostic glycolysis-related genes in HNSCC tissues and adjacent non-tumorous tissues. Results: Five glycolysis-related genes were used to construct GRGPI. The GRGPI and the nomogram model exhibited robust validity in prognostic prediction. Clinical correlation analysis revealed positive correlations between the risk score used to construct the GRGPI model and the clinical stage. Immune checkpoint analysis revealed that the risk model was associated with immune checkpoint-related biomarkers. Immune microenvironment and immune status analysis exhibited a strong correlation between risk score and infiltrating immune cells. Gene set enrichment analysis (GSEA) pathway enrichment analysis showed typical immune pathways. Furthermore, the GRGPIdel showed excellent predictive performance in ICI treatment and drug sensitivity analysis. RT-qPCR showed that compared with adjacent non-tumorous tissues, the expressions of five genes were significantly up-regulated in HNSCC tissues. Conclusion: The model we constructed can not only be used as an important indicator for predicting the prognosis of patients but also had an important guiding role for clinical treatment.

Trichostatin A ameliorates Alzheimer's disease-related pathology and cognitive deficits by increasing albumin expression and Aß clearance in APP/PS1 mice.

BACKGROUND: Alzheimer's disease (AD) is an intractable neurodegenerative disorder in the elderly population, currently lacking a cure. Trichostatin A (TSA), a histone deacetylase inhibitor, showed some neuroprotective roles, but its pathology-improvement effects in AD are still uncertain, and the underlying mechanisms remain to be elucidated. The present study aims to examine the anti-AD effects of TSA, particularly investigating its underlying cellular and molecular mechanisms. METHODS: Novel object recognition and Morris water maze tests were used to evaluate the memory-ameliorating effects of TSA in APP/PS1 transgenic mice. Immunofluorescence, Western blotting, Simoa assay, and transmission electron microscopy were utilized to examine the pathology-improvement effects of TSA. Microglial activity was assessed by Western blotting and transwell migration assay. Protein-protein interactions were analyzed by co-immunoprecipitation and LC-MS/MS. RESULTS: TSA treatment not only reduced amyloid ß (Aß) plaques and soluble Aß oligomers in the brain, but also effectively improved learning and memory behaviors of APP/PS1 mice. In vitro study suggested that the improvement of Aß pathology by TSA was attributed to the enhancement of Aß clearance, mainly by the phagocytosis of microglia, and the endocytosis and transport of microvascular endothelial cells. Notably, a meaningful discovery in the study was that TSA dramatically upregulated the expression level of albumin in cell culture, by which TSA inhibited Aß aggregation and promoted the phagocytosis of Aß oligomers. CONCLUSIONS: These findings provide a new insight into the pathogenesis of AD and suggest TSA as a novel promising candidate for the AD treatment.

Low-Dose Sirolimus Immunoregulation Therapy in Patients with Active Rheumatoid Arthritis: A 24-Week Follow-Up of the Randomized, Open-Label, Parallel-Controlled Trial.

BACKGROUND: We have reported previously the insufficient absolute number or functional defects of regulatory T cells (Tregs) in patients with rheumatoid arthritis (RA), challenging conventional unspecific immunosuppressive therapy. Sirolimus, a mTOR inhibitor, is reported to allow growth of functional Tregs; here, we investigated the efficacy of low-dose sirolimus combined with conventional immunosuppressants (sirolimus immunoregulation therapy) for RA treatment with lower side effects and better tolerance. METHODS: In this nonblinded and parallel-group trial, we randomly assigned 62 patients to receive conventional glucocorticoids and immunosuppressants with or without sirolimus at a dosage of 0.5 mg on alternate days for 24 weeks in a 2 : 1 ratio. The demographic features, clinical manifestations, and laboratory indicators including peripheral blood lymphocyte subgroups and CD4+T subsets were compared before and after the treatment. RESULTS: Finally, 37 patients in the sirolimus group and 18 in the conventional treated group completed the 6-month study. By 24 weeks, the patients with sirolimus experienced significant reduction in disease activity indicators including DAS28, ESR, and the number of tender joints and swollen joints (p < 0.001). Notably, they had a higher level of Tregs as compared with those with conventional therapy alone (p < 0.05), indicating that sirolimus could partly restore the reduced Tregs. Concomitantly, their usage of immunosuppressants for controlling disease activity was decreased as compared with the conventional group with no difference in blood routine, and liver and renal functions both before and after the treatment of sirolimus and between the two groups (p > 0.05). CONCLUSIONS: Low-dose sirolimus immunoregulatory therapy selectively upregulated Tregs and partly replaced the usage of immunosuppressants to control disease activity without overtreatment and evaluable side effect. Further study is required using a large sample of RA patients treated with sirolimus for a longer period. This trial is registered at the Chinese Clinical Trial Registry (http://www.chictr.org.cn/showproj.aspx?proj=17245).

[Bioinformatics analysis of genes related to poor prognosis of human hepatocellular carcinoma and its clinical significance].

OBJECTIVE: To screen genes associated with poor prognosis of hepatocellular carcinoma (HCC) and to explore the clinical significance of these genes. METHODS: The proper expression profile data of HCC was obtained from the Gene Expression Omnibus (GEO) database, and the differentially expressed genes (DEGs) were identified by differential expression analysis. The DAVID and String database were used for function enrichment analysis and to construct the protein-protein interaction (PPI) network respectively. The Cancer Genome Atlas (TCGA) database and the Cox Proportional Hazard Model were used for prognosis analysis of the DEGs. RESULTS: A eligible human HCC data set (GSE84402) met the requirements. A total of 1141 differentially expressed genes were identified, including 720 up-regulated and 421 down-regulated genes. The results of function enrichment analysis and PPI network performed that CDK1、CDC6、CCNA2、CHEK1、CENPE 、PIK3R1、RACGAP1、BIRC5、KIF11 and CYP2B6 were prognosis key genes. And the prognosis analysis showed that the expressions of CDC6、PIK3R1、KIF11 and RACGAP1 were increased, and the expression of CENPE was decreased, which was closely related to prognosis of HCC. CONCLUSION: CDC6、CENPE、PIK3R1、KIF11 and RACGAP1 may be closely related to poor prognosis of HCC, and can be used as molecular biomarkers for future research of HCC prognosis.

[Analysis of Unfavorable Prognosis Gene Markers in Patients with Acute Myeloid Leukemia by Multiomics].

OBJECTIVE: To analyze the molecular markers associated with occurrence, development and poor prognosis of acute myeloid leukemia (AML) by using the data of GEO and TCGA database, as well as multiomics analysis. METHODS: The transcriptome data meeting requirements were down-loaded from GEO database, the differentially expressed genes were screened by using the R language limma package, and the GO function enrichment analysis and KEGG pathway analysis were performed for differentially expressed genes, at the same time, the protein interaction network was contracted by using STRING database and cytoscape software to screen out the hub gene, then the prognosis analysis was carried out for hub gene by combination with the clinical information affected in TCGA database. RESULTS: 620 differentially expressed genes were screened out, among which 162 differentially expressed genes were up-regulated, and 458 differentially expressed genes were down-regulated. Based on the results of GO functional enrichment, the KEGG pathway enrichment and protein interaction network, CXCL4, CXCR4, CXCR1, CXCR2, CCL5 and JUN were selected as hub genes. The survival analysis showed that the high expression of CXCL4, CXCR1, and CCL5 was a risk factor for poor prognosis of patiants. CONCLUSION: CXCL4, CXCR1 and CCL5 can be used as biomarkers for the occurrence and development of AML, which relateds with the unfavorable prognosis and can provide a basis for further study.

[Prognosis-related miRNA bioinformatics screening of lung adenocarcinoma and its clinical significance].

OBJECTIVE: To investigate the prognosis-related miRNA histological features and clinical significance of lung adenocarcinoma. METHODS: Using The Cancer Genome Atlas (TCGA) data, the miRNA expression profile data of human lung adenocarcinoma were searched for differential analysis, and the prognosis-related miRNAs were screened by Cox risk regression model. The targeted miRNAs were predicted by mirwalk analysis platform, KEGG functional enrichment analysis, and finally, predict the function of prognosis-related miRNAs. RESULTS: A total of 46 differential miRNAs in lung adenocarcinoma were screened, including 19 up-regulated and 27 down-regulated. Six prognostic-related miRNAs were screened by Cox survival analysis, namely hsa-mir-21, hsa-mir-142, hsa-mir-200a high expression, hsa-mir-101, hsa-let-7c, hsa-mir-378e low expression, hsa-mir-21 and hsa-mir-378e were associated with poor prognosis in patients with lung adenocarcinoma, and the survival time was shortened significantly (P<0.05, AUC=0.618). KEGG analysis showed that the above prognosis-related miRNA targeting regulatory genes were related with immune response pathways, miRNA and cancer pathways, metabolic pathways and so on. CONCLUSIONS: Hsa-mir-21 and hsa-mir-378e are associated with poor prognosis of lung adenocarcinoma, and may be used as a molecular marker for prognosis of lung adenocarcinoma after further clinical verification.

[Relationship between multi-gene proteins in colorectal carcinoma complicated with chronic schistosomiasis: an immunohistochemical study by using tissue microarray techniques].

OBJECTIVE: To investigate the relationship between p53, COX-2, Bax, c-myc genes and colorectal carcinoma complicated with chronic schistosomiasis. METHODS: One hundred and sixty patients with colorectal carcinoma were selected and divided into two groups; a schistosomiasis group (colorectal carcinoma complicated with chronic schistosomiasis, n = 80) and a non-schistosomiasis group (colorectal carcinoma uncomplicated with chronic schistosomiasis, n = 80). The tissue microarray techniques and immunohistochemistry method were used in all the patients to detect the expressions of p53, COX-2, Bax and c-myc proteins. RESULTS: The positive rate and level of p53 protein expression in the schistosomiasis group were lower than those in the non-schistosomiasis group, but there were no significant differences between the two groups (both P > 0.05). The COX-2 protein in both groups was positive, but the positive expression level of COX-2 in the schistosomiasis group was higher than that in the nonschistosomiasis group, and there was a significant difference between the two groups (P < 0.01). The positive rate and level of Bax protein expression were not significantly different between the two groups (both P > 0.05). The positive rate of c-myc expression in the schistosomiasis group was higher than that in the non-schistosomiasis group, with a significant difference (P < 0.01), but the positive expression level was lower than that in the non-schistosomiasis group, and there was a significant difference between the two groups (P < 0.01). CONCLUSIONS: Schistosome infection may impact on the deficiency of p53 of human colorectal cancer cells. It may promote the excessive expression of COX-2 protein, which is an indirect carcinogenic factor. The expression of Bax gene has no correlation with schistosome infection. The schistosome chronic infection may cause a persistent low level expression of c-myc gene.