Kyoto Encyclopedia Of Genes And Genomes Pathway Analysis Research Articles

Identification of potential biomarkers for colorectal cancer by clinical database analysis and Kaplan-Meier curves analysis.

This study aimed to explore critical genes as potential biomarkers for the diagnosis and prognosis of colorectal cancer (CRC) for clinical utility. To identify and screen candidate genes involved in CRC carcinogenesis and disease progression, we downloaded microarray datasets GSE89076, GSE73360, and GSE32323 from the GEO database identified differentially expressed genes (DEGs), and performed a functional enrichment analysis. A protein-protein interaction network was constructed, and correlated module analysis was performed using STRING and Cytoscape. The Kaplan-Meier survival curve shows the survival of the hub genes. The expression of cyclin-dependent kinase (CDK1), cyclin B1 (CCNB1), and PCNA in tissues and changes in tumor grade were analyzed. A total of 329 DEGs were identified, including 264 upregulated and 65 downregulated genes. The functions and pathways of DEGs include the mitotic cell cycle, poly(A) RNA binding replication, ATP binding, DNA replication, ribosome biogenesis in eukaryotes, and RNA transport. Forty-seven Hub genes were identified, and biological process analysis showed that these genes were mainly enriched in cell cycle and DNA replication. Patients with mutations in CDK1, PCNA, and CCNB1 had poorer survival rates. CDK1, PCNA, and CCNB1 were significantly overexpressed in the tumor tissues. The expression of CDK1 and CCNB1 gradually decreased with increasing tumor grade. CDK1, CCNB1, and PCNA can be used as potential markers for the diagnosis and prognosis of CRC. These genes are overexpressed in colon cancer tissues and are associated with low survival rates in CRC patients.

RNA-Seq-based transcriptomics analysis during the photodynamic therapy of primary cells in secondary hyperparathyroidism.

The aim of this study was to identify changes in gene expression before and after 5-aminolevulinic acid-mediated photodynamic therapy (5-ALA-PDT) and to investigate the potential mechanism of 5-ALA-PDT based on ribonucleic acid sequencing (RNA-Seq) analysis. Secondary hyperparathyroidism (SHPT) primary cells were isolated from surgically excised specimens and exposed to laser light. The transcription profiles of SHPT primary cells were identified through RNA-Seq. Differentially expressed genes (DEGs) were identified. Enrichment of functions and signaling pathway analysis were performed based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. Quantitative real-time polymerase chain reaction (RT-qPCR) and western blot analysis were used to validate genes based on RNA-Seq results. In total, 1320 DEGs were identified, of which 1019 genes were upregulated and 301 genes were downregulated. GO and KEGG pathway analyses identified significantly enriched pathways in DEGs, including TGF beta in extracellular matrix (ECM), negative regulation of triglyceride biosynthetic process, protein heterodimerization activity, systemic lupus erythematosus, ECM-receptor interaction, focal adhesion and protein digestion and absorption. Protein-protein interaction (PPI) network analyses identified potential heat shock protein (HSP) interactions among the DEGs. Eight HSP genes were also identified that were most likely involved in 5-ALA-PDT, which were further validated by RT-qPCR and western blotting. The findings of this descriptive study reveal changes in the transcriptome profile during 5-ALA-PDT, suggesting that gene expression and mutation, signaling pathways, and the molecular network are altered in SHPT primary cells. The above findings provide new insight for further studies on the mechanisms underlying 5-ALA-PDT in SHPT.

FNBP4 is a Potential Biomarker Associated with Cuproptosis and Promotes Tumor Progression in Hepatocellular Carcinoma.

Hepatocellular carcinoma (HCC) is one of the most common malignant tumors that lacks an efficient therapeutic approach because of its elusive molecular mechanisms. This study aimed to investigate the biological function and potential mechanism of formin-binding protein 4 (FNBP4) in HCC. FNBP4 expression in tissues and cells were detected by quantitative real-time PCR (qRT‒PCR), Western blot, and immunohistochemistry (IHC). The Kaplan-Meier method was used to explore the correlation between the FNBP4 expression and clinical survival. MTT, EdU incorporation, colony formation, and Transwell assays were performed to evaluate the function of FNBP4 in cell proliferation and migration in vitro. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was used to explore the potential mechanism of FNBP4. The prognostic risk signature and nomogram were constructed to demonstrate the prognostic value of FNBP4. We found that FNBP4 was upregulated in patients with HCC and associated with poor overall survival (OS). Furthermore, knockdown of FNBP4 inhibited the proliferation and migration in HCC cells. Then, we performed a KEGG pathway analysis of the coexpressed genes associated with FNBP4 and found that FNBP4 may be associated with tumor-related signaling pathways and cuproptosis. We verified that FNBP4 could cause cell cycle progression and inactivation of the hippo signaling pathway. A prognostic risk signature containing three FNBP4-related differentially expressed cuproptosis regulators (DECRs) was established and can be used as an independent risk factor to evaluate the prognosis of patients with HCC. In addition, a nomogram including a risk score and clinicopathological factors was used to predict patient survival probabilities. FNBP4, as a potential biomarker associated with cuproptosis, promotes HCC cell proliferation and metastasis. We provide a new potential strategy for HCC treatment by targeting FNBP4.

Association of gut microbiota and SCFAs with finishing weight of Diannan small ear pigs.

Finishing weight is a key economic trait in the domestic pig industry. Evidence has linked the gut microbiota and SCFAs to health and production performance in pigs. Nevertheless, for Diannan small ear (DSE) pigs, a specific pig breed in China, the potential effect of gut microbiota and SCFAs on their finishing weight remains unclear. Herein, based on the data of the 16S ribosomal RNA gene and metagenomic sequencing analysis, we found that 13 OTUs could be potential biomarkers and 19 microbial species were associated with finishing weight. Among these, carbohydrate-decomposing bacteria of the families Streptococcaceae, Lactobacillaceae, and Prevotellaceae were positively related to finishing weight, whereas the microbial taxa associated with intestinal inflammation and damage exhibited opposite effects. In addition, interactions of these microbial species were found to be linked with finishing weight for the first time. Gut microbial functional annotation analysis indicated that CAZymes, such as glucosidase and glucanase could significantly affect finishing weight, given their roles in increasing nutrient absorption efficiency. Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthologies (KOs) and KEGG pathways analysis indicated that glycolysis/gluconeogenesis, phosphotransferase system (PTS), secondary bile acid biosynthesis, ABC transporters, sulfur metabolism, and one carbon pool by folate could act as key factors in regulating finishing weight. Additionally, SCFA levels, especially acetate and butyrate, had pivotal impacts on finishing weight. Finishing weight-associated species Prevotella sp. RS2, Ruminococcus sp. AF31-14BH and Lactobacillus pontis showed positive associations with butyrate concentration, and Paraprevotella xylaniphila and Bacteroides sp. OF04-15BH were positively related to acetate level. Taken together, our study provides essential knowledge for manipulating gut microbiomes to improve finishing weight. The underlying mechanisms of how gut microbiome and SCFAs modulate pigs' finishing weight required further elucidation.

Identification of TRPC6 as a Novel Diagnostic Biomarker of PM-Induced Chronic Obstructive Pulmonary Disease Using Machine Learning Models.

Chronic obstructive pulmonary disease (COPD) was the third most prevalent cause of mortality worldwide in 2010; it results from a progressive and fatal deterioration of lung function because of cigarette smoking and particulate matter (PM). Therefore, it is important to identify molecular biomarkers that can diagnose the COPD phenotype to plan therapeutic efficacy. To identify potential novel biomarkers of COPD, we first obtained COPD and the normal lung tissue gene expression dataset GSE151052 from the NCBI Gene Expression Omnibus (GEO). A total of 250 differentially expressed genes (DEGs) were investigated and analyzed using GEO2R, gene ontology (GO) functional annotation, and Kyoto Encyclopedia of Genes and Genomes (KEGG) identification. The GEO2R analysis revealed that TRPC6 was the sixth most highly expressed gene in patients with COPD. The GO analysis indicated that the upregulated DEGs were mainly concentrated in the plasma membrane, transcription, and DNA binding. The KEGG pathway analysis indicated that the upregulated DEGs were mainly involved in pathways related to cancer and axon guidance. TRPC6, one of the most abundant genes among the top 10 differentially expressed total RNAs (fold change ≥ 1.5) between the COPD and normal groups, was selected as a novel COPD biomarker based on the results of the GEO dataset and analysis using machine learning models. The upregulation of TRPC6 was verified in PM-stimulated RAW264.7 cells, which mimicked COPD conditions, compared to untreated RAW264.7 cells by a quantitative reverse transcription polymerase chain reaction. In conclusion, our study suggests that TRPC6 can be regarded as a potential novel biomarker for COPD pathogenesis.

Identification and analysis of key microRNAs derived from osteoarthritis synovial fluid exosomes.

Identification and analysis of key microRNAs derived from osteoarthritis synovial fluid exosomes.

Influence of reactive oxygen species responsive self-assembled nanomicelle loaded with pyroptosis inhibitor on full-thickness skin defects in diabetic rats

Objective: To investigate the influence of reactive oxygen species (ROS) responsive self-assembled nanomicelle loaded with pyroptosis inhibitor on full-thickness skin defects in diabetic rats. Methods: Experimental research methods were employed. A nucleotide-binding oligomerization domain (NOD) 1/2 inhibitor (NOD-IN-1) was encapsulated with nanomicelle polyethylene glycol-block-polypropylene sulfide (PEG-b-PPS), and the resulting product was called PEPS@NOD-IN-1. The morphology and hydration particle size of PEG-b-PPS and PEPS@NOD-IN-1 were observed by transmission electron microscope and particle size analyzer, respectively, and the encapsulation rate and drug loading rate of PEPS@NOD-IN-1 to NOD-IN-1 and the cumulative release rate of NOD-IN-1 by PEPS@NOD-IN-1 in phosphate buffer solution (PBS) alone and hydrogen peroxide-containing PBS within 40 h were measured and calculated by microplate reader, and the sample number was 3. Twenty-four male Sprague-Dawley rats aged 6-7 weeks were injected with streptozotocin to induce type 1 diabetes mellitus. Six full-thickness skin defect wounds were made on the back of each rat. The injured rats were divided into PBS group, NOD-IN-1 group, PEG-b-PPS group, and PEPS@NOD-IN-1 group with corresponding treatment according to the random number table, with 6 rats in each group. The wound healing was observed on post injury day (PID) 3, 7, and 12, and the wound healing rate was calculated. The ROS levels in wound tissue were detected by immunofluorescence method on PID 3. On PID 7, the granulation tissue thickness in wound was assessed by hematoxylin-eosin staining, the mRNA expressions of NOD1 and NOD2 were detected by real-time fluorescence quantitative reverse transcription polymerase chain reaction, and the protein expressions of NOD1, NOD2, and GSDMD-N terminals were detected by Western blotting. Six wounds from different rats in each group were taken for detection of the above indicators. Wound tissue (3 samples per group) was taken from rats in PBS group and PEPS@NOD-IN-1 group on PID 7, and transcriptome sequencing was performed using high-throughput sequencing technology platform. Differentially expressed genes (DEGs) significantly down-regulated in PEPS@NOD-IN-1 group as compared with PBS group were screened, and the enrichment analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) was performed. The DEG heatmap of the NOD-like receptor pathway, a pyroptosis-related pathway, was made. Protein-protein interaction (PPI) analysis of DEGs in heatmap was performed through the STRING database to screen key genes of PEPS@NOD-IN-1 regulating the NOD-like receptor pathway. Data were statistically analyzed with analysis of variance for repeated measurement, one-way analysis of variance, and Tukey test. Results: PEG-b-PPS and PEPS@NOD-IN-1 were in spherical structures of uniform size, with hydration particle sizes of (134.2±3.3) and (143.1±2.3) nm, respectively. The encapsulation rate of PEPS@NOD-IN-1 to NOD-IN-1 was (60±5)%, and the drug loading rate was (15±3)%. The release of NOD-IN-1 from PEPS@NOD-IN-1 in PBS alone was slow, and the cumulative release rate at 40 h was only (12.4±2.3)%. The release of NOD-IN-1 from PEPS@NOD-IN-1 in hydrogen peroxide-containing PBS within 10 h was very rapid, and the cumulative release rate at 10 h reached (90.1±3.6)%. On PID 3 and 7, the wounds of rats in the four groups were gradually healed, and the healing in PEPS@NOD-IN-1 group was better than that in the other three groups. On PID 12, the wound scab area in PBS group was large, the wound epithelialization in NOD-IN-1 group and PEG-b-PPS group was obvious, and the wound in PEPS@NOD-IN-1 group was close to complete epithelialization. Compared with those in PBS group, NOD-IN-1 group, and PEG-b-PPS group, the wound healing rates on PID 7 and 12 in PEPS@NOD-IN-1 group were significantly increased (P<0.05), the level of ROS in wound tissue on PID 3 was significantly decreased (P<0.05), the thickness of granulation tissue in wound on PID 7 was significantly thickened (P<0.05), and the mRNA expressions of NOD1 and NOD2 and the protein expressions of NOD1, NOD2, and GSDMD-N terminals in wound tissue on PID 7 were significantly decreased (P<0.05). KEGG pathway analysis showed that DEGs significantly down-regulated in PEPS@NOD-IN-1 group as compared with PBS group were significantly enriched in NOD-like receptors, hypoxia-inducible factors, mitogen-activated protein kinases, and tumor necrosis factor (TNF) pathways. In the DEG heatmap of NOD-like receptor pathway, the genes regulating pyroptosis mainly involved NOD1, NOD2, NOD-like receptor thermoprotein domain-related protein 3, Jun, signal transduction and transcriptional activator 1 (STAT1), TNF-α-induced protein 3. The PPI results showed that NOD1, NOD2, and STAT1 were the key genes of PEPS@NOD-IN-1 regulating the NOD-like receptor pathway. Conclusions: PEPS@NOD-IN-1 can down-regulate the level of local ROS in wounds and the expression of NOD1, NOD2, and GSDMD-N terminals, the key regulators of pyroptosis, thereby promoting the repair of full-thickness skin defect wounds in diabetic rats. PEPS@NOD-IN-1 can also significantly down-regulate the pyroptosis, inflammation, and hypoxia-related pathways of wounds, and regulate NOD-like receptor pathways by down-regulating key genes NOD1, NOD2, and STAT1.

Identification and analysis of lncRNA, miRNA and mRNA related to subcutaneous and intramuscular fat in Laiwu pigs.

Long non-coding RNAs (lncRNAs) regulate adipocyte differentiation and metabolism, However, their function on subcutaneous and intramuscular adipose tissues in pigs is unclear. Intramuscular fat (IMF) is an important indicator for evaluating meat quality. Breeds with high IMF content are often accompanied by high subcutaneous fat (SCF), which severely affects the meat rate of pigs. It is of great significance for porcine breeding to study the mechanism of lncRNA related to adipogenesis and lipid metabolism. We identified differentially expressed lncRNAs, miRNAs and mRNAs in subcutaneous and intramuscular adipose tissues in three female Laiwu pigs by deep RNA-sequencing(|log2foldchange|≥1, P_value ≤ 0.05). The gene expression profiles of IMF and SCF in Laiwu pigs were comparatively analyzed by Bioinformatics methods to identify key lncRNAs, miRNAs, and mRNAs associated with lipid metabolism and adipogenesis. A total of 1209 lncRNAs (DElncRNAs), 286 miRNAs (DEmiRNAs), and 1597 mRNAs (DEgenes) were differentially expressed between two types of adipose. Among them, 17 DElncRNAs and 103 target genes play a role in the co-expression network, as well as 59 DElncRNAs, 44 DEmiRNAs, and 88 DEgenes involved in ceRNA network. In GO(Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis of DElncRNAs their target genes involved in many adipogenesis and lipid metabolism biological processes and signaling pathways, such as PPAR signaling pathway, Wnt signaling pathway, MAPK signaling pathway. By constructing co-expression and ceRNAs network we found that Wnt signaling pathway play a critical regulatory role in intramuscular adipogenesis and lipid accumulation in Laiwu pigs. TCONS_00006525, TCONS_00046551 and TCONS_00000528 may target WNT5A, WNT10B and FDZ3 in co-expression network, TCONS_00026517 and other lncRNAs regulate the expression of PPARG, RXRG and SCD in ceRNA network, and were involved in Wnt signaling pathway. This study provides a theoretical basis for further understanding the post-transcriptional regulation mechanism of meat quality formation, predicting and treating diseases caused by ectopic fat.

Differential expression profiles and functional prediction of circRNA in mice with traumatic heterotopic ossification.

Traumatic heterotopic ossification (HO) is an intractable sequela incited by inflammatory insult. To date, the exact molecular mechanisms of traumatic HO formation remain unclear. Recent studies have indicated that circular RNAs (circRNAs) participate in various human skeletal diseases. Although the formation of HO recapitulates many programs during bone development and remodeling, few data are available concerning whether circRNAs could participate in this pathological osteogenesis. To investigate the differentially expressed circRNAs (DE-circRNAs) in HO formation, microarray assay was performed to analyze the circRNA expression profile in four pairs of mice HO tissues and normal tissues. Then, qRT-PCR was applied to verify the microarray data. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed the biological functions of the differentially expressed circRNAs target genes. Cytoscape software was used to construct the circRNA-miRNA-mRNA network for circRNAs with different expression levels as well as the target genes. We demonstrated that 491 circRNAs were significantly differentially expressed in mouse HO tissues by a fold-change ≥ 2 and p-value ≤ 0.05. Among them, the expressions of 168 circRNAs were increased, while 323 were decreased. The expression levels of 10 selected circRNAs were verified successfully by qRT-PCR. GO analysis exhibited that these DE-circRNAs participated in a series of cellular processes. KEGG pathway analysis revealed that multiple upregulated and downregulated pathways were closely related to the DE-circRNAs in HO mice. The circRNA-miRNA-mRNA networks demonstrated that DE-circRNAs may be involved in the pathological osteogenesis of HO through the circRNA-targeted miRNA-mRNA axis. Our study first demonstrated the expression profiles and predicted the potential functions of DE-circRNAs in mice traumatic HO, which may shed new light on the elucidation of mechanisms as well as provide novel potential peripheral biological diagnostic markers and therapeutic targets for traumatic HO.

A five-pseudouridylation-associated-LncRNA classifier for primary prostate cancer prognosis prediction.

Background: Prostate cancer (PCa) is one of the most common cancers in males around the globe, and about one-third of patients with localized PCa will experience biochemical recurrence (BCR) after radical prostatectomy or radiation therapy. Reportedly, a proportion of patients with BCR had a poor prognosis. Cumulative studies have shown that RNA modifications participate in the cancer-related transcriptome, but the role of pseudouridylation occurring in lncRNAs in PCa remains opaque. Methods: Spearman correlation analysis and univariate Cox regression were utilized to determine pseudouridylation-related lncRNAs with prognostic value in PCa. Prognostic pseudouridylation-related lncRNAs were included in the LASSO (least absolute shrinkage and selection operator) regression algorithm to develop a predictive model. KM (Kaplan-Meier) survival analysis and ROC (receiver operating characteristic) curves were applied to validate the constructed model. A battery of biological cell assays was conducted to confirm the cancer-promoting effects of RP11-468E2.5 in the model. Results: A classifier containing five pseudouridine-related lncRNAs was developed to stratify PCa patients on BCR and named the "ψ-lnc score." KM survival analysis showed patients in the high ψ-lnc score group experienced BCR more than those in the low ψ-lnc score group. ROC curves demonstrated that ψ-lnc score outperformed other clinical indicators in BCR prediction. An external dataset, GSE54460, was utilized to validate the predictive model's efficacy and authenticity. A ceRNA (competitive endogenous RNA) network was constructed to explore the model's potential molecular functions and was annotated through GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analyses. RP11-468E2.5 was picked for further investigation, including pan-cancer analysis and experimental validation. Preliminarily, RP11-468E2.5 was confirmed as a tumor promoter. Conclusion: We provide some evidence that pseudouridylation in lncRNA played a role in the development of PCa and propose a novel prognostic classifier for clinical practice.

Titanium Dioxide Nanoparticles Altered the lncRNA Expression Profile in Human Lung Cells.

Respiration is considered to be the main occupational or environmental exposure pathway of titanium dioxide nanoparticles (TiO2 NPs), and the lung is considered to be the target organ of respiratory exposure; however, the mechanism of respiratory toxicity is not fully understood. In this study, the effect of TiO2 NPs on the expression profile of long non-coding RNA (lncRNA) in bronchial epithelial cells (BEAS-2B) was investigated to understand their potential toxic mechanism. BEAS-2B cells were treated with 100 μg/mL TiO2 NPs for 48 h, then RNA sequencing was performed to screen the differential lncRNAs compared with the control group, and the enrichment pathways of the differentially expressed lncRNAs were further analyzed using the Kyoto Encyclopedia of Genes and Genomes (KEGG). The results identified a total of 45,769 lncRNAs, and 277 different lncRNAs were screened. KEGG pathway analysis showed that the targeted mRNAs of these different lncRNAs were enriched in the pyrimidine metabolism pathway. This work demonstrates that TiO2 NPs could alter the lncRNA expression profile in BEAS-2B cells, and epigenetics may play a role in the mechanism of respiratory toxicity induced by TiO2 NPs.

System Analysis Based on Pancreatic Cancer Progression Identifies BRINP2 as a Novel Prognostic Biomarker.

Pancreatic adenocarcinoma (PAAD) is a malignant tumor of the digestive system, which develops rapidly and has no obvious early symptoms. This study aims to discover the biomarkers associated with PAAD development. We obtained RNA expression of PAAD patient samples and corresponding clinical data from The cancer genome atlas (TCGA), and screened out BMP/RA-inducible neural-specific protein 2 (BRINP2) gene which is highly associated with PAAD severity. Then, gene ontology (GO) enrichment, Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis and single-sample gene set enrichment analysis (ssGSEA) analysis were performed to explore the biological functions of BRINP2. Subsequently, long non-coding RNA (lncRNAs) associated with BRINP2 were screened out via correlation analysis, and Cox regression analysis and least absolute shrinkage selection operator (LASSO) regression analysis were used to construct the risk prediction model. We further validated the expression level of BRINP2 and its associated lncRNAs in BRINP2-associated lncRNAs prognostic model in vitro. We proposed that BRINP2 might be correlated to the tumor immune microenvironment and could also be used as a biomarker for PAAD progression. GO enrichment analysis and KEGG pathway analysis showed that the prognostic model was highly correlated to immune microenvironment-related pathways. Additionally, we established a BRINP2-associated lncRNAs prognostic model consisting of three lncRNAs. We validated the expression trends of BRINP2 and its associated lncRNAs in BRINP2-associated lncRNAs prognostic model in PAAD cells with various severity of metastatic potential using the quantitative real-time PCR (qRT-PCR). Meanwhile, pRRophetic R package was employed to predict potential therapeutic drugs for BRINP2-associated lncRNAs prognostic model of PAAD. The results suggest that BRINP2 can be used as a novel prognostic biomarker for PAAD.

Abstract P028: Methylation biomarkers for colorectal cancer early detection and survival prognostics impact gene expression and link to cancer-related biological pathways

Abstract Background: DNA methylation has been previously shown to have diagnostic and predictive potential for colorectal cancer (CRC). Aim of this study was to evaluate putative methylation markers in the context of early cancer development and diagnostics as well as further investigate the biological significance of these regions. Methods: Biomarker discovery was done by whole genome bisulfite sequencing (WGBS) of 88 CRC, 48 advanced adenoma (AA) and corresponding adjacent normal tissue (NAT) samples. Short-list of significantly hypermethylated regions (DMRs) was correlated to transcriptomics data from 512 CRC patients in The Cancer Genome Atlas (TCGA) cohort. Pathway enrichment for biological pathway analysis of the DMRs was done by using Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. Survival analysis was performed using Kaplan–Meier method on sub-groups of patients divided by the methylation status of individual markers. Finally, individual marker significance of selected regions was evaluated by analyzing 26 plasma samples from early stage (stage I-IIA) CRC samples and 42 colonoscopy verified controls (CNT) with targeted methylation sequencing assay. Results: 4167 putative marker regions were identified from biomarker discovery with WGBS. Differential signal could be observed both between AA and NAT and CRC and NAT, while several of these regions were differentially methylated also between AA and CRC samples, indicating biological signal change with adenoma progression to cancer. 84 hypermethylated DMRs from several verification studies were further evaluated against transcriptome data from TCGA, where overlap for 69 genes was found. 19 of these genes showed a significant down- regulation (p&amp;lt; 0.05), indicating a link between hypermethylation and gene expression. 2 genes showed significant up-regulation (p&amp;lt; 0.05), which could indicate other epigenic processes to be in place. KEGG pathway analysis revealed that the top pathways involved were axonal guidance, ephrin receptor signaling, epithelial-mesenchymal transition and FGF signaling, which all play significant role in the context of cancer development and progression. Kaplan-Meier analysis showed significant correlation to patients 5- year survival prediction linked to 3 genes: FGF14 (p=0.025, HR = 1.75) DPY19L2P1 (p=0.012, HR = 1.86), PTPRO (p=0.046, HR = 1.63). Targeted sequencing analysis on plasma samples of patients with early stage (I-IIA) colorectal cancer and age and gender matching colonoscopy-verified controls, showed high individual marker accuracy with AUC= 0.78 for FGF14, AUC= 0.81 for DPY19L2P1 and AUC= 0.73 for PTPRO. Conclusions: Methylation markers have distinct signals in early development of CRC, with high individual accuracy for separating early-stage cancers from matching controls. These regions have impact on gene expression and can be linked to relevant biological pathways. Extending early detection potential of the markers to further prognostics and stratification, could lead to better outcomes and improved survival of the patients. Citation Format: Pol Canal Noguer, Alejandro Requena Bermejo, Francesco Mattia Mancuso, Juan Carlos Higareda, Marina Manrique López, Marko Chersicola, Pablo Pérez Martínez, Pablo Antonio Camino, Primoz Knap, Vivian Erklavec Zajec, Kristi Kruusmaa. Methylation biomarkers for colorectal cancer early detection and survival prognostics impact gene expression and link to cancer-related biological pathways. [abstract]. In: Proceedings of the AACR Special Conference: Precision Prevention, Early Detection, and Interception of Cancer; 2022 Nov 17-19; Austin, TX. Philadelphia (PA): AACR; Can Prev Res 2023;16(1 Suppl): Abstract nr P028.

Mechanism Underlying the Action of Berberine in the Treatment of Gouty Arthritis Based on Network Pharmacology

Introduction: Gouty arthritis (GA) is induced by a purine metabolism disorder and monosodium urate (MSU) crystal-related inflammation. Berberine (BBR), extracted from Coptis chinensis, ameliorates MSU-induced GA. However, the mechanisms of BBR against GA remain to be fully elucidated. This study aimed to identify the key targets and pathways mediating the effects of BBR against GA using network pharmacology. Methods: BBR and GA targets were obtained from several databases, and the network of BBR-GA common targets was visualized using Cytoscape software. Protein–protein interaction (PPI), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed using the STRING and Database for Annotation, Visualization, and Integrated Discovery databases, respectively. Animal experiments were performed to determine the outcomes of the BBB intervention. The serum levels of IL-1β, IL-8, and IL-6 were detected using enzyme-linked immunosorbent assay. Results: Thirty-three common targets (including NF-κB, RelA, MAPK1, IL-6, and IL-1β) of BBR and GA were identified, and a network of common targets between BBR and GA was constructed. PPI analysis demonstrated that IL-1β, IL-6, TNF, MAPK, and RelA are key targets with high degree values. GO and KEGG pathway analyses revealed the involvement of inflammation-related biological processes and signaling pathways, such as the NF-κB, MAPK, and TNF signaling pathways. Animal experiments demonstrated that the uric acid, IL-1β, IL-6, and IL-8 serum levels were significantly lower in the BBR group compared with those in hyperuricemic rats. Conclusions: Using systematic network analysis, potential targets mediating the effects of BBR on GA were detected. The pathways and inflammatory factors involved were identified using in vivo experiments, thus providing a reference for further basic research and clinical applications of BBR in the treatment of GA.

The mechanism study of Miao medicine Tongfengting decoction in the treatment of gout based on network pharmacology and molecular docking.

This study sought to clarify the mechanism of action of Miao medicine Tongfengting decoction in the treatment of gout through network pharmacology and molecular docking by searching for its key targets and related pathways. The active ingredients of Miao medicine Tongfengting Decoction were obtained from the TCMSP data platform, searched the relevant databases for gout-related targets,using String and Cytoscape 3.9 to build a "compound-cross-target-disease" network diagram,performed gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis in the DAVID database, and performed the docking analysis using PyMoL 2.3.0 and AutoDock. After screening, 298 main targets of the Miao medicine Tongfengting decoction for gout were identified. The target network is established, and the topology of protein-protein interaction (PPI) network is analyzed. The enrichment analysis of KEGG pathway showed that these targets were related to Pathways in cancer, PI3K Akt signaling pathway, MAPK signaling pathway and other pathways. Molecular docking showed that the target protein had good binding power with the main active components of the compound of Miao medicine Tongfengting Decoction. Miao medicine Tongfengting decoction probably regulates immune mechanism using a multi-component, multi-target, and multi-pathway strategy to reduce inflammatory response and exert its therapeutic effect on gout.

Transcriptome profiling of male and female Ascaris lumbricoides reproductive tissues

BackgroundAscaris lumbricoides causes human ascariasis, the most prevalent helminth disease, infecting approximately 1 billion individuals globally. In 2019 the global disease burden was estimated to be 754,000 DALYs and resulted in 2090 deaths. In the absence of a vaccination strategy, treatment of ascariasis has relied on anthelminthic chemotherapy, but drug resistance is a concern. The propensity for reinfection is also a major challenge to disease control; female worms lay up to 200,000 eggs daily, which contaminate surrounding environments and remain viable for years, resulting in high transmission rates. Understanding the molecular mechanisms of reproductive processes, including control of egg production, spermatogenesis, oogenesis and embryogenesis, will drive the development of new drugs and/or vaccine targets for future ascariasis control.MethodsTranscriptome profiles of discrete reproductive and somatic tissue samples were generated from adult male and female worms using Illumina HiSeq with 2 × 150 bp paired-end sequencing. Male tissues included: testis germinal zone, testis part of vas deferens, seminal vesicle and somatic tissue. Female tissues included: ovary germinal zone, ovary part of the oviduct, uterus and somatic tissue. Differentially expressed genes (DEGs) were identified from the fragments per kilobases per million reads (FPKM) profiles. Hierarchical analysis was performed to identify tissue-specific genes. Furthermore, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were employed to identify significant terms and pathways for the DEGs.ResultsDEGs involved in protein phosphorylation and adhesion molecules were indicated to play a crucial role in spermatogenesis and fertilization, respectively. Those genes associated with the G-protein-coupled receptor (GPCR) signaling pathway and small GTPase-mediated signal transduction pathway play an essential role in cytoskeleton organization during oogenesis. Additionally, DEGs associated with the SMA genes and TGF-β signaling pathway are crucial in adult female embryogenesis. Some genes associated with particular biological processes and pathways that were identified in this study have been linked to defects in germline development, embryogenesis and reproductive behavior. In the enriched KEGG pathway analysis, Hippo signaling, oxytocin signaling and tight junction pathways were identified to play a role in Ascaris male and female reproductive systems.ConclusionsThis study has provided comprehensive transcriptome profiles of discrete A. lumbricoides reproductive tissue samples, revealing the molecular basis of these functionally important tissues. The data generated from this study will provide fundamental knowledge on the reproductive biology of Ascaris and will inform future target identification for anti-ascariasis drugs and/or vaccines.Graphical

Chemokine receptor 7 mediates miRNA-182 to regulate cerebral ischemia/reperfusion injury in rats.

Chemokine receptor 7 (CXCR7) exerts protective effects on the brain. MicroRNAs (miRNAs) are involved in cerebral ischemia/reperfusion (I/R) injury, but their involvement in CXCR7-mediated brain protection is unknown. In this study, we investigated the role of miRNAs in CXCR7-mediated brain protection. CXCR7 levels in peripheral blood samples from patients with acute ischemic stroke (AIS) and ischemic penumbra area brain tissues from middle cerebral artery occlusion (MCAO) rats after recanalization were measured. An miRNA microarray analysis was performed to examine the expression of miRNAs caused by CXCR7 knockdown in ischemic penumbra area brain tissue in middle cerebral artery occlusion-reperfusion rats and to predict corresponding downstream target genes. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis revealed the most enriched pathways. A dual-luciferase reporter assay confirmed the direct regulation of miR-182 on the target gene TCF7L2. The correlation between TCF7L2 and CXCR7/miR-182 was verified using rescue assays. CXCR7 expression was upregulated in MCAO rats and mechanical thrombectomy patients with AIS compared to that in controls. The motor and sensory functions of MCAO rats with CXCR7 knockdown further decreased, and the infarct volume and cerebral edema increased. miRNA microarray data showed that seven miRNAs were differentially expressed after shRNA-CXCR7 treatment. The dual-luciferase reporter assay confirmed that miR-182 directly targeted the TCF7L2 gene. Rescue assays confirmed that TCF7L2 is downstream of CXCR7/miR-182. KEGG pathway analysis showed that the Hippo pathway may be a key pathway in CXCR7 upregulation and plays a role in protecting the brain after interventional surgery. Animal experiments have shown that CXCR7-mediated cerebral I/R injury promotes the phosphorylation of key molecules YAP and TAZ in the Hippo pathway. CXCR7 protects against cerebral I/R injury, possibly via the miR-182/TCF7L2/Hippo pathway. These results indicate that CXCR7 affects cerebral ischemia-reperfusion injury through miRNA regulation and downstream pathways.

Network pharmacology and LC-MS approachs to explore the active compounds and mechanisms of Yuanjiang decoction for treating bradyarrhythmia

BackgroundYuanjiang decoction (YJD), a traditional Chinese medicinal prescription, has been found to have a significant heart rate-increasing effect and is effective in the treatment of symptomatic bradyarrhythmia in previous studies. However, its specific components and potential mechanisms remain unclear. MethodsIn this study, we detected and identified the main compounds of YJD using liquid chromatography-mass spectrometry (LC-MS). Through the approach of network pharmacology, we predicted the core targets of the active components, bradyarrhythmia targets, and obtained potential anti-bradyarrhythmia targets of YJD. We further performed protein to protein interaction (PPI), gene ontology (GO) enrichment analyses and kyoto encyclopedia of genes and genomes (KEGG) signaling pathway analyses for core targets, and constructed network of key active ingredients-core targets of YJD. Finally, molecular docking and molecular dynamics simulation were performed for key active ingredients and core targets. ResultsThe YJD contains a total of 35 main chemical components. The key active ingredients-core targets network contains 36 nodes and 90 edges, including 20 key active ingredients and 16 core targets. The core targets in the PPI network were TP53, TNF, HRAS, PPARG, IL1B, KCNH2, SCN5A, IDH1, LMNA, ACHE, F2, DRD2, CALM1, KCNQ1, TNNI3, IDH2 and TNNT2. KEGG pathway analysis showed that YJD treatment of bradyarrhythmia mainly involves neuroactive ligand-receptor interaction, adrenergic signaling in cardiomyocytes, cAMP signaling pathway, calcium signaling pathway, cholinergic synaptic and serotonergic synapse signaling pathway. The biological processes mainly include regulation of hormone levels, regulation of cardiac contraction, chemical synaptic transmission, circadian rhythm, positive regulation of heart rate, smooth muscle contraction, response to metal ion, oxidation-reduction process, neurotransmitter transport and import across plasma membrane. Molecular docking and molecular dynamics simulation results showed that hesperidin and tetrahydropalmatine had higher affinity with DRD2 and KCNQ1, respectively. ConclusionThis study reveals the pharmacodynamic material basis of YJD and its potential multicomponent–multitarget–multipathway pharmacological effects, predicted its potential anti-bradyarrhythmia mechanism may be related to the regulation of myocardial autonomic nervous function and related ion channels. Our work demonstrates that YJD has great potential for treating bradyarrhythmias as a complementary medicine, and the results can provide a theoretical basis for the development and clinical application of YJD.

Development and validation of the prognostic model based on autophagy-associated genes in idiopathic pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a chronic progressive interstitial lung disease. Many studies suggest that autophagy may be related to disease progression and prognosis in IPF. However, the mechanisms involved have not been fully elucidated. We incorporated 232 autophagy-associated genes (AAGs) and two datasets, GSE28042 and GSE27957, from the GEO database. Univariate Cox analysis and least absolute shrinkage and selection operator (LASSO) regression were used to construct the autophagy-associated prognostic model. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to investigate the functions of these autophagy-associated genes. CIBERSORT algorithm was used to calculate the immune cell infiltration between patients in the high-risk score and low-risk score groups. Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) was performed to explore the mRNA expression of five genes in the autophagy-associated risk model. We constructed a 5-autophagy-associated genes signature based on Univariate Cox analysis and LASSO regression. In our autophagy-associated risk model, IPF patients in the high-risk group demonstrated a poor overall survival rate compared to patients in the low-risk group. For 1-, 2-, and 3-year survival rates, the AUC predictive value of the AAG signature was 0.670, 0.787, and 0.864, respectively. These results were validated in the GSE27957 cohort, confirming the good prognostic effect of our model. GO and KEGG pathway analyses enriched immune-related pathways between the high-risk and low-risk groups. And there was also a significant difference in immune cell infiltration between two groups. And the results of qRT-PCR showed that the expression levels of FOXO1, IRGM, MYC, and PRKCQ were significantly decreased in the Peripheral Blood Mononuclear Cell (PBMC) of IPF patient samples. Our study constructed and validated an autophagy-associated risk model based on MYC, MAPK1, IRGM, PRKCQ, and FOXO1. And those five genes may influence the progression of IPF by regulating immune responses and immune cells.

HIF-1 signalling pathway was identified as a potential new pathway for Icariin's treatment against Alzheimer's disease based on preclinical evidence and bioinformatics.

Aim: Alzheimer's disease (AD) is a neurodegenerative condition that is characterized by the gradual loss of memory and cognitive function. Icariin, which is a natural chemical isolated from Epimedii herba, has been shown to protect against AD. This research examined the potential mechanisms of Icariin's treatment against AD via a comprehensive review of relevant preclinical studies coupled with network pharmacology. Methods: The PubMed, Web of Science, CNKI, WANFANG, and VIP databases were used to identify the relevant studies. The pharmacological characteristics of Icariin were determined using the SwissADME and TCMSP databases. The overlapping targets of Icariin and AD were then utilized to conduct disease oncology (DO) analysis to identify possible hub targets of Icariin in the treatment of AD. The hub targets were then used for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses, and the interactions of the targets and Icariin were assessed via molecular docking and molecular dynamics simulation (MDS). Results: According to the literature review, Icariin alleviates cognitive impairment by regulating the expression of Aβ1-42, Aβ1-40, BACE1, tau, hyperphosphorylated tau, and inflammatory mediators. DO analysis revealed 35 AD-related hub targets, and the HIF-1 signalling pathway was ranked first according to the KEGG pathway analysis. Icariin effectively docked with the 35 hub targets and HIF-1α, and the dynamic binding of the HIF-1-Icariin complex within 100 ns indicated that Icariin contributed to the stability of HIF-1α. Conclusion: In conclusion, our research used a literature review and network pharmacology methods to identify the HIF-1 signalling pathway as a potential pathway for Icariin's treatment against AD.