Functional Pathway Enrichment Research Articles

The urinary microbiota composition and functionality of calcium oxalate stone formers.

Accumulated evidences indicate that dysbiosis of the urinary microbiota is associated with kidney stone formation. In the present study, we aimed to investigate the urinary microbiota composition and functionality of patients with calcium oxalate stones and compare it with those of healthy individuals. We collected bladder urine samples from 68 adult patients with calcium oxalate stones and 54 age-matched healthy controls by transurethral catheterization. 16S rRNA gene and shotgun sequencing were utilized to characterize the urinary microbiota and functionality associated with calcium oxalate stones. After further exclusion, a total of 100 subjects was finally included and analyzed. The diversity of the urinary microbiota in calcium oxalate stone patients was not significantly different from that of healthy controls. However, the urinary microbiota structure of calcium oxalate stone formers significantly differed from that of healthy controls (PERMANOVA, r = 0.026, P = 0.019). Differential representation of bacteria (e.g., Bifidobacterium) and several enriched functional pathways (e.g., threonine biosynthesis) were identified in the urine of calcium oxalate stone patients. Our results showed significantly different urinary microbiota structure and several enriched functional pathways in calcium oxalate stone patients, which provide new insight into the pathogenesis of calcium oxalate stones.

Advances in Inflammatory Bowel Disease Diagnostics: Machine Learning and Genomic Profiling Reveal Key Biomarkers for Early Detection.

This study, utilizing high-throughput technologies and Machine Learning (ML), has identified gene biomarkers and molecular signatures in Inflammatory Bowel Disease (IBD). We could identify significant upregulated or downregulated genes in IBD patients by comparing gene expression levels in colonic specimens from 172 IBD patients and 22 healthy individuals using the GSE75214 microarray dataset. Our ML techniques and feature selection methods revealed six Differentially Expressed Gene (DEG) biomarkers (VWF, IL1RL1, DENND2B, MMP14, NAAA, and PANK1) with strong diagnostic potential for IBD. The Random Forest (RF) model demonstrated exceptional performance, with accuracy, F1-score, and AUC values exceeding 0.98. Our findings were rigorously validated with independent datasets (GSE36807 and GSE10616), further bolstering their credibility and showing favorable performance metrics (accuracy: 0.841, F1-score: 0.734, AUC: 0.887). Our functional annotation and pathway enrichment analysis provided insights into crucial pathways associated with these dysregulated genes. DENND2B and PANK1 were identified as novel IBD biomarkers, advancing our understanding of the disease. The validation in independent cohorts enhances the reliability of these findings and underscores their potential for early detection and personalized treatment of IBD. Further exploration of these genes is necessary to fully comprehend their roles in IBD pathogenesis and develop improved diagnostic tools and therapies. This study significantly contributes to IBD research with valuable insights, potentially greatly enhancing patient care.

Tissue Accumulation and Biotransformation of 6PPD-Quinone in Adult Zebrafish and Its Effects on the Intestinal Microbial Community.

The pronounced lethality of N-(1,3-dimethylbutyl)-N'-phenyl-p-phenylenediamine quinone (6PPD-quinone or 6PPDQ) toward specific salmonids, while sparing other fish species, has received considerable attention. However, the underlying cause of this species-specific toxicity remains unresolved. This study explored 6PPDQ toxicokinetics and intestinal microbiota composition in adult zebrafish during a 14-day exposure to environmentally realistic concentrations, followed by a 7-day recovery phase. Predominant accumulation occurred in the brain, intestine, and eyes, with the lowest levels in the liver. Six metabolites were found to undergo hydroxylation, with two additionally undergoing O-sulfonation. Semiquantitative analyses revealed that the predominant metabolite featured a hydroxy group situated on the phenyl ring adjacent to the quinone. This was further validated by assessing enzyme activity and determining in silico binding interactions. Notably, the binding affinity between 6PPDQ and zebrafish phase I and II enzymes exceeded that with the corresponding coho salmon enzymes by 1.04-1.53 times, suggesting a higher potential for 6PPDQ detoxification in tolerant species. Whole-genome sequencing revealed significant increases in the genera Nocardioides and Rhodococcus after exposure to 6PPDQ. Functional annotation and pathway enrichment analyses predicted that these two genera would be responsible for the biodegradation and metabolism of xenobiotics. These findings offer crucial data for comprehending 6PPDQ-induced species-specific toxicity.

Oleanolic Acid Promotes the Formation of Probiotic Escherichia coli Nissle 1917 (EcN) Biofilm by Inhibiting Bacterial Motility.

Probiotic biofilms have been beneficial in the fight against infections, restoring the equilibrium of the host's gut microbiota, and enhancing host health. They are considered a novel strategy for probiotic gut colonization. In this case, we evaluated the effects of various active substances from traditional Chinese medicine on Escherichia coli Nissle 1917 (EcN) to determine if they promote biofilm formation. It was shown that 8-64 μg/mL of oleanolic acid increased the development of EcN biofilm. Additionally, we observed that oleanolic acid can effectively suppress biofilm formation in pathogenic bacteria such as Salmonella and Staphylococcus aureus. Next, we assessed the amount of EcN extracellular polysaccharides, the number of live bacteria, their metabolic activity, the hydrophobicity of their surface, and the shape of their biofilms using laser confocal microscopy. Through transcriptome analysis, a total of 349 differentially expressed genes were identified, comprising 134 upregulated and 215 downregulated genes. GO functional enrichment analysis and KEGG pathway enrichment analysis revealed that oleanolic acid functions are through the regulation of bacterial motility, the iron absorption system, the two-component system, and adhesion pathways. These findings suggest that the main effects of oleanolic acid are to prevent bacterial motility, increase initial adhesion, and encourage the development of EcN biofilms. In addition, oleanolic acid interacts with iron absorption to cooperatively control the production of EcN biofilms within an optimal concentration range. Taking these results together, this study suggests that oleanolic acid may enhance probiotic biofilm formation in the intestines, presenting new avenues for probiotic product development.

Transcriptomic Analysis of Newborn Hanwoo Calves: Effects of Maternal Overnutrition during Mid- to Late Pregnancy on Subcutaneous Adipose Tissue and Liver.

This study investigated the transcriptomic responses of subcutaneous adipose tissue (SAT) and liver in newborn Hanwoo calves subjected to maternal overnutrition during mid- to late gestation. Eight Hanwoo cows were randomly assigned to control and treatment groups. The treatment group received a diet of 4.5 kg of concentrate and 6.5 kg of rice straw daily, resulting in intake levels of 8.42 kg DMI, 5.69 kg TDN, and 0.93 kg CP-higher than the control group (6.07 kg DMI, 4.07 kg TDN, and 0.65 kg CP), with respective NEm values of 9.56 Mcal and 6.68 Mcal. Following birth, newly born calves were euthanized humanely as per ethical guidelines, and SAT and liver samples from newborn calves were collected for RNA extraction and analysis. RNA sequencing identified 192 genes that were differentially expressed in the SAT (17 downregulated and 175 upregulated); notably, HSPA6 emerged as the most significantly upregulated gene in the SAT and as the singular upregulated gene in the liver (adj-p value < 0.05). Additionally, differential gene expression analysis highlighted extensive changes across genes associated with adipogenesis, fibrogenesis, and stress response. The functional enrichment pathway and protein-protein interaction (PPI) unraveled the intricate networks and biological processes impacted by overnutrition, including extracellular matrix organization, cell surface receptor signaling, and the PI3K-Akt signaling pathway. These findings underscore maternal overnutrition's substantial influence on developmental pathways, suggesting profound cellular modifications with potential lasting effects on health and productivity. Despite the robust insights that are provided, the study's limitations (sample size) underscore the necessity for further research.

Age-Related Clusters and Favorable Immune Phenotypes in Young Breast Cancer Patients

Breast cancer (BC) patients aged <40 years at diagnosis experience aggressive disease and poorer survival compared with women diagnosed with BC at 40 to 49 years, but the age-related biology is described to little extent. Here, we explored transcriptional alterations in BC to gain better understanding of age-related tumor biology. We studied a subset of the Bergen in-house cohort (n = 127; age range, 26-49 years) and used the NanoString Breast Cancer 360 expression panel on formalin-fixed paraffin-embedded BC tissue, and publicly available global BC messenger RNA expression data (n = 204, age range, 22-49 years), to explore differentially expressed genes between the young (age <40 years) and older (age 40-49 years) patients. Unsupervised hierarchical clustering was applied to identify gene expression–based patient clusters. We applied established computational approaches to define the PAM50 subtypes, risk of recurrence scores (ROR), and risk groups and to infer the proportions of 22 immune cell types from bulk gene expression profiles of patients aged <50 years at BC diagnosis. Differentially expressed genes and gene sets were investigated using OncoEnrichR and g:Profiler to describe functional profiles and pathway enrichment. We identified 4 age-related patient clusters presenting distinct characteristics of PAM50 subtypes and ROR profiles, which demonstrated independent prognostic value when adjusted for traditional clinicopathologic variables and the known molecular subtypes. Our findings showed better survival than expected in the basal-enriched cluster 2 and in triple-negative and basal-like BC. Deconvolution analyses of immunophenotypes indicated higher levels of M0 and M1 macrophages than M2 macrophages in subsets of young BC. Our approach identifies age-based patient clusters with distinct clinicopathologic profiles, to a large extent overlapping with the PAM50 subtypes, although with independent prognostic values in multivariate survival analyses. The patient clusters provided new insight in the immune cell distribution across tumor subtypes, potentially contributing to survival differences between the clusters and the molecular subtypes and indicating age-related mechanisms improving outcome. Our study confirms the applicability of ROR as a valid prognosticator also in a young BC cohort.

Identification of potential biomarkers for pancreatic ductal adenocarcinoma: a bioinformatics analysis.

PDA is an aggressive cancer with a 5-year survival rate, which is very low. There is no effective prognosis or therapy for PDA because of the lack of target biomarkers. The objective of this article is to identify the target biomarkers for PDA using a bioinformatics approach. In this work, we have analysed the three microarray datasets from the NCBI GEO database. We used the Geo2R tool to analyse the microarray data with the Benjamini and Hochberg false discovery rate method, and the significance level cut-off was set to 0.05. We have identified 659 DEGs from the datasets. There are a total of 15 hub genes that were selected from the PPI network constructed using the STRING application. Furthermore, these 15 genes were evaluated on PDA patients using TCGA and GTEx databases in (GEPIA). The online tool DAVID was used to analyse the functional annotation information for the DEGs. The functional pathway enrichment was performed on the GO and KEGG. The hub genes were mainly enriched for cell division, chromosome segregation, protein binding and microtubule binding. Further, the gene alteration study was performed using the cBioportal tool and screened out six hub genes (ASPM, CENPF, BIRC5, TTK, DLGAP5, and TOP2A) with a high alteration rate in PDA samples. Furthermore, Kaplan-Meier survival analysis was performed on the six hub genes and identified poor-survival outcomes that may be involved in tumorigenesis and PDA development. So, this study concludes that, these six hub genes may be potential prognostic biomarkers for PDA.

Investigation of the Shared Biomarkers in Heterotopic Ossification Between Ossification of the Ligamentum Flavum and Ankylosing Spondylitis.

Bioinformatics analysis of Gene Expression Omnibus (GEO). Ossification of the ligamentum flavum (OLF) and ankylosing spondylitis (AS) represent intricate conditions marked by the gradual progression of endochondral ossification. This investigation endeavors to unveil common biomarkers associated with heterotopic ossification and explore the potential molecular regulatory mechanisms. Microarray and RNA-sequencing datasets retrieved from the Gene Expression Omnibus (GEO) repository were harnessed to discern differentially expressed genes (DEGs) within the OLF and AS datasets. Subsequently, Weighted Gene Co-expression Network Analysis (WGCNA) was implemented to pinpoint co-expression modules linked to OLF and AS. Common genes were further subjected to an examination of functional pathway enrichment. Moreover, hub intersection genes were identified using the Least Absolute Shrinkage and Selection Operator (LASSO) regression, followed by an evaluation of diagnostic performance in external OLF and AS cohorts. Lastly, an analysis of immune cell infiltration was conducted to scrutinize the correlation of immune cell presence with shared biomarkers in OLF and AS. A total of 1353 and 91 Differentially Expressed Genes (DEGs) were identified in OLF and AS, respectively. Using the Weighted Gene Co-expression Network Analysis (WGCNA), 2 modules were found to be notably significant for OLF and AS. The integrative bioinformatic analysis revealed 3 hub genes (MAB21L2, MEGF10, ISLR) as shared risk biomarkers, with MAB21L2 being the central focus. Receiver Operating Characteristic (ROC) analysis exhibited a strong diagnostic potential for these hub genes. Gene Ontology (GO) analysis indicated their involvement in the positive regulation of myoblast proliferation. Notably, MAB21L2 was singled out as the optimal common biomarker for OLF and AS. Furthermore, an analysis of immune infiltration demonstrated a correlation between MAB21L2 expression and changes in immune cells. Activated CD8 T cells were identified as shared differential immune infiltrating cells significantly linked to MAB21L2 in both OLF and AS. This study represents the first instance of identifying MAB21L2 as a prospective diagnostic marker for patients contending with OLF associated with AS. The research results indicate that the ECM-receptor interaction and the cell-cell adhesion may play a role in both disease processes. This newfound knowledge not only enhances our understanding of the pathogenesis behind spinal ligament ossification but also uncovers potential targets for therapeutic interventions.

Systematic investigation of Radix Salviae for treating diabetic peripheral neuropathy disease based on network Pharmacology.

Diabetic peripheral neuropathy (DPN) is a debilitating complication of diabetes mellitus with limited available treatment options. Radix Salviae, a traditional Chinese herb, has shown promise in treating DPN, but its therapeutic mech-anisms have not been systematically investigated. Radix Salviae (Danshen in pinin), a traditional Chinese medicine (TCM), is widely used to treat DPN in China. However, the mechanism through which Radix Salviae treats DPN remains unclear. Therefore, we aimed to explore the mechanism of action of Radix Salviae against DPN using network pharmacology. The active ingredients and target genes of Radix Salviae were screened using the TCM pharmacology database and analysis platform. The genes associated with DPN were obtained from the Gene Cards and OMIM databases, a drug-com-position-target-disease network was constructed, and a protein-protein inter-action network was subsequently constructed to screen the main targets. Gene Ontology (GO) functional annotation and pathway enrichment analysis were performed via the Kyoto Encyclopedia of Genes and Genomes (KEGG) using Bioconductor. A total of 56 effective components, 108 targets and 4581 DPN-related target genes of Radix Salviae were screened. Intervention with Radix Salviae for DPN mainly involved 81 target genes. The top 30 major targets were selected for enrichment analysis of GO and KEGG pathways. These results suggested that Radix Salviae could treat DPN by regulating the AGE-RAGE signaling pathway and the PI3K-Akt signaling pathway. Therefore, Danshen may affect DPN by regulating inflammation and apoptosis.

Analysis and Validation of the Network Pharmacology of the Mechanism of Glycyrrhetinic Acid for the Treatment of High-Altitude Pulmonary Hypertension

Glycyrrhetinic acid (GA) is a pentacyclic triterpene component in Glycyrrhize glabra L, it has demonstrated an inhibitive effect on high-altitude pulmonary hypertension (HAPH), but the molecular action is still not known. We aimed to explore the mechanism of GA for the treatment of HAPH based on network pharmacology and molecular docking method. Cell experiment validation was also conducted. The targets for GA were screened using the Swiss Target Prediction and Batman databases. The HAPH-related targets were obtained using the GeneCards and OMIM databases. The common targets for diseases and drugs were obtained using a Venn diagram. The core targets were screened using the String database. Then, a component-target-disease diagram and protein–protein interaction (PPI) network mutual assistance diagram were developed using Cytoscape3.9.1 software. GO functional and KEGG pathway enrichment analyses were conducted using the Metascape database. Finally, molecular docking of the target and its corresponding active components were performed using Autodock software. A total of 68 common targets for glycyrrhetinic acid high-altitude pulmonary hypertension were screened out. The core targets include PTGS1, MMP1, SERPINA6, and nitric oxide synthase (NOS2), involving PPAR signal pathway, human cytomegalovirus infection, IL-17 signal pathway, proteoglycans in cancer, and other pathways. The molecular docking affinity was −8.4 kcal·mol−1 in average, indicating that GA has a good binding stability with key target proteins. In the PDGF-BB-induced PASMC proliferation model, PASMC proliferation and the p-p38, p38, p-ERK1/2, and ERK1/2 protein expression were inhibited. The pharmacological mechanism of GA for the treatment of HAPH was characterized by multi-target and multi pathway. GA may serve as a promising therapeutic candidate for HAPH but still needs further in vivo/in vitro experiment.

Mechanism of Si Ni San Combined with Astragalus in Treating Hepatic Fibrosis: A Network Pharmacology and Molecular Docking Study.

Si Ni San combined with Astragalus (SNSQ) has demonstrated significant efficacy in the treatment of hepatic fibrosis (HF), as confirmed by clinical practice. However, its pharmacological mechanism remains unclear. This study employs network pharmacology to identify key targets and proteins for molecular docking. Additionally, animal experiments were conducted to validate the network pharmacology results, providing further insights into the mechanism of SNSQ in treating HF. Effective compounds of SNSQ were screened from the Traditional Chinese Medicine Systems Pharmacology (TCMSP) and Encyclopedia of Traditional Chinese Medicine (ETCM) databases. Molecular formula structures of these effective compounds were obtained from the PubChem database. Partial target proteins with a probability greater than 0.6 were sourced from the SWISS database. Uniprot IDs corresponding to these target proteins were retrieved from the SUPERPRED database. The remaining target proteins of the compounds were obtained from the Uniprot database based on the Uniprot IDs. The drug target proteins were then summarized. Target points related to HF were selected from the GeneCards and OMIM databases. Common target points were identified in the Venn diagram and imported into Cytoscape 3.9.1 software to construct the "SNSQ-effective compound-target pathway-HF" network. AutoDock software was used for molecular docking of compounds and target proteins with high-degree values. The common target points underwent GO function enrichment and KEGG pathway enrichment analysis using the DAVID database. An HF rat model was established, and serum AST and ALT activities were measured. The Hyp assay kit was utilized to detect the Hyp content in liver tissue. To the transcription levels of pro-inflammatory factors (IL-1β, TNF-α, IL-6) and anti-inflammatory factors (IL-10, TGF-β1, IL-4) in rat serum and liver.IL-1β, TNF-α, IL-10, and TGF-β1 were chosen for validation through ELISA. Western blotting and qRT-PCR were used to assess the expression of related proteins, namely NFKB1, NF-κBp65, NF-κBp50, α-SMA, and Col-1 in liver tissue. qRT-PCR was also employed to study the expression of ECM synthesis and proliferation-related genes, including Cyclin D1, TIMP1, COL1A1 in HSC-T6 cells and rat liver tissue, as well as the inhibition of the ECM-related gene MMP13 in HSC-T6 cells and rat liver tissue. A total of 16 valid compounds were predicted, with kaempferol, sitosterol, and isorhamnetin exhibiting high-degree values. KEGG enrichment analysis revealed that the target genes of SNSQ were enriched in multiple pathological pathways, with the NF-Kappa B signaling pathway being predominant. Molecular docking simulations indicated strong affinities between SNSQ's primary components-kaempferol, sitosterol, isorhamnetin-and NFKB1. Experimental results demonstrated significant reductions in AST, ALT, and Hyp levels in the SNSQ group. Pro-inflammatory factors (IL-1β, TNF-ɑ) were markedly reduced, while anti-inflammatory factors (IL-10, TGF-β1) were substantially increased. The protein expression and transcription levels of α-SMA and Col-1 were significantly decreased, whereas those of NFKB1, NF-κBp65, and NF-κBp50 were notably elevated. mRNA expression levels of Cyclin D1, TIMP1, COL1A1 in HSC-T6 cells and rat liver tissue were significantly decreased, whereas MMP13 mRNA expression level was significantly increased. Treatment of HF with SNSQ involves multiple targets and pathways, with a close association with the overexpression of NFKB1 and activation of the NF-Kappa B signaling pathway. Its mechanism is closely linked to the activation of inflammatory responses, HSC activation, and proliferation.

Huangqin-Huanglian Decoction Protects Liver against Non-alcoholic Fatty Liver Disease in High Fat-diet Mice.

Traditional Chinese medicine (TCM) has the advantage of low toxicity of natural ingredients, multiple targets and effects, and low medication costs. It has unique advantages for metabolic and chronic diseases. Huangqin-Huanglian decoction (HQHLD) is composed of Scutellariae Radix, Coptidis Rhizoma, Rehmanniae Radix, and Gentianae Radix Et Rhozima; it has great potential for the treatment of NAFLD with the modern pharmacological research and TCM theory, but there is still a relative lack of research on the potential targets and pharmacological effects of HQHLD. In this work, we have used network pharmacology to predict the targets and signaling pathways of HQHLD, and validated NAFLD-related targets using the HFD model in order to explore more therapeutic drugs and methods for NAFLD. We collected the HQHLD ingredients and NAFLD targets through TCMSP, ETCM, DisGeNET, HGMD, MalaCards, OMIM, and TTD, built ingredients-target networks by Cytoscape, and screened key ingredients in HQHLD. DAVID and Metascape databases were used for GO functional enrichment analysis and KEGG pathway enrichment analysis, respectively. Molecular docking of the key ingredients and key targets was performed by AutoDock. We verified the effect of HQHLD on high-fat diet (HFD) mice by measuring the weight, liver weight index, and the level of TG, TC, LDL-C, and HDLC. HE staining and oil-red staining were performed to detect the damage and fat accumulation in the liver. The changes in INSR, PPAR-α, PPAR-γ, TNF-α, and caspase3 were experimented with WB. With the network pharmacology analysis, we found quercetin, baicalein, sitosterol, wogonin, oroxylin-A, glycyrrhizin, hydroberberine, berberine, sesamin, and carotene to be the main ingredients in HQHLD. According to KEGG pathway analysis, INSR, AKT, JNK1, PPAR-α, PPAR-γ, and the other 16 targets are the main targets of HQHLD in the treatment of NAFLD. We took HFD mice as the in vivo model of NAFLD. Our results showed that HQHLD could reduce liver weight, and TG and LDL-C levels, and increase HDL-C level in serum. By HE and oil red staining, we found that HQHLD could protect the morphology of hepatocytes and reduce fat in the liver. We also found HQHLD to protect the liver by increasing the expression of INSR and PPAR-α, and reducing the expression of PPAR-γ, TNF-α, and caspase3 in the liver. In conclusion, our study has firstly studied the main ingredients and key targets of HQHDL in treating NAFLD by network pharmacology analysis, and preliminarily confirmed that HQHLD could alleviate NAFLD in a multi-target way by lowering fatty acids, and decreasing insulin resistance, inflammation, and apoptosis in the liver.

IGF2BP3 participates in the pathogenesis of recurrent spontaneous abortion by regulating ferroptosis

The aberrant invasive capability of trophoblast cells is widely acknowledged as a primary mechanism underlying RSA. Recently, IGF2BP3 has been implicated in various cancers due to its influence on cellular invasion and migration. However, whether IGF2BP3 involve in the occurrence of RSA and the specific functions it assumes in the development of RSA remain elusive. In our study, we firstly collected villous tissues from RSA and those with normal pregnancies individuals to performed Protein sequencing and then detected the expression of IGF2BP3 through Western blot, qRT-PCR and immunohistochemistry. Secondly, we analyzed the single-cell data (GSE214607) to assess the expression of IGF2BP3 in invasive EVT trophoblasts. Thirdly, we utilized lentivirus technology to establish HTR-8/SVneo cell lines with stable IGF2BP3 knockdown and RNA-seq analysis was employed to investigate the GO functional pathway enrichment of IGF2BP3. Meanwhile, the effect of IGF2BP3 knockdown on trophoblast cells apoptosis, migration, and ferroptosis was evaluated through functional experiments. Additionally, LPS-induced abortion animal model was constructed to evaluate IGF2BP3 expression in placental tissues. A significant downregulation of IGF2BP3 was observed in the villous tissues of RSA patient, a finding corroborated by subsequent single cell sequencing results. Furthermore, it suggested that IGF2BP3 may be involved in the migration and apoptotic processes of trophoblast cells. Mechanistic research indicated that IGF2BP3 knockdown could compromise GPX4 mRNA stability, leading to the promotion of ferroptosis. Finally, our investigation observed the down-regulation of IGF2BP3 expression in placental villous tissues of an LPS-induced abortion animal model. Our findings revealed that IGF2BP3 was downregulated in the villous tissues of RSA patients. Mechanically, down-regulation of IGF2BP3 may induce RSA by promoting GPX4-mediated ferroptosis and inhibiting trophoblast invasion and migration. Our study may provide new targets and research directions for the pathogenesis of RSA.

Transcriptomic Analysis of Melatonin-Mediated Salt Stress Response in Germinating Alfalfa

Salt stress poses a significant threat to crop yields worldwide. Melatonin (MT), an endogenous hormone synthesized in plants, has emerged as a crucial player in plant responses to various abiotic stresses, including drought, salinity, heat, and cold. However, the precise molecular mechanisms underlying MT-mediated abiotic stress responses remain incompletely understood. To elucidate the key genes and pathways involved in MT-mediated alleviation of salt stress, we conducted physiological, biochemical, and transcriptomic analyses on alfalfa seedlings. Our results demonstrated that alfalfa seedlings treated with melatonin exhibited higher germination rates, longer bud lengths, and greater fresh weights compared to those subjected to salt stress alone. Furthermore, the levels of malondialdehyde (MDA) and superoxide anion (O2−) were reduced, while the activities and contents of antioxidant enzymes such as superoxide dismutase (SOD), peroxidase (POD), and glutathione (GSH) increased in response to melatonin treatment. Transcriptome analysis revealed 2181 differentially expressed genes (DEGs) in the salt-treated group, with 780 upregulated and 1401 downregulated genes. In contrast, the MT-treated group exhibited 4422 DEGs, including 1438 upregulated and 2984 downregulated genes. Functional annotation and pathway enrichment analysis indicated that DEGs were primarily involved in the biosynthesis of flavonoids, isoflavones, plant hormones, glutathione (GSH), soluble sugars, and other substances, as well as in ABC transporter and MAPK signaling pathways. Notably, the MT-treated group showed greater enrichment of DEGs in these pathways, suggesting that MT mitigates salt stress by modulating the expression of genes related to phytohormones and antioxidant capacity. Overall, our findings provide valuable insights into the molecular mechanisms underlying MT-mediated salt tolerance in alfalfa, with important implications for breeding salt-tolerant alfalfa and other crops.

P180 RNA sequencing of regulatory T cells sorted from Rheumatoid Arthritis peripheral blood

Abstract Background/Aims In the presence of an inflammatory environment, regulatory T cells (Tregs) have the ability to adopt an effector T helper(Th)-like phenotype, enabling them to produce pro-inflammatory cytokines, particularly IFN-γ (resembling Th1-like Tregs) and IL-17(resembling Th17-like Tregs). This adaptability of Tregs may be driven by the increased activation of effector pathways. In order to unveil distinct gene expression patterns while simultaneously characterizing Tregs based on their phenotype, we did RNA sequencing of sorted Tregs from the peripheral blood(PB) of individuals with rheumatoid arthritis(RA), osteoarthritis(OA), and healthy controls(HC). Methods Phenotyping of Tregs, Tregs expressing CXCR3(Th1-like Treg), CCR6(Th17-like Treg) from PB and synovial fluid(SF) of RA(n = 20 each), OA(n = 20 each) and HC (n = 20) was done using flow cytometry. Cell frequencies between the groups were compared using Student t test or one-way ANOVA (Tukey's post hoc Test) as appropriate; p &amp;lt; 0.05 was considered as statistically significant. RNA was extracted from the sorted PB Tregs of RA, OA and HC(n = 3 each) and was sequenced using Illumina Novaseq 6000 platform. Sample clustering was done using principal component analysis generated using DeSeq2. Genes were categorized as upregulated or downregulated based on their log2 fold change cut off of 2. Significantly upregulated and downregulated differentially expressed genes (DEGs) were shortlisted using a false discovery rate(FDR) corrected p-value cut-off of ≤ 0.05. Heatmaps and volcano plots were generated using R package. Functional pathway enrichment analysis was done using GO and KEGG software. Results Proportion of Tregs in the PB of RA was significantly lower than OA and HC. This difference was not observed in the SF. Both RA PB and SF exhibited higher percentages of Th1 and Th17-like Tregs when compared to OA and HC (Table 1A).DEG analysis revealed altered gene regulation in RA Tregs compared to Tregs from HC and OA. Compared to HC, genes related to TCR signaling, cytokine signaling, and chemokine signaling pathways were enriched in RA. In contrast, compared to OA, an upregulation of pathway genes such as NR4A1(PI3K-Akt signaling), and TNFRSF1B(TNF signaling) was noted in RA. Furthermore, genes associated with Th1 differentiation and activation (IFNGR1, IL12RB1, STAT4) and genes associated with Th17 profile(RORC, STAT3, IL6R, IL17B) were upregulated in RA relative to HC and but not in OA compared to to HC.(Table 1B). Additionally, genes involved in glucose metabolism, purine metabolism, oxidative phosphorylation, and pyrimidine metabolism exhibited distinct patterns of regulation in RA. Conclusion Our Treg phenotyping data shows that there is higher conversion of Tregs into Th1 and Th17 like Tregs in RA and transcriptome data indicates cytokine(IFNγ, IL-6, IL-12 &amp; IL-17) mediated transcriptional changes likely drove the conversion of Tregs towards Th1 and Th17 phenotype. Our data also suggests a possible metabolic re-programming of Tregs towards an inflammatory phenotype affecting its stability. Disclosure V. Kommoju: None. C. mariaselvam: None. S. Vembar: None. S. Bulusu: None. C. Kavadichanda: None. M. Thabah: None. V. Negi: None.

Study on the regulation of trophoblast activity by abnormally expressed hsa_circ_0024838/miR-543/HIF1A in patients with gestational diabetes mellitus

IntroductionThis study aimed to screen circRNAs involved in gestational diabetes mellitus (GDM)-related macrosomia. One differentially expressed circRNA (DEC), hsa_circ_0024838, was further tested for its potential role and mechanism in trophoblasts. MethodsDECs in GDM were selected through GSE182737 and GSE194119. The targets were predicted for DECs and microRNAs (miRNAs), to complete the construction of the circRNA-miRNA-gene network. Functional annotation and related biological pathway enrichment analysis were performed on the target genes of miRNAs in the network. Subsequently, the expression levels of hsa_circ_0024838, miR-543, and HIF1A mRNA were identified by real-time quantitative real-time PCR (RT-qPCR) in GDM patients. Trophoblast activity was assessed via CCK-8 assay, apoptosis assay, and Matrigel invasion assay. Finally, interactions between miR-543 and either hsa_circ_0024838 or HIF1A were confirmed using dual-luciferase reporter assays. ResultsA GDM-related circRNA-miRNA-genes interaction network was constructed, consisting of 35 circRNAs, 46 miRNAs, and 122 target genes. Functional enrichment revealed that the enriched pathways were involved in GDM. Hsa_circ_0024838 and HIF1A mRNA expression levels were upregulated in GDM, while miR-543 expression levels were downregulated. A significant positive correlation between hsa_circ_0024838 and newborn weight was observed. Both hsa_circ_0024838 and HIF1A possessed binding sites for miR-543. Overexpressing hsa_circ_0024838 in high-glucose (HG)-cultured trophoblasts can partially reverse HG-induced reduction in trophoblast cell proliferation/migration and increase apoptosis. But this reversal can be negated by co-transfection with miR-543 mimics. The effects of miR-543 can be counteracted by HIF1A. DiscussionHsa_circ_0024838 can regulate the expression of HIF1A by interacting with miR-543. This regulates the HIF1A signaling pathway and enhance vitality in trophoblast cells.

Longitudinal profiling of human androgenotes through single-cell analysis unveils paternal gene expression dynamics in early embryo development.

How do transcriptomics vary in haploid human androgenote embryos at single cell level in the first four cell cycles of embryo development? Gene expression peaks at the fourth cell cycle, however some androcytes exhibit unique transcriptional behaviors. The developmental potential of an embryo is determined by the competence of the oocyte and the sperm. However, studies of the contribution of the paternal genome using pure haploid androgenotes are very scarce. This study was performed analyzing the single-cell transcriptomic sequencing of 38 androcytes obtained from 10 androgenote bioconstructs previously produced in vitro (de Castro et al., 2023). These results were analyzed through different bioinformatics software such as g: Profiler, GSEA, Cytoscape, and Reactome. Single cell sequencing was used to obtain the transcriptomic profiles of the different androcytes. The results obtained were compared between the different cycles studied using the DESeq2 program and functional enrichment pathways using g: Profiler, Cytoscape, and Reactome. A wave of paternally driven transcriptomic activation was found during the third-cell cycle, with 1128 upregulated and 225 downregulated genes and the fourth-cell cycle, with 1373 upregulated and 286 downregulated genes, compared to first-cell cycle androcytes. Differentially expressed routes related to cell differentiation, DNA-binding transcription, RNA biosynthesis and RNA polymerase II transcription regulatory complex, and cell death were found in the third and fourth with respect to the first-cell cycle. Conversely, in the fourth cell cycle, 153 downregulated and 332 upregulated genes were found compared with third cell cycle, associated with differentially expressed processes related to E-box binding and zinc finger protein 652 (ZNF652) transcription factor. Further, significant overexpression of LEUTX, PRAMEF1, DUXA, RFPL4A, TRIM43, and ZNF675 found in androgenotes, compared to biparental embryos, highlights the paternal contributions to zygote genome activation. All raw sequencing data are available through the Gene Expression Omnibus (GEO) under accessions number: GSE216501. Extrapolation of biological events from uniparental constructs to biparental embryos should be done with caution. Maternal and paternal genomes do not act independently of each other in a natural condition. The absence of one genome may affect gene transcription of the other. In this sense, the haploid condition of the bioconstructs could mask the transcriptomic patterns of the single cells. The results obtained demonstrated the level of involvement of the human paternal haploid genome in the early stages of embryo development as well as its evolution at the transcriptomic level, laying the groundwork for the use of these bioconstructs as reliable models to dispel doubts about the genetic role played by the paternal genome in the early cycles of embryo development. This study was funded by Instituto de Salud Carlos III (ISCIII) through the project 'PI22/00924', co-funded by European Regional Development Fund (ERDF); 'A way to make Europe'. F.D. was supported by the Spanish Ministry of Economy and Competitiveness through the Miguel Servet program (CPII018/00002). M.J.E. was supported by Instituto de Salud Carlos III (PI19/00577 [M.J.E.]) and FI20/00086. P.dC. was supported by a predoctoral grant for training in research into health (PFIS PI19/00577) from the Instituto de Salud Carlos III. All authors declare having no conflict of interest with regard to this trial.

Intestinal cell diversity and treatment responses in a parasitic nematode at single cell resolution

BackgroundParasitic nematodes, significant pathogens for humans, animals, and plants, depend on diverse organ systems for intra-host survival. Understanding the cellular diversity and molecular variations underlying these functions holds promise for developing novel therapeutics, with specific emphasis on the neuromuscular system’s functional diversity. The nematode intestine, crucial for anthelmintic therapies, exhibits diverse cellular phenotypes, and unraveling this diversity at the single-cell level is essential for advancing knowledge in anthelmintic research across various organ systems.ResultsHere, using novel single-cell transcriptomics datasets, we delineate cellular diversity within the intestine of adult female Ascaris suum, a parasitic nematode species that infects animals and people. Gene transcripts expressed in individual nuclei of untreated intestinal cells resolved three phenotypic clusters, while lower stringency resolved additional subclusters and more potential diversity. Clusters 1 and 3 phenotypes displayed variable congruence with scRNA phenotypes of C. elegans intestinal cells, whereas the A. suum cluster 2 phenotype was markedly unique. Distinct functional pathway enrichment characterized each A. suum intestinal cell cluster. Cluster 2 was distinctly enriched for Clade III-associated genes, suggesting it evolved within clade III nematodes. Clusters also demonstrated differential transcriptional responsiveness to nematode intestinal toxic treatments, with Cluster 2 displaying the least responses to short-term intra-pseudocoelomic nematode intestinal toxin treatments.ConclusionsThis investigation presents advances in knowledge related to biological differences among major cell populations of adult A. suum intestinal cells. For the first time, diverse nematode intestinal cell populations were characterized, and associated biological markers of these cells were identified to support tracking of constituent cells under experimental conditions. These advances will promote better understanding of this and other parasitic nematodes of global importance, and will help to guide future anthelmintic treatments.

The m6A methyltransferase METTL5 promotes neutrophil extracellular trap network release to regulate hepatocellular carcinoma progression.

Hepatocellular carcinoma (HCC) is one of the most common malignant tumors worldwide, it has a poor prognosis due to its highly invasive and metastatic nature. Consequently, identifying effective prognostic markers and potential therapeutic targets has been extensively investigated. METTL5, an 18S rRNA methyltransferase, is abnormally high in HCC. But its biological function and prognostic significance in HCC remain largely unelucidated. This study aimed to investigate the role of METTL5 in HCC progression, and elucidate its possible molecular mechanisms in HCC via transcriptome sequencing, providing new insights for identifying new HCC prognostic markers and therapeutic targets. The METTL5 expression in HCC and paracancerous tissues was analyzed using HCC immunohistochemical microarrays and bioinformatic retrieval methods to correlate METTL5 with clinicopathological features and survival prognosis. We constructed a METTL5 knockdown hepatocellular carcinoma cell line model and an animal model to determine the effect of METTL5 on hepatocellular carcinoma progression. Subsequently, RNA sequencing was performed to analyze the molecular mechanism of METTL5 in HCC based on the sequencing results, and relevant experiments were performed to verify it. We found that METTL5 expression was elevated in hepatocellular carcinoma tissues and correlated with poor patient prognosis, and in the analysis of clinicopathological features showed a correlation with TNM staging. In hepatocellular carcinoma cell lines with knockdown of METTL5, the malignant biological behavior was significantly reduced both invitro and invivo. Based on the sequencing results as well as the results of GO functional enrichment analysis and KEGG pathway enrichment analysis, we found that METTL5 could promote the generation and release of neutrophil extracellular capture network (NETs) and might further accelerate the progression of HCC. The m6A methyltransferase METTL5 is overexpressed in hepatocellular carcinoma (HCC) and correlates with poor prognosis. METTL5 accelerates malignant progression of HCC by promoting generation and release of the neutrophil extracellular traps (NETs) network, providing new insights for clinical biomarkers and immunotherapeutic targets in HCC prognosis.

Waterpipe smoke condensate induces epithelial-mesenchymal transformation and promotes metastasis of oral cancer by FOXD1 expression

Background/purposeSmoking is a major contributor to global oral cancer cases, necessitating urgent intervention. FOXD1, involved in developmental processes and various cancers, shows promise as a prognostic marker in oral squamous cell carcinoma (OSCC). This study investigates the impact of waterpipe smoke condensate (WPSC) on OSCC, focusing on FOXD1 role in inducing epithelial-mesenchymal transition (EMT) and metastasis. MethodsThe study involved using OSCC cells treated with WPSC to evaluate their proliferation, colony formation, gene expression, and protein levels. The researchers also explored the clinical relevance of their findings using online databases to analyze FOXD1 expression in cancer tissues and its correlation with clinicopathological features and patient survival. Additionally, in silico tools were employed for functional analysis, pathway enrichment, and network exploration. ResultsThe study found that WPSC increased the expression of FOXD1 in OSCC cells, which led to increased cell growth. The study also showed that FOXD1 plays a critical role in the EMT process induced by WPSC, as evidenced by changes in the expression of EMT-related genes and proteins. Clinical analysis revealed that FOXD1 was significantly associated with more aggressive tumor features and poorer prognosis in cancer patients. ConclusionThe study highlights FOXD1 as a key player in OSCC pathogenesis and a potential prognostic marker and therapeutic target, particularly when influenced by WPSC exposure. Further research is needed to explore FOXD1 molecular mechanisms and clinical implications to enhance OSCC treatment strategies.