SEARCH · Search PubMed

Results for “transcriptomic”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

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

At least 343 records · Page 19Linked to original sources

PubMed2025

GSK3B inhibition partially reverses brain ethanol-induced transcriptomic changes in C57BL/6J mice: Expression network co-analysis with human genome-wide association studies.

Alcohol use disorder (AUD) is a chronic behavioral disease with greater than 50% of its risk due to complex genetic contributions. Existing pharmacological and behavioral treatments for AUD are minimally effective and underutilized. Animal model behavioral genetics and human genome-wide association studies have begun to identify individual genes contributing to the progressive compulsive consumption of ethanol that occurs with AUD, promising possible new therapeutic targets. Our laboratory has previously identified Gsk3b as a central member in a network of ethanol-responsive genes in mouse prefrontal cortex, which altered ethanol consumption with genetic manipulation and was also significantly associated with risk for alcohol dependence in human genome-wide association studies. Here we perform detailed brain RNA sequencing transcriptomic studies to characterize a highly specific and clinically available GSK3B pharmacological inhibitor, tideglusib, as a possible therapeutic for clinical trials on treatment of AUD. A model of chronic intermittent ethanol consumption was used to study gene expression changes in prefrontal cortex and nucleus accumbens in the presence or absence of tideglusib treatment. Multivariate analysis of differentially expressed genes showed that tideglusib largely reversed ethanol- induced expression changes for two prominent clusters of genes in both prefrontal cortex and nucleus accumbens. Bioinformatic analysis showed these genes to have prominent roles in neuronal functioning and synaptic activity. Additionally, mouse brain differential gene expression data was analyzed together with human protein-protein interaction and genome-wide association studies on AUD to derive networks responding to tideglusib and relevant to human genetic risk for alcohol dependence. These studies identified discrete networks significantly enriched with genes provisionally associated with AUD, and provide key information on central hubs of such networks. Together these studies document tideglusib as a major modulator of chronic ethanol consumption-evoked brain gene expression signatures, and identify possible new targets for therapeutic modulation of AUD.

Journal Article↗

PubMed2025

Multi-ancestry genome-wide and transcriptome-wide association analyses identified new risk loci and genes for inflammatory bowel disease.

To advance genetic understanding of inflammatory bowel disease (IBD), we conducted genome-wide association meta-analyses of 63,415 IBD cases of European and East Asian descendants and identified 90 previously unknown risk loci. Integrating multi-ancestry transcriptome-wide association studies (TWAS), cell type-specific TWAS, alternative splicing (AS-WAS), and alternative polyadenylation (APA-WAS) analyses using RNA-seq data from normal colon tissues of 707 European and 364 East Asian individuals, we uncovered 506 high-confidence IBD risk genes, including 384 not previously reported. These genes converge on immune regulation, microbial interaction, and other pathways central to IBD pathogenesis, with over half showing transcriptional dysregulation supported by single-cell and spatial omics analyses. Notably, 46 risk genes are targeted by 225 drugs that have been approved or in Phase II/III trials, including sulfasalazine already used in IBD therapy. Our study findings deepen the understanding of IBD genetics and support the development of precision medicine for its prevention and treatment.

GWAS↗

PubMed2026

Dual-contrastive learning for spatial domain identification in spatial transcriptomics with STAMGC.

Spatial transcriptomics (STs) have become a valuable approach for understanding the growth and development of organisms. Despite the recent emergence of numerous ST models, accurately identifying spatial domains remains challenging owing to the trade-off between preserving local details and reducing noise. Here, we introduce STAMGC, which is a dual-contrastive learning framework built upon graph convolutional networks. This model leverages regional and topological contrastive learning to jointly optimize the model, effectively reducing the noise in spatial domain identification and enhancing the extraction of detailed features. In this study, Gaussian smoothing, originally developed in the image processing field, is introduced to process ST data, providing a foundation for region-level contrastive learning by mitigating spatial discontinuities of gene expression signals. Experimental results indicate that STAMGC outperforms existing methods across multiple data sets according to comprehensive evaluations. Furthermore, STAMGC not only identifies finer structures in the mouse brain but also brings new discoveries for human breast cancer research.

Journal Article↗

PubMed2026

Comparative Transcriptomics Reveals Shared Downstream Pathways in Craniofacial Pathology.

Treacher Collins syndrome and Nager syndrome are craniofacial developmental disorders caused by defects in ribosome biogenesis and RNA splicing, respectively, yet they exhibit overlapping abnormalities affecting neural crest cell-derived craniofacial structures. To investigate shared downstream pathogenic mechanisms, we performed a comparative transcriptomic analysis of zebrafish polr1c and sf3b4 mutant models from our previous studies. Comparative analysis identified 17 shared differentially expressed genes (DEGs) between polr1c and sf3b4 mutants, with the majority of shared genes dysregulated in the same direction, indicating a coordinated rather than random transcriptional response. Gene ontology analysis identified ATP-dependent protein folding chaperone activity as the only shared molecular function, driven in part by upregulation of hsp90aa1.2, indicating a common proteostasis response. Because chaperone activity is linked to extracellular matrix (ECM) protein processing, we cross-referenced DEGs from both mutants against the curated zebrafish matrisome. Three of the 17 shared DEGs (serpinh1b, il11a, and lepa) were matrisome-associated and upregulated in both mutants. Serpinh1b, a collagen-specific chaperone, was strongly expressed in craniofacial cartilage and mesenchymal populations during pharyngeal arch development and exhibited nearly identical fold changes in both mutants. Il11a is of particular interest because its receptor, IL11RA, is known to be associated with human craniosynostosis, suggesting potential relevance to craniofacial development. Together, it is possible to hypothesize that shared chaperone-associated transcriptional changes, together with altered ECM-related gene expression, may contribute to polr1c- and sf3b4-associated craniofacial disorders, warranting further functional validation.

Extracellular Matrix↗

PubMed2004

Function-informed transcriptome analysis of Drosophila renal tubule.

BACKGROUND: Comprehensive, tissue-specific, microarray analysis is a potent tool for the identification of tightly defined expression patterns that might be missed in whole-organism scans. We applied such an analysis to Drosophila melanogaster Malpighian (renal) tubule, a defined differentiated tissue. RESULTS: The transcriptome of the D. melanogaster Malpighian tubule is highly reproducible and significantly different from that obtained from whole-organism arrays. More than 200 genes are more than 10-fold enriched and over 1,000 are significantly enriched. Of the top 200 genes, only 18 have previously been named, and only 45% have even estimates of function. In addition, 30 transcription factors, not previously implicated in tubule development, are shown to be enriched in adult tubule, and their expression patterns respect precisely the domains and cell types previously identified by enhancer trapping. Of Drosophila genes with close human disease homologs, 50 are enriched threefold or more, and eight enriched 10-fold or more, in tubule. Intriguingly, several of these diseases have human renal phenotypes, implying close conservation of renal function across 400 million years of divergent evolution. CONCLUSIONS: From those genes that are identifiable, a radically new view of the function of the tubule, emphasizing solute transport rather than fluid secretion, can be obtained. The results illustrate the phenotype gap: historically, the effort expended on a model organism has tended to concentrate on a relatively small set of processes, rather than on the spread of genes in the genome.

Animals↗

PubMed2025

Comparative genomics and transcriptomics of the Spiroplasma glossinidia strain sGff reveal insights into host interaction and trypanosome resistance in Glossina fuscipes fuscipes.

Tsetse (Glossina spp.) are vectors of African trypanosomes, the causative agents of Human and African Animal trypanosomiases, diseases that remain significant medical and socioeconomic challenges in sub-Saharan Africa. In addition to trypanosomes, tsetse harbor both obligate and facultative symbiotic bacteria that can influence vector competence and reproductive biology. One such facultative symbiont, Spiroplasma glossinidia, infects several tsetse species within the Palpalis subgroup. In Glossina fuscipes fuscipes (Gff), the Spiroplasma glossinidia strain sGff induces a trypanosome-refractory phenotype and negatively impacts reproductive fitness by reducing female fecundity. However, the mechanisms behind these Spiroplasma-derived phenotypes remain poorly understood. Here, we report successful in vitro cultivation of sGff and present complete genomes from three sources: in vitro cultured sGff and sGff isolated from both laboratory-maintained and wild-caught (Uganda) Gff flies. Comparative genomic analyses revealed a high degree of similarity in gene content and synteny among these sGff samples, confirming that they represent isolates of the same strain. Phylogenomic analyses placed sGff within the Spiroplasma poulsonii clade. The sGff genome is highly dynamic, containing numerous mobile genetic elements. Additionally, in silico annotations indicate that sGff relies on its host for both lipids and carbohydrates and produces several toxins, all of which could be implicated in the observed trypanosome refractory phenotype. Finally, comparative transcriptomic analysis of sGff from host hemolymph versus in vitro culture provided insights into potential factors relevant to host-symbiont interactions. Our findings provide a foundation for understanding the nutritional dialogue between sGff and its host and identify symbiotic products that may contribute to trypanosome resistance. Furthermore, the establishment of an in vitro culture system for sGff represents a significant resource for future functional studies with potential implications for vector control.

Glossina fuscipes fuscipes↗

PubMed2026

Genome-wide cis-expression Quantitative Trait Loci (eQTL) and transcriptomic signals reveal distinct molecular regulation across correlated feed efficiency traits.

INTRODUCTION: Feed efficiency (FE) is a complex trait which determines livestock production profitability, yet the molecular mechanisms behind it remain unclear. This study investigated the blood transcriptomic profile of lambs, alongside genotype data with the aim to uncover the genetic basis of FE traits such as absolute dry matter intake (DMIabsolute), DMI adjusted for body size (DMIadjusted), average daily live weight gain (ADG), and residual feed intake (RFI). MATERIALS AND METHODS: Bulk RNA-Seq and genotype data were analysed using three complementary approaches: differential gene expression (DGE) analysis, weighted gene co-expression network analysis (WGCNA), and cis-expression Quantitative Trait Loci (cis-eQTL) mapping. These methods were used independently to identify genes and regulatory networks associated with FE traits and to investigate evidence supporting multi-trait candidate gene selection. RESULTS: DGE analysis revealed 2, 24, 85 and 4 differentially expressed genes for DMIabsolute, DMIadjusted, ADG, and RFI (Padjusted < 0.05), functionally enriched in sensory perception, ATP-dependent chromatin remodeling, Notch signaling and immune response pathways. 9 gene modules significantly associated with the FE traits (P ≤ 0.05) with correlations ranging from r = -0.56 to 0.49, were identified using WGCNA. Single nucleotide polymorphism (SNP)-level cis-eQTL analysis identified 93 eSNPs associated with 74 genes (false discovery rate (FDR) < 0.05), while permutation-derived gene level analysis identified 280 eGenes (FDR < 0.2, empirical P < 0.03). Across the three analyses, applying thresholds of DGE (Padjusted < 0.05), WGCNA (correlation, P ≤ 0.05), and cis-eQTL gene-level significance (empirical P < 0.05), multiple overlapping genes were identified including DNMT3A, KANSL1, NCOR1 for DMIadjusted, ACOX2, FANCF, CIMIP2B, LOC101115106, ARMH2, LOC132657496 for ADG, and LOC114114576 for RFI representing regulators of variations in FE. DISCUSSION: The integration of DGE, WGCNA, and cis-eQTL analyses identified key genes and regulatory mechanisms associated with variation in FE traits. These results highlight that integrated multi-trait candidate gene identification approaches can reveal key genes that lower feed intake while maintaining animal growth, supporting breeding strategies aimed at improving efficiency and long-term economic sustainability in sheep.

average daily gain (ADG)↗

PubMed2026

From transcriptomic profiling to precision oncology: a bibliometric analysis of RNA sequencing in acute myeloid leukemia.

BACKGROUND: RNA sequencing (RNA-seq) has become an important tool for investigating the molecular heterogeneity of acute myeloid leukemia (AML); however, the global development and thematic evolution of this field remain inadequately characterized. OBJECTIVE: To map the global landscape of AML RNA-seq research and identify major knowledge domains, emerging themes, and temporal changes in research priorities. METHODS: Publications indexed in the Web of Science Core Collection and Scopus between January 1, 2007, and August 18, 2025, were retrieved. After database filtering, merging, and deduplication, 3,460 articles and reviews were included. CiteSpace, VOSviewer, the bibliometrix R package, and Microsoft Excel were used to analyze publication trends, collaboration networks, co-citation structures, keyword evolution, and citation bursts. RESULTS: Publication output increased steadily, accelerating after 2014. China contributed the largest number of publications (n = 547, 15.8%), whereas the United States had the highest total citation count. Major publication outlets spanned hematology, oncology, genomics, and molecular biology. Co-citation analysis identified prominent themes involving next-generation sequencing, gene mutations, KMT2A rearrangements, epigenetic dysregulation, leukemia-initiating cells, drug resistance, biomarkers, T-cell biology, and single-cell sequencing. Earlier literature emphasized sequencing technologies, gene expression profiling, and molecular alterations, whereas recent publications show increasing representation of cellular heterogeneity, single-cell transcriptomics, drug resistance, biomarker applications, immune-related research, and computational interpretation. CONCLUSION: While molecular characterization remains foundational, AML RNA-seq research has broadened to encompass increasingly prominent cellular, functional, computational, and translational dimensions. This study provides a structured overview of the field; nevertheless, bibliometric prominence should not be interpreted as direct evidence of clinical utility.

RNA sequencing↗

PubMed2026

Integrative cross-tissue transcriptome-wide association and metabolomic analysis reveals novel genetic risk loci for aortic aneurysm.

BACKGROUND: Aortic aneurysm (AA) is a life-threatening cardiovascular condition with a strong genetic component, however, its molecular mechanisms remain poorly understood. Although genome-wide association studies (GWAS) have identified numerous risk loci, most prior studies have investigated genetic and metabolic factors separately, leaving the causal pathways from genetic variants to disease largely unexplored. METHODS: We established an integrative framework combining cross-tissue transcriptome-wide association studies (TWAS) with metabolomic mediation analysis. First, we integrated GWAS data from FinnGen R12 with multi-tissue expression quantitative trait loci (eQTL) data from Genotype-Tissue Expression Project (GTEx) V8, then performed cross-tissue TWAS using the Unified Test for MOlecular SignaTures (UTMOST) and single-tissue validation with the Functional Summary-based Imputation (FUSION) to prioritize susceptibility genes. Second, we applied Mendelian randomization (MR), colocalization, and Fine-mapping Of CaUsal gene Sets (FOCUS) to assess causality and identify high-confidence genes. Third, we performed metabolite mediation analysis to uncover metabolic pathways linking genetic variants to disease risk. Finally, we validated key findings in mouse models of thoracic aortic aneurysm (TAA) and abdominal aortic aneurysm (AAA) using Quantitative Real-Time Reverse Transcription Polymerase Chain Reaction (RT-qPCR) and Western blotting. RESULTS: We identified multiple novel susceptibility genes for AA and its subtypes. Key genes included ADH family members (ADH1A, ADH1B, ADH4, ADH6) and ZNF827, which showed cross-subtype associations with strong colocalization evidence in vascular tissues. Metabolite mediation analysis revealed significant pathways involving N-acetylphenylalanine and methionine sulfoxide. Functional enrichment revealed distinct biological mechanisms: AA and AAA were primarily associated with metabolic pathways, whereas TAA-related genes were enriched in developmental and contractile processes. PheWAS indicated no significant off-target associations. Critically, experimental validation in mouse models confirmed significant upregulation of ZNF827 in TAA and ADH6 in AAA at both mRNA and protein levels, corroborating the genetic predictions. CONCLUSION: This integrated cross-omics analysis identifies novel genetic loci and, crucially, uncovers specific nutrient-related metabolic pathways that mediate genetic risk. These findings provide a mechanistic basis for future nutritional and metabolic intervention studies in AA and its subtypes.

MAGMA↗

PubMed2026

A minimal three-arm oral regimen for healthspan: mechanistic alignment with transcriptomic signals from a large parental-lifespan GWAS.

A large genome-wide association study of parental lifespan was reported in 2019. A later transcriptome-wide association study (TWAS) based on those summary statistics identified a set of transcriptional programs associated with longer genetically predicted survival, including increased brain NAD + salvage, especially NMNAT2, reduced glucose-stimulated insulin secretion, a shift toward synaptic pruning with less broad plasticity, and a glial pattern characterized by relatively greater microglial and lower astrocytic signatures, with only weak pan-tissue senescence signals. Building on those directional findings, this short communication proposes a minimal three-arm oral regimen with unequal evidentiary weight: first, the Cheung Glutamatergic Regimen, consisting of low-dose dextromethorphan potentiated by a CYP2D6 inhibitor together with piracetam and L-glutamine, as an exploratory adjunct aimed at preserving residual functional connectivity; second, daily nicotinamide mononucleotide and N-acetylcysteine with pulsed senolytics for NAD + salvage and senescence modulation; and third, GLP-1 receptor agonism for metabolic reprogramming. The NAD+/senescence arm is the primary mechanistic anchor, GLP-1 receptor agonism provides secondary metabolic support, and the glutamatergic arm is exploratory. Each arm targets a separate node within the pruning-plasticity-metabolic triad. The regimen is fully oral, uses conservative dosing, and draws on prior therapeutic or human-exposure data, although the proposed combination has no established safety profile. Although direct combination data are lacking and the foundational TWAS remains a preprint, the components show plausible but uneven mechanistic alignment with the TWAS signals and may justify carefully designed, safety-focused pilot evaluation.

GLP-1↗

PubMed2026

Transcriptomic Profiling of Canine Testicular Leydig Cell Tumors Uncovers Key Upregulated Gene Pathways.

Total RNA was isolated from sections of healthy testes and Leydig cell tumors of mixed-breed dogs using TMA Master II device. The RNA-seq libraries were sequenced on the Illumina platform. Following differential expression analysis, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Set Enrichment Analysis (GSEA) were applied with quality control obtained using FastQC and Trimmomatic. This analysis revealed 1500 transcripts, including 928 upregulated and 168 downregulated genes. The results demonstrated that a significant proportion of these differentially expressed genes are directly involved in the control of sex steroid production (CYP11A1, STAR, and 3β-HSD3B1) or tube formation, angiogenesis, and extracellular matrix remodeling in interstitial cells (ESM1, FGG, and VEGFA). Moreover, we identified the upregulation of transcripts responsible for neurotransmitter or neuroendocrine signaling (SLC6A4, GRIN2C, GABRB3) and cholesterol metabolism and its regulation (GPX3, MSMO1, DHCR24). These genes were strongly associated with the phosphatidylinositol-3-kinase (PI3K)-Protein Kinase B (Akt) cascade and extracellular matrix interactions, features shared with various malignancies. Alterations in estrogen and relaxin signaling appear to be distinctive, understudied mechanisms specific to canine Leydig cell tumors. Concurrently, downregulated genes (e.g., DMRTC2, SEMA3C, ALOX12) were linked with cell differentiation, signaling and immunoregulatory pathway suppression involved in tumorigenesis. A complex transcriptomic profile of canine Leydig cell tumors was developed, revealing a conserved oncogenic core shared in some aspects with human malignancies alongside unique species-specific alterations. Findings seem to be useful for identifying novel diagnostic biomarkers and targeted therapies in veterinary oncology, establishing canine reproductive tissues as a valuable comparative biomedical model for research in human.

Leydig cell tumor↗

PubMed2026

Integrated Metabolomic and Transcriptomic Analysis Reveals Tissue-Specific Secondary Metabolic Differentiation and Indole Alkaloid Accumulation in Evodia rutaecarpa.

Evodia rutaecarpa is a valuable medicinal plant, yet its non-medicinal tissues remain largely underexplored. Here, we integrated ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS)-based widely targeted metabolomics and RNA sequencing (RNA-seq) transcriptomics to systematically profile the metabolic and transcriptional landscapes of roots, stems, leaves, and flowers of Evodia rutaecarpa (Juss.) Benth. Our aim was to characterize tissue-specific metabolic differentiation and its underlying transcriptional regulatory mechanisms. Metabolomic analysis, employing principal component analysis (PCA) and orthogonal partial least squares-discriminant analysis (OPLS-DA) with robust model parameters (R2Y > 0.9, Q2 > 0.5), identified 3090 differential metabolite features (variable importance in projection, VIP > 1.0; p < 0.05) across the four tissues, which exhibited distinct tissue-specific clustering patterns. Integrated Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis and weighted gene co-expression network analysis (WGCNA) revealed that roots specifically accumulated quinolone alkaloids and flavonoid glycosides, accompanied by the coordinated upregulation of genes involved in flavonoid and phenylpropanoid biosynthetic pathways. In contrast, stems, leaves, and flowers were enriched in indole alkaloids (evodiamine and rutaecarpine) and volatile oil precursors, with concurrent upregulation of genes involved in tryptophan metabolism and indole alkaloid biosynthesis (e.g., tryptophan decarboxylase, TDC; s N-methyltransferase, NMT). Notably, leaves and flowers displayed particularly high accumulation levels of these bioactive alkaloids, suggesting their potential as alternative sources for industrial and pharmaceutical applications. WGCNA further identified multiple transcription factors and structural gene modules tightly correlated with evodiamine accumulation, offering promising candidate regulators for future biosynthetic pathway engineering. Collectively, this multi-omics integration study systematically elucidates the tissue-partitioned secondary metabolism of Evodia rutaecarpa (Juss.) Benth. and provides a solid scientific foundation for full-plant resource utilization, targeted development of non-medicinal tissues, and future metabolic engineering of indole alkaloid production.

Evodia rutaecarpa↗

PubMed2026

Integrative transcriptomic, spatial and functional-genomic analysis identifies a UFMylation-related vascular-stromal program and prioritizes WWTR1 in glioblastoma.

Glioblastoma (GBM) contains spatially organized stress-adaptive and vascular niches. Because transcript abundance does not measure UFM1 conjugation, we asked whether a UFMylation-related transcriptional axis identifies a reproducible tissue program and alters candidate prioritization. In 518 unique primary TCGA-GBM tumors profiled on the Affymetrix HT Human Genome U133A array, weighted gene co-expression network analysis of 8,000 variable genes yielded 12 modules. The 278-gene green module ranked first across nine prespecified traits (mean |r|=0.637). Direct overlap comprised 1/3 measurable UFMylation-core, 5/19 ER-stress/UPR, and 2/15 proteostasis genes; after excluding overlapping genes, correlations with the green eigengene remained significant (r = 0.373, 0.831, 0.639, and 0.699 for UFMylation-core, ER-stress/UPR, proteostasis, and composite scores, respectively). The green score was associated with overall survival per standard-deviation increase (HR 1.17, 95% CI 1.07-1.28), although clinical adjustment attenuated the estimate. In a 10-sample single-cell dataset, sample-level scores were higher in pericytes and endothelial cells than in malignant cells. Donor-aware IvyGAP analysis supported regional organization, whereas one Visium section showed stronger concordance with ER-stress/UPR and mesenchymal scores than with the UFMylation-core score. CellChat indicated pathway-selective rather than global remodeling of inferred vascular communication. Layer ablation moved WWTR1 from rank 48 using WGCNA alone to rank 4 overall and rank 1 among non-common-essential genes after cross-platform integration. These findings define an ER-stress/mesenchymal-weighted, UFMylation-related vascular-stromal transcriptional association and nominate WWTR1 for experimental testing.

Humans↗

PubMed2026

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↗

PubMed2026

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‑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↗

PubMed2023

GBMdeconvoluteR accurately infers proportions of neoplastic and immune cell populations from bulk glioblastoma transcriptomics data.

BACKGROUND: Characterizing and quantifying cell types within glioblastoma (GBM) tumors at scale will facilitate a better understanding of the association between the cellular landscape and tumor phenotypes or clinical correlates. We aimed to develop a tool that deconvolutes immune and neoplastic cells within the GBM tumor microenvironment from bulk RNA sequencing data. METHODS: We developed an IDH wild-type (IDHwt) GBM-specific single immune cell reference consisting of B cells, T-cells, NK-cells, microglia, tumor associated macrophages, monocytes, mast and DC cells. We used this alongside an existing neoplastic single cell-type reference for astrocyte-like, oligodendrocyte- and neuronal progenitor-like and mesenchymal GBM cancer cells to create both marker and gene signature matrix-based deconvolution tools. We applied single-cell resolution imaging mass cytometry (IMC) to ten IDHwt GBM samples, five paired primary and recurrent tumors, to determine which deconvolution approach performed best. RESULTS: Marker-based deconvolution using GBM-tissue specific markers was most accurate for both immune cells and cancer cells, so we packaged this approach as GBMdeconvoluteR. We applied GBMdeconvoluteR to bulk GBM RNAseq data from The Cancer Genome Atlas and recapitulated recent findings from multi-omics single cell studies with regards associations between mesenchymal GBM cancer cells and both lymphoid and myeloid cells. Furthermore, we expanded upon this to show that these associations are stronger in patients with worse prognosis. CONCLUSIONS: GBMdeconvoluteR accurately quantifies immune and neoplastic cell proportions in IDHwt GBM bulk RNA sequencing data and is accessible here: https://gbmdeconvoluter.leeds.ac.uk.

Humans↗

PubMed2026

Transcriptome Analysis Reveals Key Drought-Stress-Responsive Genes in Two Bermudagrass Genotypes.

Drought inhibits grass development and survival. However, molecular-based studies on drought tolerance mechanisms in bermudagrass (Cynodon dactylon) remain scarce. Therefore, a drought-resistant bermudagrass (Tianshui) and a drought-sensitive (Zhengzhou) genotype were selected and subjected to 28 days of 50% drought stress. Leaves were sampled for RNA sequencing. Under drought stress, 2410 differentially expressed genes (DEGs) were discovered in which 1214 were upregulated (tolerant vs. sensitive) and 1196 downregulated. Kyoto Encyclopedia of Genes and Genomes (KEGG) indicated that these specific DEGs are notably present in hormonal signal transduction pathways, flavonoids biogenesis, carbohydrate metabolic processes, abscisic acid-mediated pathways, MAPK signaling, and gluconeogenesis. Additionally, the plant hormone signal transduction pathway is predominantly linked to abscisic acid signal transduction, and many other plant hormones were also drought-responsive. The study specifically targeted genes associated with the antioxidant enzyme system, with a particular emphasis on responsive TFs such as MYB, bHLH, bZIP, GRAS, and WRKY. This study establishes the theoretical framework and identifies gene sources for the genetic enhancement and breeding of bermudagrass in the future.

Cynodon↗

PubMed2026

Long-read transcriptomics corrects Trichomonas vaginalis intron annotations and refines transcript-end features.

BACKGROUND: Trichomonas vaginalis causes the most prevalent non-viral sexually transmitted infection worldwide. Despite its large genome (181.5 Mb; 36,310 predicted protein-coding genes in NYU_TvagG3_2), intron annotations remain limited and inconsistently validated. A recent short-read RNA-seq study reported 63 putative active introns, but short reads can misassign splice boundaries and cannot resolve complete transcript structures. METHODS: We integrated Oxford Nanopore direct RNA sequencing (DRS), ONT cDNA long-read sequencing, and Illumina RNA-seq to refine intron annotations, transcript-end features, and UTR boundaries in T. vaginalis. Candidate introns were validated by targeted PCR and Sanger sequencing, and representative splicing events were further assessed using public SRA datasets. RESULTS: Starting from 31 historically annotated introns, motif-guided long-read screening and orthogonal validation identified 17 additional validated introns, increasing the curated set to 48 confirmed introns. Among these 17 events, three were previously unrecognized in the current NYU_TvagG3_2 reference annotation. We also corrected five reported loci, including two false-positive introns, two splice-coordinate misannotations, and one gene-sequence error. DRS further supported transcript termination site mapping, UAAA polyadenylation-signal profiling relative to poly(A) addition sites, and single-molecule poly(A)-tail estimation. StringTie mixed-mode assemblies provided updated UTR boundaries for intron-bearing transcripts and transcripts without curated introns. CONCLUSIONS: This study provides a rigorously validated, long-read-refined resource of intron annotations, UTR boundaries, and UAAA-guided transcript-end features for T. vaginalis, together with a reproducible workflow for non-model protists. These refinements improve the current reference annotation and support future studies of functional genomics, parasite biology, pathogenesis, and diagnostic development.

Trichomonas vaginalis↗