Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “RNA biomarkers”

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.

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

383 records · Page 2Linked to original sources
PubMed2026

Epigenetics and In Silico Transcriptome Analysis of Pediatric Acute Myeloid Leukemia.

Pediatric acute myeloid leukemia (AML) is a heterogeneous hematologic malignancy that accounts for about 15%-20% of childhood leukemias. Despite therapeutic advances, relapses remain common, and survival for high-risk patients is below 60%. Unlike adult AML, pediatric AML displays distinct genetic mutations, including FLT3-ITD, NPM1, KMT2A rearrangements, and core-binding factors (CBF) fusions, as well as extensive epigenetic dysregulation. Aberrant DNA methylation, histone modifications, and altered non-coding RNA expressions disrupt hematopoietic differentiation and activate oncogenic transcriptional networks. Recent advances in silico transcriptomic analysis have transformed the study of pediatric AML by integrating gene expression and epigenetic data to identify molecular drivers and regulatory networks. Computational RNA-seq pipelines and pathway analyses have highlighted key epigenetic regulators, including DNMT3A, TET2, and HDACs, as potential therapeutic targets. Multi-omics approaches combining transcriptomic, methylomic, and chromatin accessibility data are increasingly used to define biomarkers for diagnosis, prognosis, and therapeutic response. This review provides a comprehensive overview of the molecular and epigenetic landscape of pediatric AML, emphasizing the power of in silico transcriptome analysis to uncover disease mechanisms, refine patient stratification, and guide the development of precision-based epigenetic therapies aimed at improving long-term outcomes in children with AML.

Humans↗
PubMed2026

Human endogenous retroviruses leading to autoimmune diseases.

Human endogenous retroviruses (HERVs) comprise approximately 8% of the human genome and were long regarded as inert remnants of ancestral retroviral infections. Increasing evidence indicates that HERVs are active genomic elements capable of influencing transcriptional programs, modulating immune responses, and contributing to disease pathogenesis. Under physiological conditions, HERV expression is tightly controlled by epigenetic mechanisms; however, infections, chronic inflammation, aging, and diverse environmental stimuli can promote HERV reactivation. HERV-derived RNAs and proteins engage innate immune sensors and trigger antiviral-like responses through mechanisms of viral mimicry, leading to activation of type I interferon and other inflammatory pathways. HERV dysregulation has been associated with disease-relevant immune pathways. This review summarizes recent advances linking HERVs to autoimmune disease pathogenesis and discusses their potential translational relevance as biomarkers and therapeutic targets.

Humans↗
PubMed2026

A machine learning-derived and functionally validated circadian rhythm signature predicts clinical outcomes and in silico drug sensitivity in colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) displays considerable heterogeneity in clinical outcomes, highlighting the need for reliable prognostic biomarkers. While the aberrant expression of circadian rhythm-related genes has been implicated in cancer pathogenesis, its comprehensive role in CRC progression and predicted therapeutic vulnerabilities remains inadequately characterized. METHODS: Bulk and single-cell RNA-sequencing data were integrated from multiple CRC cohorts. A circadian rhythm signature (CRS) was developed through machine learning algorithms and validated for prognostic value. Comprehensive analyses of tumor microenvironment, genomic alterations, and drug sensitivity were performed. Furthermore, the biological function of the core gene, BHLHE40, was validated in CRC cell lines through CCK-8, EdU, and wound healing assays. RESULTS: Single-cell analysis demonstrated an elevated expression signature of circadian rhythm-related genes in dendritic cells. The optimized CRS, comprising 14 circadian rhythm-related genes, successfully categorized patients into high- and low-risk groups. Patients with a high CRS showed markedly poorer overall survival and computationally inferred immunosuppressive features, including reduced CD8+ T cell infiltration and increased M2 macrophage polarization. Genomic analysis revealed enhanced mutation burden in TP53 and alterations in RTK-RAS/WNT pathways. Notably, in vitro assays confirmed that BHLHE40 is significantly overexpressed in CRC cells. Knockdown of BHLHE40 markedly inhibited tumor cell proliferation and migration. Drug sensitivity profiling identified bexarotene and SMER-3 as potential therapeutic options for high-CRS patients. A nomogram integrating CRS with clinical parameters demonstrated superior predictive accuracy for 1-, 3-, and 5-year survival. CONCLUSIONS: The CRS represents a promising prognostic biomarker that reflects tumor immune status and genomic features, providing valuable insights for personalized treatment strategies in CRC.

Circadian rhythm↗
PubMed2026

Genome-Wide Impact of Human DBR1 Depletion on RNA Processing Networks Reveal a Connection Between Pre-mRNA Splicing, mRNA Surveillance and Stress Granule Dynamics.

The RNA lariat debranching enzyme DBR1 is essential for intron turnover and RNA metabolism, yet its broader impact on transcriptome regulation remains incompletely defined. To elucidate the consequences of DBR1 depletion, we performed transcriptome-wide RNA sequencing of DBR1-knockdown and wild-type HEK293 cells. Differential expression analysis revealed widespread perturbations in pathways linked to RNA splicing, mRNA surveillance, translational control, and stress-granule biology. Many of the most significantly altered transcripts encode splicing factors and RNA quality-control components, underscoring DBR1's influence on post-transcriptional regulation. Alternative splicing analysis showed changes across multiple event types, with exon skipping accounting for >50% of events, followed by mutually exclusive exons, alternative 5' and 3' splice sites, and retained introns, indicating that DBR1 depletion induces pervasive splicing defects. Direct spliceosome inhibition using isoginkgetin (blocks tri-snRNP recruitment) and pladienolide B (targets SF3B1) reproduced the DBR1-KD mis-splicing patterns of cell signaling genes and factors involved in RNA metabolism, supporting a functional link between DBR1 activity and alternative splicing. Notably, DBR1 knockdown revealed a subset of transcripts that are both NMD-sensitive and enriched within stress granules. Consistent with this observation, G3BP1 immunopurification and confocal microscopy further support a role for DBR1 and UPF1 in stress-granule dynamics, suggesting that these factors may participate at distinct stages to influence mRNA fate under stress conditions. Together, these findings indicate that DBR1 functions beyond lariat RNA turnover as a common regulator of RNA processing, transcriptome stability, and stress granule homeostasis, revealing intricate crosstalk between RNA splicing and RNA quality control pathways in human cells.

Humans↗
PubMed2026

Whole-transcriptome RNA sequencing and ceRNA network analyses provide novel insights into the antibacterial immune response of Hippocampus abdominalis against Vibrio harveyi.

Long non-coding RNAs (lncRNAs) stand as newly-arisen molecular types that exert regulatory effects, able to operate as competitive endogenous RNAs (ceRNAs) to engage microRNAs (miRNAs) in interaction, resulting in the recovery of target mRNA expression and activity. Increasing evidences indicate that the ceRNA network affects various biological processes in mammals, including development, cellular differentiation, metabolism, immune response, and disease pathogenesis. In teleost fish, the lncRNA-miRNA-mRNA regulatory networks have been reported occasionally. However, up to now, the roles of lncRNAs in the big-belly seahorse (Hippocampus abdominalis) remains unclear. In this study, we reported for the first time, via whole-transcriptome RNA sequencing, the lncRNA mediated ceRNA regulatory network in Vibrio harveyi-infected H. abdominalis. A total of 4197 differentially expressed mRNAs (DE-mRNAs), 1317 DE-lncRNAs, and 183 DE-miRNAs were identified. Furthermore, the crosstalk between miRNAs and lncRNAs as well as between miRNAs and mRNAs was inferred based on the negative correlations between miRNAs and their target lncRNAs/mRNAs. A core immune associated lncRNA-miRNA-mRNA putative regulatory network was thus constructed, comprising 211 lncRNA-miRNA and 224 mRNA-miRNA pairs. In conclusion, our findings provide an integrative overview of the ceRNA regulatory networks on the underlying immune responses to V. harveyi infection in the big-belly seahorse, and offer a solid theoretical foundation for the comparative immunological research of teleost fish.

Animals↗
PubMed2026

Circular RNAs in amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder characterized by the progressive loss of motor neurons, with most cases lacking a clear genetic basis. Emerging evidence highlights the involvement of non-coding RNAs, particularly circular RNAs (circRNAs), in disease onset and progression. Here, we investigated circRNAs implicated in ALS and related motor neuron diseases (MNDs). Here, we provide a general overview of circular RNA metabolism and cellular functions. We then present our systematic literature review that identified ALS-associated circRNAs, followed by in silico analyses of 15 circular RNA candidates that were selected based on the most compelling data regarding ALS. Our results revealed that several circular RNAs regulate ALS-related genes, such as unfolded protein response, oxidative stress, cell cycle regulation, and apoptosis. Protein-RNA interaction analysis further showed that ALS-related circRNAs can sponge 20 RNA-binding proteins. Additionally, molecular docking analysis demonstrated that ALS-associated FUS variants significantly alter its binding affinity to circular RNAs. RNA-seq data from ALS patients confirmed significant alterations in the expression of host genes of ALS-related circRNAs and hub proteins in ALS-affected CNS tissues. Collectively, our findings identify circRNAs as potential key contributors to ALS pathogenesis.

Amyotrophic Lateral Sclerosis↗
PubMed2026

RNA dysregulation as a determinant of aging and neurodegenerative vulnerability.

In the nervous system, aging causes deterioration of cellular and molecular processes that are associated with declines in cognition, sensory perception, and motor coordination. Aging is also the strongest risk factor for neurodegenerative disease, yet the mechanisms by which aging predisposes neurons to dysfunction remain incompletely understood. While genomic instability, proteostasis decline, mitochondrial dysfunction, and chronic inflammation have dominated prevailing models, recent evidence highlights RNA dysregulation as a central component of age-associated decline. In this review, we summarize recent findings suggesting that aging progressively erodes RNA regulatory fidelity through alterations in RNA-binding protein abundance, localization, biophysical behavior, and RNA interactions. We argue that age-dependent RNA dysregulation represents an important mechanism that converges with genetic risk to drive neuronal vulnerability and neurodegeneration.

RNA dysregulation↗
PubMed2026

A contextual activity score (CAS) for inferring ADAR-associated transcriptional activity across RNA-seq, single-cell, and spatial transcriptomics.

BACKGROUND AND OBJECTIVE: Adenosine-to-inosine RNA editing, catalyzed by Adenosine Deaminases Acting on RNA (ADARs), is a widespread modification involved in neural function, immune regulation, and cancer. The Alu Editing Index (AEI) is the standard metric to estimate ADAR activity but requires raw sequencing reads and is poorly suited for single-cell and spatial transcriptomic data. This study aimed to develop an alternative framework for inferring ADAR-associated transcriptional activity from gene expression data across diverse transcriptomic technologies. METHODS: We developed the Contextual Activity Score (CAS), a framework based on transcriptional signatures from ADAR perturbation experiments. Context-specific signatures were generated for human neurons, mouse neurons, and cancer models to infer ADAR1 and ADAR2 activity. CAS was computed from normalized gene expression matrices using regulon-based enrichment analysis. Performance was evaluated by comparing with the Alu Editing Index across bulk RNA sequencing datasets, simulated sequencing depths, and library preparation protocols. RESULTS: CAS showed strong concordance with the Alu Editing Index across multiple datasets, while remaining robust to reduced sequencing depth and different library protocols. Unlike the Alu Editing Index, CAS can be applied to single-cell and spatial transcriptomic data and enables the independent assessment of ADAR2 activity. In cancer and neuronal contexts, CAS captured biologically meaningful variations in ADAR-associated transcriptional activity at sample, cell-type, and spatial levels. CONCLUSION: CAS provides a scalable approach applicable across multiple RNA-seq protocols for estimating ADAR-associated transcriptional activity using gene expression data. This method, implemented in an open-source R package for broad adoption, expands the ability to study ADAR-associated transcriptional activity across transcriptomic modalities where direct editing quantification is challenging, such as single-cell and spatial transcriptomics.

Adenosine Deaminase↗
PubMed2026

O'nyong-nyong virus adaptive mutations in non-structural protein 1 and 3 enhance RNA replication and overcome FHL1 requirement.

Arthritogenic alphaviruses, like o'nyong-nyong virus (ONNV), cause debilitating musculoskeletal diseases and are geographically expanding. To predict their emergence, we seek to better understand evolutionary mechanisms that enable changes in virus tropism. Here, we identify adaptive mutations in the ONNV non-structural proteins (nsPs) that arose during cellular serial passaging and enabled ONNV to infect non-permissive Lunet cells. Using shotgun proteomics, we show that this human hepatoma cell line lacks the four-and-a-half-LIM domain protein 1 (FHL1), an essential host factor in ONNV RNA replication. Individual single nucleotide mutations in the nsP1 ring-aperture membrane-binding and oligomerization domain, the nsP3 macrodomain, and the nsP3 opal stop codon overcome FHL1 deficiency in Lunet cells by enhanced RNA replication. These findings demonstrate how subtle genomic changes in nsPs can profoundly influence alphavirus replication and tropism.

LIM Domain Proteins↗
PubMed2026

Comprehensive analysis of mRNA-microRNA-lncRNA expression profiles in post-traumatic elbow heterotopic ossification using RNA sequencing and experimental validation.

BACKGROUND: This study aimed to profile the molecular signatures of post-traumatic elbow heterotopic ossification (HO) to identify key regulators and potential therapeutic targets. METHODS: Total RNA from post-traumatic elbow HO tissues (n=4) and normal bone tissues (n=6) was subjected to high-throughput sequencing to identify differentially expressed mRNAs (DEGs), microRNAs (DEMs), and lncRNAs (DELs). Bioinformatics analyses included Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, protein-protein interaction network construction, and transcription factor (TF)-microRNA-mRNA network analysis. The expression trends of four most upregulated and four most downregulated DEGs were validated by real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR). RESULTS: We identified 2,138 DEGs, 40 DEMs, and 905 DELs. DEGs were significantly enriched in biological process "bone mineralization," cellular component "plasma membrane," molecular function "integrin binding," and pathways including PI3K-Akt, NF-κB, JAK-STAT, and TNF signaling pathways. Hub genes with high connectivity included MMP9, IL6, MMP3, CTSK, and BGLAP. Integrated network analysis highlighted the transcription factor JUN and key microRNAs (hsa-miR-124-3p, hsa-miR-548c-3p, and hsa-miR-135b). The qRT-PCR results confirmed the expression trends of selected DEGs. CONCLUSIONS: This study, for the first time, profiled the differentially expressed mRNAs, microRNAs, and lncRNAs in post-traumatic elbow HO using high-throughput RNA sequencing. These findings provide valuable insights into the molecular mechanisms of HO following elbow trauma. The identified hub genes (MMP9, IL6, MMP3, CTSK, and BGLAP), key TF (JUN), and key microRNAs (hsa-miR-124-3p, hsa-miR-548c-3p, and hsa-miR-135b) may serve as potential therapeutic targets for preventing and treating post-traumatic elbow HO.

Humans↗
PubMed2026

Stage-specific ROMO1 in rheumatoid arthritis: predictive immune insights into the MIF pathway and HLA-DR/IL2RA axis via integrated GWAS, transcriptomic, single-cell, and spatial profiling.

Emerging evidence links reactive oxygen species modulator 1 (ROMO1), a key mitochondrial ROS regulator, to rheumatoid arthritis (RA) pathogenesis. However, its exact mechanism remains elusive given the conflicting evidence about its specific function. We used a four-level integrative framework combining multi-omics data and literature‑supported mechanistic inference. At the genetic level, Mendelian randomization (MR) was performed to explore potential causal relationships between ROMO1, IL2RA, HLA-DR, MIF, and RA risk, followed by differential expression analysis and machine learning-based feature selection to identify key mROS genes. The temporal expression dynamics of ROMO1 were assessed in RA progression. At the cellular and tissue levels, we integrated single-cell RNA sequencing and spatial transcriptomics to map cell-type-specific expression and synovial localization of ROMO1-related immune cells and pathways. Finally, our multi-omics findings were contextualized with literature-supported mechanistic inference. (1) MR results were consistent with a potential protective effect of ROMO1 on RA (OR = 0.52) and its potential regulation of risk factors IL2RA (OR = 0.46) and HLA-DR (OR = 0.40). Conversely, IL2RA (OR = 1.42), HLA-DR (OR = 1.88), and MIF (OR = 1.17) were positively associated with RA risk. Additionally, ROMO1 was identified as a top candidate diagnostic predictor with stage-specific dynamics: downregulated in the early but upregulated in the late/remission stages. (2) Single-cell RNA sequencing showed ROMO1's cell-specific expression in CD14+ HLA-DR+ CD74+ monocytes and CD4+ IL2RA+ T cells. Cell communication analysis further suggested that these cells may participate in MIF pathway regulation. Spatial transcriptomics subsequently identified that ROMO1-related cells localized to synovial pathological regions, with MIF pathway changes correlated with RA progression. (3) Finally, literature-supported mechanistic inference suggests that ROMO1 may modulate mROS levels to promote anti-inflammatory M2 macrophage polarization, which could theoretically contribute to reduced systemic inflammation and the alleviation of multi-organ decline in RA. This integrated multi-omics investigation, supported by literature-based mechanistic inference, suggests ROMO1 as a stage-dependent biomarker candidate and potential immune regulator in RA.

Humans↗
PubMed2026

CNOT1 is a potential YTHDF2 target that orchestrates maternal mRNA decay and zygotic genome activation during goat embryogenesis.

Timely and efficient degradation of maternal mRNA is essential for early embryonic development, which occurs from fertilization through the initiation of zygotic genome activation (ZGA). Yet, the regulatory mechanisms governing this process remain poorly characterized. In the present study, we investigated the function of CCR4-NOT transcription complex subunit 1 (CNOT1) during goat embryogenesis. We found that CNOT1 was upregulated during mammalian ZGA, and that its knockdown led to developmental arrest and a marked reduction in blastocyst formation. Moreover, CNOT1 knockdown impaired nascent RNA activity, resulting in 814 upregulated and 1014 downregulated genes, which were enriched for RNA splicing, regulation of chromosome organization, and RNA localization. RNA splicing analysis revealed differential splicing events in 2959 genes, of which 259 were downregulated following CNOT1 knockdown. Notably, CNOT1 was predicted to crosstalk with the m6A reader YTHDF2. Knockdown of YTHDF2 resulted in CNOT1 downregulation at the 8-cell stage in goats and increased transcription levels around polyadenylation sites during ZGA in mice. Together, these findings indicate that CNOT1 is a potential YTHDF2 target that orchestrates maternal mRNA decay and ZGA during goat embryogenesis. Our work provides new insight into the complex regulatory landscape underlying ZGA and may inform strategies to improve the efficiency of goat embryogenesis.

Animals↗
PubMed2026

Ramu stunt virus genome reveals previously unreported segments and nucleocapsid domain duplication in Mechlorovirus.

Ramu stunt virus (RmSV), a member of the genus Mechlorovirus within the family Phenuiviridae, was previously described as a six-segmented RNA virus infecting sugarcane. In this study, we re-examined type material and additional isolates using high-throughput sequencing and RT-PCR validation, revealing that RmSV possesses a nine-segmented genome, making it the largest reported in the Phenuiviridae. This expanded architecture includes duplicated RNA segments (RNA 2a and RNA 2b) encoding nucleocapsid-like proteins and two novel segments (RNA 7 and RNA 8). Comparative analysis showed that RNA 2a and 2b share about 84% amino acid identity, while RNA 5 encodes a third nucleocapsid homolog, indicating unprecedented domain redundancy. Structural modeling confirmed that all three nucleocapsid proteins maintain a conserved fold despite low sequence identity, with electrostatic mapping suggesting differential RNA-binding potential. Additionally, RNA 6 encodes a hypothetical protein structurally similar to the rice stripe virus disease-specific S-protein, implicating a role in symptom development. Transcript abundance analysis revealed RNA 6 as the most highly expressed segment across isolates. These findings revise the genomic composition of RmSV, highlight mechanisms of genome plasticity and adaptive evolution in plant-infecting bunyaviruses, and underscore practical implications for diagnostic assay design, resistance breeding, and biosecurity surveillance.

Genome, Viral↗
PubMed2026

Comparative transcriptome analysis reveals ncRNA-mediated regulatory networks associated with muscle crispiness in grass carp.

Non-coding RNAs (ncRNAs) have been demonstrated to be involved in muscle development and to function as key regulators. However, the molecular mechanism underlying muscle crispiness in grass carp (GC) remains poorly understood, and whether these ncRNAs are involved in its regulation is still unknown. In the current investigation, differentially expressed (DE) RNAs (including lncRNAs, circRNAs, miRNAs, and mRNAs) were identified; concomitantly, target genes prediction was conducted, and functional and signaling pathway enrichment analyses were performed. Pathways related to muscle crispiness were identified, and the competitive endogenous RNA (ceRNA) (lncRNA/circRNA-miRNA-mRNA) regulatory network was further constructed. The results showed that a total of 126 DE-lncRNAs, 17 DE-circRNAs, 329 DE-miRNAs, and 442 DE-mRNAs were identified in muscle tissues of both the GC and crisp grass carp (CGC). GO and KEGG enrichment analyses revealed that target genes of DE-ncRNAs were significantly enriched in signaling pathways, including structural constituents of muscle, apoptosis, oxidative phosphorylation, and regulation of actin cytoskeleton, suggesting that these pathways may be involved in muscle texture remodeling. Subsequently, DE-RNAs enriched in related pathways were identified, and a core ceRNA regulation network comprising 3 lncRNAs, 4 circRNAs, 3 miRNAs, and 17 mRNAs was constructed. Additionally, 10 DE-RNAs from randomly selected groups were validated by qRT-PCR. Our findings not only provide scientific evidence elucidating the molecular mechanisms underlying muscle crispiness in GC but also establish a foundation for studying changes in muscle textural qualities across other fish species.

Animals↗
PubMed2026

Meta-PseU: A meta-classifier for robust prediction of RNA pseudouridine modification sites from long sequences.

BACKGROUND AND OBJECTIVES: Pseudouridine (Ψ) represents one of the most abundant and conserved RNA modifications. Ψ provides an additional hydrogen-bond donor that enhances RNA structural stability and modulates translation. It participates in diverse biological processes, including RNA-protein interactions, splicing, translational control, and stress responses. Aberrant pseudouridylation is implicated in cancer, neurodegenerative disorders, and autoimmune diseases. Despite its biological importance, experimental identification of Ψ sites remains time-consuming and costly, limiting the feasibility of transcriptome-wide profiling. Computational approaches have therefore become essential complements to experimental techniques. However, state-of-the-art machine-learning and deep-learning predictors often suffer from limited generalizability due to small training datasets. To overcome these issues, we aim at constructing new long-sequence datasets and developing a novel Ψ site predictor. METHODS: New long-sequence datasets were constructed as benchmarks for RNA Ψ-site prediction. The Ψ modification sites in RMBase 3.0 were mapped to the reference genomes across three species of human, mouse, and yeast, and the RNA sequences with a length of 201 were generated by extending the upstream and downstream from the mapped, central sites. To eliminate sequence redundancy, the sequences were clustered using CD-HIT with a 70% sequence identity threshold. We developed Meta-PseU, a logistic regression-based meta-classifier that considered 118 machine learning and deep learning classifiers. The datasets and programs are freely accessible at https://github.com/kuratahiroyuki/MetaPseU. RESULTS: By optimizing model configuration, we proposed the Meta-PseU model stacking 32 machine learning and deep learning classifiers out of 118 classifiers. Meta-PseU substantially improved model generalizability, overcoming a key limitation of existing approaches. It greatly outperformed state-of-the-art predictors and achieved increasing accuracy with increasing sequence length. CONCLUSIONS: Long-sequence datasets were newly constructed as benchmarks for RNA Ψ-site prediction. Meta-PseU offers a new framework for robust Ψ-site identification by using long sequences.

Pseudouridine↗