PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “transcriptomic profiling”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16Linked to original sources

Integrative multi-omics analyses suggest a candidate microbial metabolite-associated host gene network in ulcerative colitis.

Ulcerative colitis (UC) is associated with gut microbial dysbiosis, but the host molecular alterations potentially linked to microbially derived metabolites remain incompletely understood. We integrated Mendelian randomization (MR), microbial metabolite annotation, computational target prediction, colonic transcriptomics, network analysis, and machine learning. MiBioGen microbiome GWAS data were used as exposures and FinnGen Release 12 ULCERENTER as the outcome. Metabolites linked to MR-prioritized taxa were retrieved from GutMGene, and human targets were predicted using SwissTargetPrediction and SEA. UC-related genes were defined by integrating differential expression analysis and WGCNA and then intersected with predicted metabolite targets. MR prioritized one family and eight genera showing nominal genetically supported associations with UC, but none remained significant after Benjamini-Hochberg FDR correction. Three prioritized genera were linked to 15 microbe-metabolite records, corresponding to 13 unique metabolites; nine were retained for target prediction, yielding 277 unique predicted human targets. Transcriptomic analysis identified 1,530 DEGs and a 312-gene MEgrey60 module, with 273 overlapping genes, producing 1,569 unique UC-related genes. Their intersection with the 277 predicted targets yielded 47 candidate genes. Enrichment analyses highlighted mainly metabolic and lipid-related processes. Random Forest showed the highest mean AUC across the two independent external benchmarking cohorts, and SHAP prioritized EPHX1, HSD17B2, IGFBP5, and MMP10. IBDome analysis showed inflammation-associated expression differences in these genes. This study provides a genomics-informed, hypothesis-generating framework that prioritizes candidate microbe-metabolite-host relationships in UC for future experimental validation.

Humans↗

De novo clustering of large long-read transcriptome datasets with isONclust3.

MOTIVATION: Long-read sequencing techniques can sequence transcripts from end to end, greatly improving our ability to study the transcription process. Although there are several well-established tools for long-read transcriptome analysis, most are reference-based. This limits the analysis of organisms without high-quality reference genomes and samples or genes with high variability (e.g. cancer samples or some gene families). In such settings, analysis using a reference-free method is favorable. The computational problem of clustering long reads by region of common origin is well-established for reference-free transcriptome analysis pipelines. Such clustering enables large datasets to be split roughly by gene family and, therefore, an independent analysis of each cluster. There exist tools for this. However, none of those tools can efficiently process the large amount of reads that are now generated by long-read sequencing technologies. RESULTS: We present isONclust3, an improved algorithm over isONclust and isONclust2, to cluster massive long-read transcriptome datasets into gene families. Like isONclust, isONclust3 represents each cluster with a set of minimizers. However, unlike other approaches, isONclust3 dynamically updates the cluster representation during clustering by adding high-confidence minimizers from new reads assigned to the cluster and employs an iterative cluster-merging step. We show that isONclust3 yields results with higher or comparable quality to state-of-the-art algorithms but is 10-100 times faster on large datasets. Also, using a 256 Gb computing node, isONclust3 was the only tool that could cluster 37 million PacBio reads, which is a typical throughput of the recent PacBio Revio sequencing machine. AVAILABILITY AND IMPLEMENTATION: https://github.com/aljpetri/isONclust3.

Algorithms↗

Spatial transcriptomics-aided localization for single-cell transcriptomics with STALocator.

Single-cell RNA-sequencing (scRNA-seq) techniques can measure gene expression at single-cell resolution but lack spatial information. Spatial transcriptomics (ST) techniques simultaneously provide gene expression data and spatial information. However, the data quality of the spatial resolution or gene coverage is still much lower than the quality of the single-cell transcriptomics data. To this end, we develop a ST-Aided Locator for single-cell transcriptomics (STALocator) to localize single cells to corresponding ST data. Applications on simulated data showed that STALocator performed better than other localization methods. When applied to the human brain and squamous cell carcinoma data, STALocator could robustly reconstruct the relative spatial organization of critical cell populations. Moreover, STALocator could enhance gene expression patterns for Slide-seqV2 data and predict genome-wide gene expression data for fluorescence in situ hybridization (FISH) and Xenium data, leading to the identification of more spatially variable genes and more biologically relevant Gene Ontology (GO) terms compared with the raw data. A record of this paper's transparent peer review process is included in the supplemental information.

Single-Cell Analysis↗

Genome-wide DNA methylation and transcriptome sequencing analyses of lens tissue in an age-related mouse cataract model.

DNA methylation is known to be associated with cataracts. In this study, we used a mouse model and performed DNA methylation and transcriptome sequencing analyses to find epigenetic indicators for age-related cataracts (ARC). Anterior lens capsule membrane tissues from young and aged mice were analyzed by MethylRAD-seq to detect the genome-wide methylation of extracted DNA. The young and aged mice had 76,524 and 15,608 differentially methylated CCGG and CCWGG sites, respectively. The Pearson correlation analysis detected 109 and 33 differentially expressed genes (DEGs) with negative methylation at CCGG and CCWGG sites, respectively, in their promoter regions. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses showed that DEGs with abnormal methylation at CCGG sites were primarily associated with protein kinase C signaling (Akap12, Capzb), protein threonine kinase activity (Dmpk, Mapkapk3), and calcium signaling pathway (Slc25a4, Cacna1f), whereas DEGs with abnormal methylation at CCWGG sites were associated with ribosomal protein S6 kinase activity (Rps6ka3). These genes were validated by pyrosequencing methylation analysis. The results showed that the ARC group (aged mice) had lower Dmpk and Slc25a4 methylation levels and a higher Rps6ka3 methylation than the control group (young mice), which is consistent with the results of the joint analysis of differentially methylated and differentially expressed genes. In conclusion, we confirmed the genome-wide DNA methylation pattern and gene expression profile of ARC based on the mouse cataract model with aged mice. The identified methylation molecular markers have great potential for application in the future diagnosis and treatment of ARC.

Animals↗

An anti-androgen resistance-related gene signature acts as a prognostic marker and increases enzalutamide efficacy via PLK1 inhibition in prostate cancer.

BACKGROUND: Anti-androgen resistance remains a major clinical challenge in the treatment of prostate cancer (PCa), leading to disease progression and treatment failure. Despite extensive research on resistance mechanisms, a reliable prognostic model for predicting patient outcomes and guiding therapeutic strategies is still lacking. This study aimed to develop a novel gene signature related to anti-androgen resistance and evaluate its prognostic and therapeutic implications. METHODS: Anti-androgen resistance-related differentially expressed genes (ARRDEGs) were identified through transcriptomic analysis of enzalutamide- and dual enzalutamide abiraterone-resistant PCa cell lines from the GEO database. Functional enrichment analysis was performed to determine the biological roles of these genes. A prognostic gene signature was developed using univariate Cox regression, LASSO, and multivariate Cox regression models. The model was validated in independent PCa cohorts from The Cancer Genome Atlas (TCGA). Additionally, we assessed the correlation between the signature, immune infiltration, immune checkpoint expression, and drug sensitivity. The efficacy of PLK1 inhibition combined with enzalutamide was further explored using in vitro and in vivo experiments. RESULTS: We identified 304 ARRDEGs, from which three key genes (LMNB1, SSPO, and PLK1) were selected to construct a prognostic signature. This gene signature effectively stratified PCa patients into high- and low-risk groups, with the high-risk group exhibiting shorter recurrence-free survival and distinct immune characteristics. High-risk patients demonstrated elevated immune checkpoint expression (B7H3, CTLA-4, B7-1, and TIGIT), increased M2 macrophage infiltration, and enhanced sensitivity to chemotherapy and targeted therapy. Mechanistically, PLK1 inhibition potentiated the antitumor effect of enzalutamide by downregulating SLC7A11 and inducing ferroptosis, providing a potential therapeutic strategy to overcome anti-androgen resistance. CONCLUSION: We established a novel ARRDEGs-based prognostic signature that predicts PCa progression and response to chemotherapy and targeted therapy. The integration of this signature with immune profiling and drug sensitivity analysis provides a valuable tool for precision oncology in PCa. Our findings highlight the potential of PLK1 inhibition as a therapeutic strategy to enhance enzalutamide efficacy and overcome resistance.

Humans↗

Liver transcriptome analysis revealed multiple immune processes and lipid metabolism pathways involved in the defense response of the turbot (Scophthalmus maximus) against Aeromonas salmonicida.

Aeromonas salmonicida is a significant pathogen causing notable economic losses in Scophthalmus maximus aquaculture. This study utilized Illumina sequencing technology to examine the transcriptional response characteristics of S. maximus liver at 24 h following A. salmonicida infection. A total of 2363 differentially expressed genes (DEGs) were identified when compared to the negative control group. The immunity-related Toll-like receptor signaling pathway, NOD-like receptor signaling pathway, as well as metabolism-related PPAR signaling pathway and insulin signaling pathway, were notably enriched. Significant differences exist in the expression of key genes within the PPAR pathway, particularly cd36, acsl4a, pparαa, and plin2, all of which mediate the interaction between lipid metabolism and the immune response. These results offer valuable insights into the immunometabolic regulatory mechanism of S. maximus response to A. salmonicida infection.

Animals↗

Histology-Based Virtual RNA Inference Identifies Pathways Associated With Metastasis Risk in Colorectal Cancer.

Colorectal cancer (CRC) remains a major health concern, with >150,000 new diagnoses and >50,000 deaths annually in the United States, underscoring an urgent need for improved screening, prognostication, disease management, and therapeutic approaches. The tumor microenvironment (TME)-comprising cancerous and immune cells interacting within the tumor's spatial architecture-plays a critical role in disease progression and treatment outcomes, reinforcing its importance as a prognostic marker for metastasis and recurrence risk. However, traditional methods for TME characterization, such as bulk transcriptomics and multiplex protein assays, lack sufficient spatial resolution. Although spatial transcriptomics (ST) allows for the high-resolution mapping of whole transcriptomes at near-cellular resolution, current ST technologies (eg, Visium and Xenium) are limited by high costs, low throughput, and issues with reproducibility, preventing their widespread application in large-scale molecular epidemiology studies. In this study, we refined and implemented virtual RNA inference (VRI) to derive ST-level molecular information directly from hematoxylin and eosin (H&E)-stained tissue images. Our VRI models were trained on the largest matched CRC ST data set to date, comprising 45 patients and >300,000 Visium spots from primary tumors. Using state-of-the-art deep learning models (UNI, ResNet-50, Vision Transformer, and Vision Mamba), we achieved a median Spearman's correlation coefficient of 0.546 between predicted and measured spot-level expression. As validation, VRI-derived gene signatures linked to specific tissue regions (tumor, interface, submucosa, stroma, serosa, muscularis, and inflammation) showed strong concordance with signatures generated via direct ST, and VRI performed accurately in estimating cell-type proportions spatially from H&E slides. In an expanded CRC cohort controlling for tumor invasiveness and clinical factors, we further identified VRI-derived gene signatures significantly associated with key prognostic outcomes, including metastasis status. Although certain tumor-related pathways are not fully captured by histology alone, our findings highlight the ability of VRI to infer a wide range of "histology-associated" biological pathways at near-cellular resolution without requiring ST profiling. Future efforts will extend this framework to expand TME phenotyping from standard H&E tissue images, with the potential to accelerate translational CRC research at scale.

Humans↗

Higher Expression of HPV16 Derived E7_LI Transcript Observed in Men With HIV and Recurrent Anal Cancer.

Squamous cell carcinoma of the anus (SCCA) or anal cancer (AC) is an understudied cancer with a high occurrence rate in people with HIV (PWH), especially men having sex with men (MSM). Furthermore, AC recurs in approximately one-fourth of patients who undergo standard care with chemoradiation therapy (CRT). Using bulk RNA sequencing data of AC obtained from 12 patients with non-recurrent (NR, N&#x2009;=&#x2009;9) or recurrent (R, N&#x2009;=&#x2009;3) cancer, we previously showed upregulated expression of key immune genes in the NR compared to the R group. Although the main causative agent of AC is high-risk human papillomavirus (HPV), association of host and viral RNA transcript expression contributing to AC recurrence has not been extensively studied. The objective of the current study was to determine whether enrichment of specific HPV genotypes and/or HPV gene expression patterns differentiate the two groups and if any specific viral (HPV) and host (human) immune mediators correlate with each other. Using bulk RNA sequencing data and VIRTUS 2, we detected viral RNA reads mapping to seven high-risk and six low-risk HPV types, of which the high-risk HPV16 observed in 83% (10/12) AC tumors (7/9 NR and 3/3 R). Rate of all HPV genomes trended toward a decrease in NR AC isolates and correlation between HPV types was more commonly observed in low-risk ones. Analysis of HPV 16 gene expression profile showed a significantly lower positivity rate for a polycistronic transcript encoding for E7^L1 in the NR group (1/9, NR vs. 3/3, R, p&#x2009;<&#x2009;0.05). An unbiased correlation analysis of HPV-human transcript expression showed a direct correlation between HPV transcripts and human genes involved in cell growth. The data also identified human transcripts showing an inverse correlation with HPV gene expression. These included genes involved in negative regulation of growth, proliferation, and immune response. Taken together, these data indicate that concurrent analyses of viral and host factors in the same tumor can identify potential new therapeutic targets to ameliorate cancer recurrence post-treatment.

Humans↗

iSoMAs: Finding isoform expression and somatic mutation associations in human cancers.

Aberrant alternative splicing, prevalent in cancer, impacts various cancer hallmarks involving proliferation, angiogenesis, and invasion. Splicing disruption often results from somatic point mutations rewiring functional pathways to support cancer cell survival. We introduce iSoMAs (iSoform expression and somatic Mutation Association), an efficient computational pipeline leveraging principal component analysis technique, to explore how somatic mutations influence transcriptome-wide gene expression at the isoform level. Applying iSoMAs to 33 cancer types comprising 9,738 tumor samples in The Cancer Genome Atlas, we identified 908 somatically mutated genes significantly associated with altered isoform expression across three or more cancer types. Mutations linked to differential isoform expression occurred through both cis- and trans-acting mechanisms, involving well-known oncogenes/suppressor genes, RNA binding protein and splicing factor genes. With wet-lab experiments, we verified direct association between TP53 mutations and differential isoform expression in cell cycle genes. Additional iSoMAs genes have been validated in the literature with independent cohorts and/or methods. Despite the complexity of cancer, iSoMAs attains computational efficiency via dimension reduction strategy and reveals critical associations between regulatory factors and transcriptional landscapes.

Humans↗

PoweREST: Statistical power estimation for spatial transcriptomics experiments to detect differentially expressed genes between two conditions.

Recent advancements in spatial transcriptomics (ST) have significantly enhanced biological research in various domains. However, the high cost for current ST data generation techniques restricts the large-scale application of ST. Consequently, maximization of the use of available resources to achieve robust statistical power for ST data is a pressing need. One fundamental question in ST analysis is detection of differentially expressed genes (DEGs) under different conditions using ST data. Such DEG analyses are performed frequently, but their power calculations are rarely discussed in the literature. To address this gap, we developed PoweREST, a power estimation tool designed to support the power calculation for DEG detection with 10X Genomics Visium data. PoweREST enables power estimation both before any ST experiments and after preliminary data are collected, making it suitable for a wide variety of power analyses in ST studies. We also provide a user-friendly, program-free web application that allows users to interactively calculate and visualize study power along with relevant parameters.

Gene Expression Profiling↗

Partners in root nodule symbiosis respond uniquely to heavy metal stresses in a host genotype-dependent manner.

The mutualistic symbiosis between legume roots and soil rhizobia culminates in the formation of root nodules, where nitrogen is fixed. Root nodule symbiosis is inhibited by heavy metal stress. In this study, we investigated the relative responses of the symbiotic partners to a non-essential heavy metal cadmium (Cd) and an essential heavy metal zinc (Zn) stress and identified patterns in gene expression. We performed dual transcriptomics in nodules, using the Medicago truncatula-Sinorhizobium meliloti symbiotic system. Phenotypes were measured in the wild-type Medicago truncatula and a mutant in an ABC transporter gene (Mtabcg36), which showed compromised nodule formation in control conditions and further after heavy metal treatment. We observed that the rhizobia were particularly sensitive to Zn in mutant nodules. The greatest degree of differential gene expression in the host plant were observed under Cd and Zn treatments in wild-type nodules. Most Cd-regulated host genes were also differentially regulated by Zn, revealing little discernment between an essential and a non-essential ion under increased exposure. Furthermore, the host response to both the stresses affected auxin and iron homeostasis genes in a host genotype-dependent manner. Our results suggested impaired cadmium export from the mutant nodules. These results have potential implications in agricultural management systems and bioremediation strategies.

Symbiosis↗

Differential expression of neuronal function genes follows a tissue-specific temporal dynamic during Deformed Wing Virus infection in honey bees.

Deformed Wing Virus type A (DWV-A) is one of the primary threats to honeybees (Apis mellifera), significantly impacting their nervous system physiology and behavior. While its neurotropic nature is well-recognized, the temporal dynamics of the neuronal transcriptomic response following oral infection, the natural transmission route, remains poorly understood. In this study, we analyzed gene expression in the heads of worker bees orally inoculated with DWV-A over a 16-day time course (1, 4, 7, 10, 13, and 16 days post-inoculation). RNA-seq analysis at day 10 identified 147 differentially expressed genes associated with different biological processes that are critical to the organism, including cellular metabolism and neuronal activity. RT-qPCR validation revealed a persistent downregulation of key genes related to glutamatergic system (eaat-2, neto, and kainate) and sensory perception-related genes in the antennae. Notably, the simultaneous co-expression of nurse-associated and forager-associated marker genes suggests that DWV-A infection induces an asynchrony in behavioral maturation. Our findings demonstrate that DWV-A disrupts neuronal homeostasis and peripheral sensory perception in a tissue-specific and time-dependent manner, providing a molecular framework to understand the behavioral impairment and the loss of coordination at the colony level.

Animals↗

Multi-omics unveils seasonal remodeling and metabolic crosstalk between testis and abdominal fat body in a non-amplexus stream frog Nanorana taihangnica (Anura: Dicroglossidae).

BACKGROUND: Energy allocation between reproduction and survival represents a fundamental life-history challenge for animals in seasonal environments. Using integrated transcriptomics and metabolomics, we investigated Nanorana taihangnica (Anura: Dicroglossidae), a non-amplexus stream frog endemic to China, to elucidate the seasonal morphological and molecular coordination between the testis and abdominal fat body. RESULTS: Morphological analysis showed that fat body adipocyte cross-sectional area minimized at the end of the breeding season but rapidly recovered thereafter, while testicular volume continued declining post-breeding and only recovered during the non-breeding period. During breeding season, multi-omics analyses revealed that the fat body enhanced fatty acid oxidation, upregulated histidine-carnosine metabolism, activated NAD+ metabolism and FOXO3-mediated antioxidative responses to mitigate metabolic stress, and regulated adipocyte survival and apoptosis via sphingolipid signaling. Seasonal testicular development was centrally regulated by the mTOR signaling pathway, whose activity integrated autophagy levels, NAD+ availability, and aspartate metabolism to coordinate spermatogonial proliferation and spermatogenesis. CONCLUSIONS: This study demonstrates that N. taihangnica optimizes seasonal energy storage, allocation, and reproductive investment through molecular and metabolic crosstalk between the fat body and testis, providing empirical insights into the physiological integration of life-history strategies in animals inhabiting fluctuating environments.

Animals↗

Nonadditive gene expression and reduced homoeolog expression bias in an intraspecific hexaploid wheat hybrid.

BACKGROUND: Intraspecific hybridization in allopolyploid plants can generate additive and nonadditive changes in gene expression through interactions between divergent parental genomes. However, how it simultaneously affects gene expression and the relative expression of homoeologs in higher-order polyploids is less well understood. To study this, we sequenced seedling leaf transcriptomes and profiled gene body methylation in two hexaploid wheat (Triticum aestivum L.) cultivars and their F&#x2081; hybrids. RESULTS: Although only 4.3% of genes differed in expression between the parents, 22.3% deviated from mid-parent expression in the hybrids, with many showing transgressive expression. 32.1% of triads contained at least one homoeolog that deviated from mid-parent expression, and all three homoeologs deviated in 11% of triads, substantially more than expected by chance. Triads in which all three homoeologs were overexpressed also showed reduced differences in expression among homoeologs. Greater parental divergence in relative homoeolog expression was associated with nonadditive expression. Genes lacking gene body methylation were also more likely to show dominant or transgressive expression, whereas gene body methylation was associated with more balanced homoeolog expression and additive or conserved expression. CONCLUSIONS: Intraspecific hybridization in hexaploid wheat, even without a change in ploidy, was associated with widespread nonadditive gene expression and altered relative homoeolog expression within triads. These responses were associated with parental differences in homoeolog expression and the absence of gene body methylation. Although our findings are limited to seedling leaves from a single intraspecific cross, they provide a basis for testing the generality of these patterns across tissues, developmental stages, and genetic backgrounds.

Triticum↗

Comparative transcriptomic analysis of the gills and hepatopancreas of freshwater-cultured Litopenaeus vannamei under chronic nitrite stress.

To investigate the differences in molecular responses between the gills and hepatopancreas of freshwater-cultured Litopenaeus vannamei under chronic nitrite stress, a 30-day chronic stress experiment was conducted with a control group and a stress group. Transcriptomic analysis of the gills and hepatopancreas was performed using Illumina sequencing; differentially expressed genes (DEGs) were identified, and GO, KEGG, GSEA, PPI, and RT-qPCR validation were carried out. The results showed that 196 DEGs (161 up-regulated and 35 down-regulated) were identified in the gills, and 287 DEGs (199 up-regulated and 88 down-regulated) in the hepatopancreas, with only 18 DEGs shared between the two tissues. DEGs in the gills were enriched in oxidoreductase activity, glycerophospholipid metabolism, and tyrosine metabolism; DEGs in the hepatopancreas were enriched in lipid transporter activity, phagosome, ECM-receptor interaction, and riboflavin metabolism. GSEA revealed significant suppression of the mTOR pathway in the gills and the Polycomb complex pathway in the hepatopancreas. PPI network analysis identified hub genes P5CS and eEF2 in the gills, and PER, TUBB1, SHMT, and TUBB4B in the hepatopancreas. RT-qPCR validation was consistent with the RNA-seq results (R2&#xa0;=&#xa0;0.764). This study indicates that, under chronic nitrite stress, the gill response is centered on redox regulation and inhibition of growth metabolism, whereas the hepatopancreas response primarily involves lipid transport, cytoskeletal remodeling, and phagosome activation. The two tissues synergistically adapt through fundamental biosynthetic and motor protein pathways. This research provides molecular evidence for deciphering the nitrite tolerance mechanisms in freshwater-cultured shrimp.

Animals↗

Identification of Critical Genes for Recurrent Aphthous Ulcer by Transcriptome Data Analysis and Mendelian Randomization.

PURPOSE: Recurrent aphthous ulcer (RAU) is a common oral mucosal disorder with a poorly understood etiology, significantly affecting patients' quality of life. This study aims to investigate critical genes linked to RAU and explore their biological mechanisms using transcriptomic data and Mendelian randomization (MR) analysis. MATERIALS AND METHODS: RAU-related gene expression data from the GEO database (GSE37265) were analyzed to identify differentially expressed genes (DEGs). A two-sample MR approach was used to assess the causal impact of expression quantitative trait loci (eQTL) on RAU. Critical genes were identified by intersecting DEGs with significant MR findings. GO and KEGG pathway enrichment analyses were performed, along with GSEA and immune cell infiltration analysis, to investigate the functions and mechanisms of these genes in RAU. RESULTS: A total of 184 differentially expressed genes (DEGs) were identified, while 339 RAU-associated genes were screened through MR analysis. Cross-validation further identified 7 critical genes. Among these, CCR1, ERP27, HCK, MICB, and SLC2A3 showed protective associations with RAU risk, whereas CD177 and IFITM1 were positively associated with increased risk. Enrichment analysis revealed that these genes are involved in specific biological processes, including cell migration, immune response, and metabolic regulation, which are closely linked to RAU pathogenesis. CONCLUSION: This systematic study comprehensively investigates the critical causative genes underlying RAU, emphasizing the intricate relationships between immune regulation and metabolic disturbances in its pathology. These findings lay a solid foundation for the development of novel biomarkers and may inform future research on targeted therapeutic strategies for RAU.

Stomatitis, Aphthous↗

Inflammatory pathways and immune dysregulation in pediatric postoperative septic shock: A study integrating transcriptomics, machine learning and molecular docking.

This study elucidates the molecular and immune regulatory mechanisms of pediatric postoperative septic shock. Transcriptomic data were obtained from the Gene Expression Omnibus database. Differentially expressed genes were identified using the limma package, and gene co-expression modules were constructed using Weighted Gene Co-expression Network Analysis. Functional enrichment was performed via gene set enrichment analysis, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes analyses. Immune cell infiltration was assessed using ESTIMATE and CIBERSORT. Mendelian randomization was applied to explore causal relationships between gene expression and septic shock. Feature genes were selected using machine learning algorithms, and a diagnostic nomogram model was constructed. Finally, molecular docking analysis was performed to screen and evaluate the binding affinity of traditional Chinese medicine monomers to core target proteins. A total of 1331 differentially expressed genes were identified, and the turquoise module was strongly correlated with septic shock. Enrichment analysis revealed significant activation of IL-6/JAK/STAT3, TNF-&#x3b1;/NF-&#x3ba;B, and PI3K/Akt/mTOR pathways. Immune infiltration analysis indicated suppressed immune scores and imbalances in neutrophils, macrophages, T cells, and B cells. Mendelian randomization confirmed causal associations for 6 genes, including PIM3. The predictive model based on feature genes demonstrated high diagnostic performance. Molecular docking suggested that quercetin and astramembrannin I could stably bind PIM3. This study systematically identified core genes, dysregulated immune pathways, and candidate small-molecule interventions in pediatric septic shock, providing novel insights for early diagnosis and targeted therapy.

Humans↗

Integrative machine learning and transcriptomic analysis reveals molecular mechanisms underlying low survival rate in larval Chinese Bahaba (Bahaba taipingensis).

Chinese Bahaba (Bahaba taipingensis) is a Class I protected marine fish endemic to China. Low larvae survival during artificial breeding severely hinder population recovery. To investigate the molecular mechanism of high mortality in larval fish, this study performed RNA-seq on liver from naturally deceased (ND) and mass-dead (MD) individuals, combined with least absolute shrinkage and selection operator (LASSO) regression and random forest (RF) algorithms to screen for core signature genes. A total of 873 differentially expressed genes (DEGs) were identified, including 112 upregulated and 761 downregulated genes. GO and KEGG enrichment analyses revealed significant enrichment in amino acid metabolism disorders, one&#x2011;carbon folate pool impairment, PPAR signaling abnormalities, ECM-receptor interaction, focal adhesion pathway, indicating widespread metabolic suppression accompanied by extracellular matrix remodeling and signaling disturbances in the livers of MD fish. MAD pre-filtering combined with dual machine learning algorithms yielded 18 robust core signature genes, among which SLC38A4, MMP1, FADD, FKBP5, and APOB were consistently identified as high-frequency core genes by both algorithms. SLC38A4 exhibited the highest importance score in the RF model and was significantly downregulated, making it the primary molecule distinguishing ND from MD phenotypes. ROC curve analysis showed that both models achieved an AUC of 1.000 (95% CI lower bound: 0.610), confirming the precise discriminatory ability of the core genes. GSEA further demonstrated significant enrichment of this core gene set in ND samples. This study provides the first systematic elucidation of the molecular mechanisms underlying liver dysfunction in low survival rate B. taipingensis, characterized by amino acid transport impairment, metabolic reprogramming, and structural remodeling, offering theoretical foundations for health assessment, early mortality risk warning, and artificial breeding conservation of this species.

Animals↗