Gene Ontology Enrichment Analysis Research Articles

Circular RNA Expression Profile in Cervical Cancer and Construction of the Circular RNA‑MicroRNA‑Messenger RNA Regulatory Network

Introduction. Cervical cancer (CC) remains the most common cancer in women worldwide. However, effective and specific biomarkers for the diagnosis and prognosis of cervical cancer are yet to be found. In recent years, the potential of circular RNAs (circRNAs) as new diagnostic, prognostic and therapeutic tools has received much attention. The current study involved an in-depth bioinformatics research to explore the circRNA-microRNA (miRNA)-messenger RNA (mRNA) regulatory network in order to identify important molecular processes and biological pathways supposedly associated with CC. Materials and methods. The study collected data on the expression of circRNA (GSE102686), miRNA (GSE30656) and mRNA of target genes (GSE9750), based on the Gene Expression Omnibus (GEO) database, in squamous cell carcinoma of the cervix samples and normal squamous epithelium of the cervix, dividing them into study and control groups. Protein-protein interaction (PPI), Gene Ontology (GO) analysis, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed to further understand the function of circRNAs for their target genes. Results. A total of 105 differentially expressed circular RNAs (DECs), 144 differentially expressed microRNAs (DEMs), and 539 differentially expressed target genes (DEGs) were identified for cervical cancer. Concurrently, functional enrichment analysis of GO and KEGG pathways was performed for DEGs. Subsequently, searching databases for circRNA, miRNA and mRNA target genes, as well as PPI network analysis and functional enrichment revealed 3 DECs with significantly high expression levels (hsa_circ_0000745, hsa_circ_0084927 and hsa_circ_0002762), 6 DEMs with reduced expression levels (hsa -miR-145, hsa-miR-876-3p, hsa-miR-1229, hsa-miR-182, hsa-miR-520h and hsa-miR-1252) and 9 key genes such as ANGPT2, COL11A1, MEST, KIF20A, CLN6, FNDC3B, USP18, DLGAP5 and CXCL9, suggesting a potentially significant role in cervical cancer. Conclusion.Understanding the circRNA-miRNA-mRNA regulatory network is of great importance for evaluating the oncogenesis of CC, as well as discoverying new circRNAs as the main regulatory molecules in this network. This is considered to be a new direction in the diagnosis and targeted therapy of cervical cancer.

Identification of core genes shared by ischemic stroke and myocardial infarction using an integrated approach.

Both ischemic stroke (IS) and myocardial infarction (MI) are caused by vascular occlusion that results in ischemia. While there may be similarities in their mechanisms, the potential relationship between these 2 diseases has not been comprehensively analyzed. Therefore, this study explored the commonalities in the pathogenesis of IS and MI. Datasets for IS (GSE58294, GSE16561) and MI (GSE60993, GSE61144) were downloaded from the Gene Expression Omnibus database. Transcriptome data from each of the 4 datasets were analyzed using bioinformatics, and the differentially expressed genes (DEGs) shared between IS and MI were identified and subsequently visualized using a Venn diagram. A protein-protein interaction (PPI) network was constructed using the Interacting Gene Retrieval Tool database, and identification of key core genes was performed using CytoHubba. Gene Ontology (GO) term annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of the shared DEGs were conducted using prediction and network analysis methods, and the functions of the hub genes were determined using Metascape. The analysis revealed 116 and 1321 DEGs in the IS and MI datasets, respectively. Of the 75 DEGs shared between IS and MI, 56 were upregulated and 19 were downregulated. Furthermore, 15 core genes - S100a12, Hp, Clec4d, Cd163, Mmp9, Ormdl3, Il2rb, Orm1, Irak3, Tlr5, Lrg1, Clec4e, Clec5a, Mcemp1, and Ly96 - were identified. GO enrichment analysis of the DEGs showed that they were mainly involved in the biological functions of neutrophil degranulation, neutrophil activation during immune response, and cytokine secretion. KEGG analysis showed enrichment in pathways pertaining to Salmonella infection, Legionellosis, and inflammatory bowel disease. Finally, the core gene-transcription factor, gene-microRNA, and small-molecule relationships were predicted. These core genes may provide a novel theoretical basis for the diagnosis and treatment of IS and MI.

Time-course study of genetic changes in periodontal ligament regeneration after tooth replantation in a mouse model

This research focused on analyzing gene expression changes in the periodontal ligament (PDL) after tooth re-plantation to identify key genes and pathways involved in healing and regeneration. Utilizing a mouse model, mRNA was extracted from the PDL at various intervals post-replantation for RNA sequencing analysis, spanning from 3 to 56 days. The results revealed significant shifts in gene expression, particularly notable on day 28, supported by hierarchical clustering and principal component analysis. Gene ontology (GO) enrichment analysis highlighted an upregulation in olfactory receptor and G protein-coupled receptor signaling pathways at this time point. These findings were validated through reverse transcription-quantitative PCR (RT-qPCR), with immunochemical staining localizing olfactory receptor gene expression to the PDL and surrounding tissues. Moreover, a scratch assay indicated that olfactory receptor genes might facilitate wound healing in human PDL fibroblasts. These results underscore the importance of the 28-day post-transplant phase as a potential “tipping point” in PDL healing and regeneration. In conclusion, this research sheds light on the potential role of olfactory receptor genes in PDL regeneration, providing a foundation for developing new therapeutic approaches in tooth replantation and transplantation, with broader implications for regenerative medicine in oral health.

The role of N-acetylcysteine on adhesion and biofilm formation of Candida parapsilosis isolated from catheter-related candidemia.

Objectives. Anti-fungal agents are increasingly becoming less effective due to the development of resistance. In addition, it is difficult to treat Candida organisms that form biofilms due to a lack of ability of drugs to penetrate the biofilms. We are attempting to assess the effect of a new therapeutic agent, N-acetylcysteine (NAC), on adhesion and biofilm formation in Candida parapsilosis clinical strains. Meanwhile, to detect the transcription level changes of adhesion and biofilm formation-associated genes (CpALS6, CpALS7, CpEFG1 and CpBCR1) when administrated with NAC in C. parapsilosis strains, furthermore, to explore the mechanism of drug interference on biofilms.Hypothesis/Gap statement. N-acetylcysteine (NAC) exhibits certain inhibitory effects on adhesion and biofilm formation in C. parapsilosis clinical strains from CRBSIs through: (1) down-regulating the expression of the CpEFG1 gene, making it a highly potential candidate for the treatment of C. parapsilosis catheter-related bloodstream infections (CRBSIs), (2) regulating the metabolism and biofilm -forming factors of cell structure.Methods. To determine whether non-antifungal agents can exhibit inhibitory effects on adhesion, amounts of total biofilm formation and metabolic activities of C. parapsilosis isolates from candidemia patients, NAC was added to the yeast suspensions at different concentrations, respectively. Reverse transcription was used to detect the transcriptional levels of adhesion-related genes (CpALS6 and CpALS7) and biofilm formation-related factors (CpEFG1 and CpBCR1) in the BCR1 knockout strain, CP7 and CP5 clinical strains in the presence of NAC. To further explore the mechanism of NAC on the biofilms of C. parapsilosis, RNA sequencing was used to calculate gene expression, comparing the differences among samples. Gene Ontology (GO) enrichment analysis helps to illustrate the difference between two particular samples on functional levels.Results. A high concentration of NAC reduces the total amount of biofilm formation in C. parapsilosis. Following co-incubation with NAC, the expression of CpEFG1 in both CP7 and CP5 clinical strains decreased, while there were no significant changes in the transcriptional levels of CpBCR1 compared with the untreated strain. GO enrichment analysis showed that the metabolism and biofilm-forming factors of cell structure were all regulated after NAC intervention.Conclusions. The non-antifungal agent NAC exhibits certain inhibitory effects on clinical isolate biofilm formation by down-regulating the expression of the CpEFG1 gene, making it a highly potential candidate for the treatment of C. parapsilosis catheter-related bloodstream infections.

O-134 Down regulation of Clock gene in testicle may affect energy metabolism of sertoli cell

Abstract Study question Whether decrease expression of Clock gene in the testis could affect the energy metabolism of Sertoli cells and what is the possible mechanism? Summary answer Down-regulating of Clock gene could regulate the expression of PIK3CB affect AMPK signaling pathways, significantly weakening the energy supply of sertoli cells in mouse testicular. What is known already The causes of Nonobstructive azoospermia (NOA) and other male infertility diseases are diverse and complex, including spermatogenesis and germ cell maturation at different stages of meiosis, the expression level of Clock gene maybe is one reason of which. In the testis of NOA patients, expression level of Clock gene was significantly decreased compared with Obstructive azoospermia (OA), meanwhile, the decrease of Clock gene expression can result in the decrease of male mouse fetal number. The decrease of Clock gene expression cause some kind of male infertility problem is confirmed, but the process and mechanism is not clear. Study design, size, duration Differential expression genes related to low expression of Clock and pathways related to spermatogenesis were screened out through transcriptome analysis on testicular tissue samples from patients with NOA(n = 24) and OA(n = 17). The model was established by transfecting mouse testicular Sertoli cell line (TM4 cells) with lentivirus containing Clock knock-down gene and energy metabolism effects were detected by Seahorse tests. Adifluorescein reporter gene was constructed to detect the binding of CLOCK and candidated differential expressed gene. Participants/materials, setting, methods Normal testis sequencing data were obtained from GTEX database as control. Combined with sequencing data, the differentially expressed genes were obtained as candidate gene, which expression was verified by RT-qPCR, Western-blotting. Energy metabolism were detected by Seahorse mitochondrial-stress test, glycolysis-rate test. Adifluorescein reporter gene for the promoter regions of candidate gene was constructed and transfected into TM4 cells to detect whether CLOCK could regulate candidate gene by binding to the E-box of the promoter region. Main results and the role of chance After GO(GeneOntology) and KEGG(Kyoto Encyclopedia of Genes and Genomes) enrichment analysis, PIK3CB which involved in AMPK signaling pathway was selected for further verification and molecular mechanism study. In Seahorse tests, the basal respiratory capacity, maximal respiratory capacity OCR, spare respiratory capacity and ATP-related respiratory capacity of TM4 cells were significantly decreased after Clock gene knockdown, but there were no significant differences in glycolytic rate, basal glycolysis, compensatory glycolysis level and glycolytic rate. RT-qPCR and Western-blotting confirmed that after knock down of Clock gene, the expression level of Pik3cb gene was down-regulated. That imply when Clock gene is decreased, it regulates the expression of PIK3CB affects AMPK signaling pathways, weakening the ability of mouse testicular sertoli cells to supply energy to various levels of spermatogenic cells and sperm cells, which may be one reason for the occurrence of nonobstructive azoospermia. During this process, there were significant changes in the mitochondrial energy supply of testicular supporting cells, but no significant changes in glycolysis. Finally, the dual luciferin report showed that the Clock gene could bind to the promoter regions of Pik3cb to regulate the expression of PIK3CB, and that the promoter of PIK3CB was activated by the transcription factor CLOCK. Limitations, reasons for caution Although we identified that Clock gene could bind to the promoter regions of Pik3cb to regulate the expression of PIK3CB and affect AMPK signaling pathways, weakening the ability of mouse testicular sertoli cells to supply energy, the available expression level and whether it will be work in vivo is unknown. Wider implications of the findings The results imply that the main energy changes affect the occurrence and development of spermatogonia and primary spermatocytes near the base of the seminiferous tubules are related to the influence of Clock on genes involved in regulating mitochondrial energy supply. It need more further studies. Trial registration number not applicable

Insights into the pathophysiology and response of persistent spinal pain syndrome type 2 to spinal cord stimulation: a human genome-wide association study

BackgroundSpinal cord stimulation (SCS) provides pain relief for some patients with persistent spinal pain syndrome type 2 (PSPS 2), but the precise mechanisms of action and prognostic factors for a...

Weighted Gene Co-expression Network Analysis (WGCNA) of Wnt Signaling Related to Periodontal Ligament Formation: A Bioinformatics-Based Analysis.

Introduction The Wnt signaling pathway is crucial for tooth development, odontoblast differentiation, and dentin formation. It interacts with epithelial cadherin (E-cadherin) and beta-catenin in tooth development and periodontal ligament (PDL) formation. Dysregulation of Wnt signaling is linked to periodontal diseases, requiring an understanding of therapeutic interventions. Weighted gene co-expression network analysis (WGCNA) can identify co-expressed gene modules. Our study aims to identify hub genes in WGCNA analysis of Wnt signaling-based PDL formation. Methods The study used a microarray datasetGSE201313 from the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus to analyze the impact of DMP1 expression on XLH dental pulp cell differentiation and PDL formation. The standardized dataset was used for WGCNA analysis, which generated a co-expression network by calculating pairwise correlations between genes and constructing an adjacency matrix. The topological overlap matrix (TOM) was transformed into a hierarchical clustering tree and then cut into modules or clusters of highly interconnected genes. The module eigengene (ME) was calculated for each module, and the genes within this module were identified as hub genes. Gene ontology (GO) and KEGG pathway enrichment analysis were performed to gain insights into the biological functions of the hub genes. The integrated Differential Expression and Pathway analysis (iDEP) tool (http://bioinformatics.sdstate.edu/idep/; South Dakota State University, Brookings, USA) was used for WGCNA analysis. Results The study used the WGCNA package to analyze 1,000 differentially expressed genes, constructing a gene co-expression network and generating a hierarchical clustering tree and TOM. The analysis reveals a scale-free topology fitting index R2 and mean connectivity for various soft threshold powers, with an R2 value of 5. COL6A1, MMP3, BGN, COL1A2, and FBN2 are hub genes implicated in PDL development. Conclusion The study identified key hub genes, including COL6A1, MMP3, BGN, and FBN2, crucial for PDL formation, tissue remodeling, and cell-matrix interactions, guiding future therapeutic strategies.

Rapid Mining of Candidate Genes for the Branchless Phenotype in Watermelon by Bulked Segregant Analysis Using Whole-genome Resequencing

Lateral branching is an important agronomic trait in horticulture plants. The aim of this project was to reveal the genetic mechanism, map the gene localization, and predict candidate genes of watermelon lateral branching. An F2 segregating population was derived from a cross between the multibranched maternal inbred line M6 and the branchless paternal inbred line N7. Two DNA pools were constructed using 20 multibranched plants and 20 branchless plants from the F2 population. Whole-genome resequencing was performed for the DNA pools (25×) and the parents (30×) to identify the genomic region associated with lateral branching. Candidate genes were predicted based on the gene annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses; then, quantitative validation of these genes was performed. The results showed that the clean reads of four samples yielded 64,295,076 to 81,658,958 bp, with sufficient genome coverage and high quality. Based on single-nucleotide polymorphism and insertions/deletions association analyses, the candidate genes were mapped to a 2.01-Mb region on chromosome 4 (22,958,925–24,971,894 bp) containing 182 annotated genes. During the KEGG and GO enrichment analyses, these genes were annotated to 10 cellular components, 10 molecular functions, and 12 biological processes. Eight candidate genes responsible for the branchless phenotype in watermelon were identified: Cla97C04G076340, Cla97C04G075820, Cla97C04G076060, Cla97C04G076250, Cla97C04G076280, Cla97C04G076380, Cla97C04G076830, and Cla97C04G075950. These genes were involved in amino acid biosynthesis and catabolism, TCP transcription factor activity, and regulation of flower development. This study offers valuable insights into the molecular mechanisms governing the branchless phenotype in watermelon. These candidate genes serve as potential targets for gene cloning and marker-assisted selection of watermelon cultivars without lateral branches.

Characteristics, clinical significance, and cancer immune interactions of lipid metabolism in prostate cancer.

The relationship between lipid metabolism, immune response, and immunotherapy in prostate cancer (PCa) is closely intertwined, and targeted intervention in lipid metabolism may facilitate the success of anticancer immunotherapy. This research attempted to explore effective immunotherapy for PCa. We obtained RNA sequencing (RNA-seq) data for PCa patients from the UCSC Xena platform. Data analysis of differentially expressed genes (DEGs) was performed using package limma in R. Then, DEGs were subjected to enrichment analysis of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The Human Protein Atlas (HPA) database was conducted to validate the protein expression of the up-regulated lipid metabolism related genes (LMRGs) between PCa tissues and normal prostate tissues. And then we identified critical transcription factors (TFs), LMRGs and miRNA by constructing a regulatory network of TF-gene-miRNA. Furthermore, we determined the high and low groups based on the score of lipid metabolism enrichment. The hallmark gene sets were derived from gene expression profiles using the gene set variation analysis (GSVA) R package. Finally, we conducted immune infiltration analysis and drug sensitivity analysis. Immune response and lipid metabolism have undergone significant changes in PCa and paracancerous tissues compared to normal tissues. A total of 21 LMRGs were differentially up-regulated in PCa. The TF-gene-miRNA network showed that PLA2G7, TWIST1, and TRIB3 may be the key genes that elevated lipid metabolism in PCa. The high group had more infiltration of B cell memory, macrophage M0, macrophage M1, and myeloid dendritic cell resting, and the low group had more infiltration of B cell plasma, monocyte, myeloid dendritic cell activated, and mast cell resting. The majority of checkpoint genes exhibited high expression levels in the low group. Lipid metabolism was remarkedly correlated with drug sensitivity. The analysis of lipid metabolism and related genes has revealed a complex regulatory mechanism that has a significant influence on immune response, immunotherapy, and medication guidance for patients with PCa. Prostate cancer (PCa); lipid metabolism; cancer immune; RNA sequencing (RNA-seq).

Identification of hub genes and key modules in laryngeal squamous cell carcinoma.

Laryngeal squamous cell carcinoma (LSCC) is the prominent cancer in head and neck, which greatly affects life quality of patients. The pathogenesis of LSCC is not clear. Presently, the LSCC treatments include chemotherapy, surgery and radiotherapy; however, these methods have poor efficacy in patients with recurrent and persistent cancer. Therefore, the study identified the hub genes accompanied with LSCC, which may be a potential therapeutic target in the future. We extracted whole transcriptome high-throughput sequencing (HTS) LSCC data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases and calculate differentially expressed genes (DEGs) between LSCC and normal samples using statistical software RStudio. Through weighted gene co-expression network analysis (WGCNA), enrichment examination of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and Gene Ontology (GO) functions, and examination of protein-protein interaction (PPI) network, we obtained network hub genes and validated the hub genes prognostic value and expression levels of protein. Through analysis of differential gene expression, from the GEO and TCGA databases 2,139 and 2,774 DEGs were obtained, respectively, 13 and 15 modules were screened from TCGA-LSCC and GSE127165 datasets by WGCNA, respectively. The most significant positive and negative correlation modules in the WGCNA and DEG lists were overlapped, and overall 36 co-expressed overlapping genes were retrieved. Through enrichment analysis of GO and KEGG, it was found that the gene functions were highly concentrated in cell junction assembly, basement membrane, extracellular matrix (ECM) structural constituent etc., and the pathways were mainly concentrated in ECM receptor interaction, focal adhesion, small cell lung cancer, and toxoplasmosis. Through analysis of PPI network analysis, 10 network hub genes (SNAI2, ITGA6, LAMB3, LAMC2, CAV1, COL7A1, GJA1, EHF, OAT, and GPT) were obtained. Finally, survival analysis and protein expression validation of these genes confirmed that low OAT expression and high CAV1 expression remarkably influenced the survival of patient's prognosis with LSCC. We recognized the hub genes and key modules nearly associated to LSCC and these genes were validated by survival analysis and the database of Human Protein Atlas (HPA), which is of high importance for unveiling the pathogenesis of LSCC and probing for new precise biological marker and potential therapeutic targets.

Integrated approach of network pharmacology, molecular docking, and clinical observations in evaluating the efficacy and safety of Bufei Huoxue capsules for pulmonary hypertension associated with chronic obstructive pulmonary disease.

Chronic obstructive pulmonary disease (COPD) is a persistent and progressive disorder characterized by airway or alveolar abnormalities, commonly leading to pulmonary hypertension (PH). This clinical observational study investigates the therapeutic mechanisms of Bufei Huoxue capsules (BHC) in treating PH in patients with COPD-linked PH (COPD-PH) using network pharmacology and molecular docking methods, and assesses the therapeutic efficacy and safety of BHCs. The active compounds and their target proteins in BHCs were sourced from the Traditional Chinese Medicine Systems Pharmacology database, with additional target proteins derived from the GeneCards and OMIM databases. An active network was constructed using Cytoscape 3.7.1, and interaction networks were established. Intersecting targets underwent Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis using the Metascape database. Network pharmacology and molecular docking studies demonstrated favorable binding affinities of BHC active ingredients, such as quercetin, bavachalcone, and isobavachin, for key targets including PTGS1, ESR1, and PTGS2. Gene Ontology enrichment analysis highlighted the involvement of these targets in processes such as the positive regulation of locomotion, the transmembrane receptor protein tyrosine kinase signaling pathway, and peptidyl-tyrosine phosphorylation. KEGG pathway analysis indicated their roles in pathways related to cancer, AGE-RAGE signaling in diabetic complications, and prostate cancer. BHCs exhibit therapeutic effects on COPD-PH through multi-component, multi-target, and multi-pathway interactions. This clinical observational study confirms the efficacy and safety of BHCs in improving cardiac and pulmonary functions, enhancing exercise tolerance, and elevating the quality of life in patients with COPD-PH.

Identification of immune-associated biomarkers of diabetes nephropathy tubulointerstitial injury based on machine learning: a bioinformatics multi-chip integrated analysis

BackgroundDiabetic nephropathy (DN) is a major microvascular complication of diabetes and has become the leading cause of end-stage renal disease worldwide. A considerable number of DN patients have experienced irreversible end-stage renal disease progression due to the inability to diagnose the disease early. Therefore, reliable biomarkers that are helpful for early diagnosis and treatment are identified. The migration of immune cells to the kidney is considered to be a key step in the progression of DN-related vascular injury. Therefore, finding markers in this process may be more helpful for the early diagnosis and progression prediction of DN.MethodsThe gene chip data were retrieved from the GEO database using the search term ' diabetic nephropathy ‘. The ' limma ' software package was used to identify differentially expressed genes (DEGs) between DN and control samples. Gene set enrichment analysis (GSEA) was performed on genes obtained from the molecular characteristic database (MSigDB. The R package ‘WGCNA’ was used to identify gene modules associated with tubulointerstitial injury in DN, and it was crossed with immune-related DEGs to identify target genes. Gene ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed on differentially expressed genes using the ‘ClusterProfiler’ software package in R. Three methods, least absolute shrinkage and selection operator (LASSO), support vector machine recursive feature elimination (SVM-RFE) and random forest (RF), were used to select immune-related biomarkers for diagnosis. We retrieved the tubulointerstitial dataset from the Nephroseq database to construct an external validation dataset. Unsupervised clustering analysis of the expression levels of immune-related biomarkers was performed using the ‘ConsensusClusterPlus ‘R software package. The urine of patients who visited Dongzhimen Hospital of Beijing University of Chinese Medicine from September 2021 to March 2023 was collected, and Elisa was used to detect the mRNA expression level of immune-related biomarkers in urine. Pearson correlation analysis was used to detect the effect of immune-related biomarker expression on renal function in DN patients.ResultsFour microarray datasets from the GEO database are included in the analysis : GSE30122, GSE47185, GSE99340 and GSE104954. These datasets included 63 DN patients and 55 healthy controls. A total of 9415 genes were detected in the data set. We found 153 differentially expressed immune-related genes, of which 112 genes were up-regulated, 41 genes were down-regulated, and 119 overlapping genes were identified. GO analysis showed that they were involved in various biological processes including leukocyte-mediated immunity. KEGG analysis showed that these target genes were mainly involved in the formation of phagosomes in Staphylococcus aureus infection. Among these 119 overlapping genes, machine learning results identified AGR2, CCR2, CEBPD, CISH, CX3CR1, DEFB1 and FSTL1 as potential tubulointerstitial immune-related biomarkers. External validation suggested that the above markers showed diagnostic efficacy in distinguishing DN patients from healthy controls. Clinical studies have shown that the expression of AGR2, CX3CR1 and FSTL1 in urine samples of DN patients is negatively correlated with GFR, the expression of CX3CR1 and FSTL1 in urine samples of DN is positively correlated with serum creatinine, while the expression of DEFB1 in urine samples of DN is negatively correlated with serum creatinine. In addition, the expression of CX3CR1 in DN urine samples was positively correlated with proteinuria, while the expression of DEFB1 in DN urine samples was negatively correlated with proteinuria. Finally, according to the level of proteinuria, DN patients were divided into nephrotic proteinuria group (n = 24) and subrenal proteinuria group. There were significant differences in urinary AGR2, CCR2 and DEFB1 between the two groups by unpaired t test (P < 0.05).ConclusionsOur study provides new insights into the role of immune-related biomarkers in DN tubulointerstitial injury and provides potential targets for early diagnosis and treatment of DN patients. Seven different genes ( AGR2, CCR2, CEBPD, CISH, CX3CR1, DEFB1, FSTL1 ), as promising sensitive biomarkers, may affect the progression of DN by regulating immune inflammatory response. However, further comprehensive studies are needed to fully understand their exact molecular mechanisms and functional pathways in DN.

Genome-wide methylation landscape during somatic embryogenesis in Medicago truncatula reveals correlation between Tnt1 retrotransposition and hyperactive methylation regions.

Medicago truncatula is a model legume for fundamental research on legume biology and symbiotic nitrogen fixation. Tnt1, a retrotransposon from tobacco, was used to generate insertion mutants in M. truncatula R108. Approximately 21 000 insertion lines have been generated and publicly available. Tnt1 retro-transposition event occurs during somatic embryogenesis (SE), a pivotal process that triggers massive methylation changes. We studied the SE of M. truncatula R108 using leaf explants and explored the dynamic shifts in the methylation landscape from leaf explants to callus formation and finally embryogenesis. Higher cytosine methylation in all three contexts of CG, CHG, and CHH patterns was observed during SE compared to the controls. Higher methylation patterns were observed in assumed promoter regions (~2-kb upstream regions of transcription start site) of the genes, while lowest was recorded in the untranslated regions. Differentially methylated promoter region analysis showed a higher CHH methylation in embryogenesis tissue samples when compared to CG and CHG methylation. Strong correlation (89.71%) was identified between the differentially methylated regions (DMRs) and the site of Tnt1 insertions in M.truncatula R108 and stronger hypermethylation of genes correlated with higher number of Tnt1 insertions in all contexts of CG, CHG, and CHH methylation. Gene ontology enrichment and KEGG pathway enrichment analysis identified genes and pathways enriched in the signal peptide processing, ATP hydrolysis, RNA polymerase activity, transport, secondary metabolites, and nitrogen metabolism pathways. Combined gene expression analysis and methylation profiling showed an inverse relationship between methylation in the DMRs (regions spanning genes) and the expression of genes. Our results show that a dynamic shift in methylation happens during the SE process in the context of CG, CHH and CHG methylation, and the Tnt1 retrotransposition correlates with the hyperactive methylation regions.

Effects of Kojic Acid-mediated Sonodynamic Therapy as a Matrix Metalloprotease-9 Inhibitor against Oral Squamous Cell Carcinoma: A Bioinformatics Screening and In Vitro Analysis.

Oral squamous cell carcinoma (OSCC) is a type of cancer that is responsible for a significant amount of morbidity and mortality worldwide. Researchers are searching for promising therapeutic methods to manage this cancer. In this study, an in silico approach was used to evaluate the activity of sonodynamic therapy (SDT) based on the use of Kojic acid as a sonosensitizer to inhibit matrix metalloprotease-9 (MMP-9) in OSCC. The three-dimensional structure of MMP-9 was predicted and validated by computational approaches. The possible functional role of MMP-9 was determined in terms of Gene Ontology (GO) enrichment analysis. In silico, molecular docking was then performed to evaluate the binding energies of Kojic acid with MMP-9, and ADME parameters and toxicity risks were predicted. The pharmacokinetics and drug-likeness properties of Kojic acid were assessed. Moreover, after the determination of the cytotoxicity effect of Kojic acid-mediated SDT, the change of mmp-9 gene expression was assessed on OSCC cells. The results of the study showed that Kojic acid could efficiently interact with MMP-9 protein with a strong binding affinity. Kojic acid obeyed Lipinski's rule of five without violation and exhibited drug-likeness. The cytotoxic effects of Kojic acid and ultrasound waves on the OSCC cells were dose-dependent, and the lowest expression level of the mmp-9 gene was observed in SDT. Overall, Kojic acid-mediated SDT as an MMP-9 inhibitor can be a promising adjuvant treatment for OSCC. The study highlights the potential of In silico approaches to evaluate therapeutic methods for cancer treatment.

An improved bacterial mRNA enrichment strategy in dual RNA sequencing to unveil the dynamics of plant-bacterial interactions

BackgroundDual RNA sequencing is a powerful tool that enables a comprehensive understanding of the molecular dynamics underlying plant-microbe interactions. RNA sequencing (RNA-seq) poses technical hurdles in the transcriptional analysis of plant-bacterial interactions, especially in bacterial transcriptomics, owing to the presence of abundant ribosomal RNA (rRNA), which potentially limits the coverage of essential transcripts. Therefore, to achieve cost-effective and comprehensive sequencing of the bacterial transcriptome, it is imperative to devise efficient methods for eliminating rRNA and enhancing the proportion of bacterial mRNA. In this study, we modified a strand-specific dual RNA-seq method with the goal of enriching the proportion of bacterial mRNA in the bacteria-infected plant samples. The enriched method involved the sequential separation of plant mRNA by poly A selection and rRNA removal for bacterial mRNA enrichment followed by strand specific RNA-seq library preparation steps. We assessed the efficiency of the enriched method in comparison to the conventional method by employing various plant-bacterial interactions, including both host and non-host resistance interactions with pathogenic bacteria, as well as an interaction with a beneficial rhizosphere associated bacteria using pepper and tomato plants respectively.ResultsIn all cases of plant-bacterial interactions examined, an increase in mapping efficiency was observed with the enriched method although it produced a lower read count. Especially in the compatible interaction with Xanthmonas campestris pv. Vesicatoria race 3 (Xcv3), the enriched method enhanced the mapping ratio of Xcv3-infected pepper samples to its own genome (15.09%; 1.45-fold increase) and the CDS (8.92%; 1.49-fold increase). The enriched method consistently displayed a greater number of differentially expressed genes (DEGs) than the conventional RNA-seq method at all fold change threshold levels investigated, notably during the early stages of Xcv3 infection in peppers. The Gene Ontology (GO) enrichment analysis revealed that the DEGs were predominantly enriched in proteolysis, kinase, serine type endopeptidase and heme binding activities.ConclusionThe enriched method demonstrated in this study will serve as a suitable alternative to the existing RNA-seq method to enrich bacterial mRNA and provide novel insights into the intricate transcriptomic alterations within the plant-bacterial interplay.

Identification and validation of endoplasmic reticulum stress-related genes that enhance immunotherapy in colon cancer.

Endoplasmic reticulum stress (ERS)-related genes are related to tumor growth, metastasis, and immunotherapy response. In this paper, we tried to identify ERS-related genes related to immunotherapy in colon cancer. ERS-related genes were downloaded from the Molecular Signatures Database (MSigDB) and GeneCards websites. Normal and tumor samples of the colon were obtained from The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression Project (GTEx), and Gene Expression Omnibus (GEO) databases. A risk model based on gene coefficients was constructed by using the least absolute shrinkage and selection operator (LASSO) regression. The inherent biological process differences between risk groups were explored by Gene Ontology (GO) and gene set enrichment analysis (GSEA). ESTIMATE and single-sample GSEA (ssGSEA) algorithms were used to analyze the correlation between tumor microenvironment (TME) and immune checkpoint and risk score. The semi-inhibitory concentration (IC50) values of chemotherapeutic drugs between risk groups were calculated to evaluate the sensitivity of immunotherapy. The pathway analysis showed that the ERS risk model was relevant to biosynthesis and metabolism. Consistent clustering based on the ERS-related differentially expressed genes (DEGs) demonstrated that the samples divided into three clusters had significant clinicopathological differences. A risk model consisting of six ERS-related genes was established. The model was verified on GSE39582 and GSE17536 testing datasets. The results showed that ERS risk model was significantly related to TME and immune checkpoint, and these genes enhanced the immunotherapy ability of colon cancer. We established a risk model with ERS-related genes (PMM2, STC2, EIF2AK1, HSPA1A, SLC8A1, KCNQ1), which enhance the sensitivity of immunotherapy for colon cancer. These may provide a new perspective for the treatment of colon cancer.

A comprehensive pan-cancer analysis revealing SPAG6 as a novel diagnostic, prognostic and immunological biomarker in tumor.

There have been studies on the role of sperm-associated antigen 6 (SPAG6) in cytoskeleton formation and growth cone stability, but it is also unknown how spag6 affect tumor growth and development. The aim of this study was to clarify the role of SPAG6 in pan-cancer, with some findings about thyroid carcinoma (THCA) validated through experiments. We examined the role of SPAG6 in pan-cancer, with the data being collected from databases. Further analysis was conducted to assess its correlations with prognosis, gene heterogeneity, stemness, and tumor immunity. The interacting proteins of SPAG6 were also identified, and gene ontology enrichment analysis was performed to determine its biological function. We preliminarily confirmed the role of SPAG6 via in vitro experiments and immunofluorescence staining. This study found that SPAG6 expression was differentially expressed in cancers and at various tumor stages and grades. In stomach and esophageal carcinoma (STES), stomach adenocarcinoma (STAD), kidney renal clear cell carcinoma (KIRC), lung squamous cell carcinoma (LUSC), and adrenocortical carcinoma (ACC), SPAG6 expression was correlated with gender. SPAG6 expression was also found to be correlated with prognostic value, with low expression being associated with poor prognosis. Furthermore, SPAG6 expression was positively linked with immune-related cells in HNSC, chemokine receptors in LUSC, and immune checkpoint genes in THCA. Furthermore, SPAG6 overexpression suppressed the malignant phenotypes of THCA cells, manifested by slower proliferation and decreased migration. The different SPAG6 expression in THCA led to different malignant phenotypes, which are involved in the upregulation of DNA repair, MYC targets, peroxisome, and G2M checkpoint. SPAG6 plays a significant role as an oncogene and can be used as a marker to predict the prognosis of cancer. SPAG6 influences both the tumor immune infiltration and microenvironment, making it a promising immunotherapeutic target for tumor therapy.

Identification of Potential New Genes Related to the SREBP Pathway in Xanthophyllomyces dendrorhous.

The sterol regulatory element-binding protein (SREBP) pathway is an integral cellular mechanism that regulates lipid homeostasis, in which transcriptional activator SREBPs regulate the expression of various genes. In the carotenogenic yeast Xanthophyllomyces dendrorhous, Sre1 (the yeast SREBP homolog) regulates lipid biosynthesis and carotenogenesis, among other processes. Despite the characterization of several components of the SREBP pathway across various eukaryotes, the specific elements of this pathway in X. dendrorhous remain largely unknown. This study aimed to explore the potential regulatory mechanisms of the SREBP pathway in X. dendrorhous using the strain CBS.cyp61- as a model, which is known to have Sre1 in its active state under standard culture conditions, resulting in a carotenoid-overproducing phenotype. This strain was subjected to random mutagenesis with N-methyl-N'-nitro-N-nitrosoguanidine (NTG), followed by a screening methodology that focused on identifying mutants with altered Sre1 activation phenotypes. Single-nucleotide polymorphism (SNP) analysis of 20 selected mutants detected 5439 single-nucleotide variants (SNVs), narrowing them down to 1327 SNPs of interest after a series of filters. Classification based on SNP impact identified 116 candidate genes, including 49 genes with high impact and 68 genes with deleterious moderate-impact mutations. BLAST, InterProScan, and gene ontology enrichment analyses highlighted 25 genes as potential participants in regulating Sre1 in X. dendrorhous. The key findings of this study include the identification of genes potentially encoding proteins involved in protein import/export to the nucleus, sterol biosynthesis, the ubiquitin-proteasome system, protein regulatory activities such as deacetylases, a subset of kinases and proteases, as well as transcription factors that could be influential in SREBP regulation. These findings are expected to significantly contribute to the current understanding of the intricate regulation of the transcription factor Sre1 in X. dendrorhous, providing valuable groundwork for future research and potential biotechnological applications.

Integrated Analysis of Transcriptome Profiles and lncRNA-miRNA-mRNA Competing Endogenous RNA Regulatory Network to Identify Biological Functional Effects of Genes and Pathways Associated with Johne's Disease in Dairy Cattle.

Paratuberculosis or Johne's disease (JD), a chronic granulomatous gastroenteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP), causes huge economic losses and reduces animal welfare in dairy cattle herds worldwide. At present, molecular mechanisms and biological functions involved in immune responses to MAP infection of dairy cattle are not clearly understood. Our purpose was to integrate transcriptomic profiles and competing endogenous RNA (ceRNA) network analyses to identify key messenger RNAs (mRNAs) and regulatory RNAs involved in molecular regulation of peripheral blood mononuclear cells (PBMCs) for MAP infection in dairy cattle. In total, 28 lncRNAs, 42 miRNAs, and 370 mRNAs were identified by integrating gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. In this regard, we identified 21 hub genes (CCL20, CCL5, CD40, CSF2, CXCL8, EIF2AK2, FOS, IL10, IL17A, IL1A, IL1B, IRF1, MX2, NFKB1, NFKBIA, PTGS2, SOCS3, TLR4, TNF, TNFAIP3, and VCAM1) involved in MAP infection. Furthermore, eight candidate subnets with eight lncRNAs, 29 miRNAs, and 237 mRNAs were detected through clustering analyses, whereas GO enrichment analysis of identified RNAs revealed 510, 22, and 11 significantly enriched GO terms related to MAP infection in biological process, molecular function, and cellular component categories, respectively. The main metabolic-signaling pathways related to MAP infection that were enriched included the immune system process, defense response, response to cytokine, leukocyte migration, regulation of T cell activation, defense response to bacterium, NOD-like receptor, B cell receptor, TNF, NF-kappa B, IL-17, and T cell receptor signaling pathways. Contributions of transcriptome profiles from MAP-positive and MAP-negative sample groups plus a ceRNA regulatory network underlying phenotypic differences in the intensity of pathogenicity of JD provided novel insights into molecular mechanisms associated with immune system responses to MAP infection in dairy cattle.

Network Pharmacology Based Elucidation of Molecular Mechanisms of Laoke Formula for Treatment of Advanced Non-Small Cell Lung Cancer.

To explore the specific pharmacological molecular mechanisms of Laoke Formula (LK) on treating advanced non-small cell lung cancer (NSCLC) based on clinical application, network pharmacology and experimental validation. Kaplan-Meier method and Cox regression analysis were used to evaluate the survival benefit of Chinese medicine (CM) treatment in 296 patients with NSCLC in Tianjin Medical University Cancer Institute and Hospital from January 2011 to December 2015. The compounds of LK were screened using the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform, and the corresponding targets were performed from Swiss Target Prediction. NSCLC-related targets were obtained from Therapeutic Target Database and Comparative Toxicogenomics Database. Key compounds and targets were identified from the compound-target-disease network and protein-protein interaction (PPI) network analysis, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analysis were used to predict the potential signaling pathways involved in the treatment of advanced NSCLC with LK. The binding affinities between key ingredients and targets were further verified using molecular docking. Finally, A549 cell proliferation and migration assay were used to evaluate the antitumor activity of LK. Western blot was used to further verify the expression of key target proteins related to the predicted pathways. Kaplan-Meier survival analysis showed that the overall survival of the CM group was longer than that of the non-CM group (36 months vs. 26 months), and COX regression analysis showed that LK treatment was an independent favorable prognostic factor (P=0.027). Next, 97 components and 86 potential targets were included in the network pharmacology, KEGG and GO analyses, and the results indicated that LK was associated with proliferation and apoptosis. Moreover, molecular docking revealed a good binding affinity between the key ingredients and targets. In vitro, A549 cell proliferation and migration assay showed that the biological inhibition effect was more obvious with the increase of LK concentration (P<0.05). And decreased expressions of nuclear factor κB1 (NF-κB1), epidermal growth factor receptor (EGFR) and AKT serine/threonine kinase 1 (AKT1) and increased expression of p53 (P<0.05) indicated the inhibitory effect of LK on NSCLC by Western blot. LK inhibits NSCLC by inhibiting EGFR/phosphoinositide 3-kinase (PI3K)/AKT signaling pathway, NFκB signaling pathway and inducing apoptosis, which provides evidence for the therapeutic mechanism of LK to increase overall survival in NSCLC patients.