Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “RNA sequencing analysis”

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.

At least 55 records · Page 3Linked to original sources
PubMed2026

The H3K27me3 reader GmLHP1 impairs Phytophthora sojae resistance by repressing ethylene precursor accumulation in soybean.

Phytophthora root rot, caused by Phytophthora sojae, is a devastating soilborne disease of soybean (Glycine max). However, the epigenetic regulation of soybean responses to P. sojae remains incompletely understood. Here, using genetic, molecular and biochemical approaches, we characterized the functions of LIKE HETEROCHROMATIN PROTEIN 1 (GmLHP1) and its downstream regulatory network. We demonstrated that GmLHP1, as a reader of H3K27me3, negatively regulates soybean resistance to P. sojae. GmLHP1 binds to H3K27me3 peptides in vitro and colocalizes with H3K27me3 marks genome-wide in vivo. The integrated chromatin immunoprecipitation sequencing and RNA sequencing analysis identified the ethylene biosynthesis pathway gene 1-AMINO-CYCLOPROPANE-1-CARBOXYLATE SYNTHASE 18 (GmACS18) as being enriched for H3K27me3 and bound by GmLHP1, leading to its transcriptional downregulation. Notably, GmLHP1 associates with the GmACS18 promoter by directly binding to AATTAA motifs and recognizing H3K27me3 marks. Moreover, GmACS18 enhances defense against P. sojae by accumulating the ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC). Further analysis unveiled that recognition of H3K27me3 by GmLHP1 is essential for regulating soybean resistance to P. sojae through repressing GmACS18 transcription and decreasing ACC accumulation. Our findings reveal a novel epigenetic regulatory mechanism in which the H3K27me3 reader GmLHP1 blocks soybean resistance to P. sojae by repressing ethylene precursor ACC accumulation.

ACC accumulation↗
PubMed2025

The Tartary Buckwheat FtMYB46-FtNRAMP3 Module Enhances Plant Lead and Cadmium Tolerance.

The presence of toxic heavy metals lead (Pb) and cadmium (Cd) in polluted soil damage crop production and consequently harms human and livestock health. Tartary buckwheat (Fagopyrum tataricum) is a potential model plant for heavy metal phytoremediation because of its valuable characteristics of high heavy metal tolerance and abundant biomass production. Here, we report that the Tartary buckwheat FtMYB46-FtNRAMP3 module enhances plant Pb and Cd tolerance. RNA sequencing analysis showed that Pb treatment specifically induced expression of FtNRAMP3, a member of the NRAMP (Natural Resistance-Associated Macrophage Protein) transporter gene family. Further cytological and biochemical analysis revealed that FtNRAMP3 was localised to the plasma membrane and significantly contributed to increased tolerance to Pb and Cd in yeast cells. Consistently, transgenic overexpression of FtNRAMP3 in Arabidopsis significantly increased plant tolerance to Pb and Cd applications, reducing Pb concentration but increasing Cd concentration in the overexpression transgenic plants. Subsequent yeast one-hybrid and electrophoretic mobility shift assays showed that the transcription factor FtMYB46 directly binds to the FtNRAMP3 promoter. Further, FtMYB46 promoted FtNRAMP3 expression and increased plant Pb and Cd tolerance. Overall, this study demonstrates the important role of the FtMYB46-FtNRAMP3 module and its potential value in the phytoremediation of Pb and Cd stress.

Cadmium↗
PubMed2025

NAM and CUC3 boundary genes maintain shoot apical meristem viability and suppress the development of axillary shoot in rice seedlings.

Cell division and differentiation within the shoot apical meristem (SAM) are essential for the morphogenesis of aboveground plant organs. This study reveals that the boundary genes OsNAM and OsCUC3 collaboratively maintain SAM activity. Loss of function in both OsNAM and OsCUC3 during the fourth leaf stage reduced SAM size, with the osnam oscuc3 mutant exhibiting abnormal leaf number and morphology. Furthermore, OsNAM and OsCUC3 inhibited the growth of axillary shoots. In the osnam oscuc3 mutant, the number of new leaves decreased, while buds in the coleoptile and the axil of the first leaf developed into tillers. Since OsNAM and OsCUC3 are involved in regulating both SAM activity and the growth of lateral shoots, we examined their expression patterns at the base of the main shoot. β-Glucuronidase (GUS) reporter activity and GFP reporter lines demonstrated that OsNAM and OsCUC3 have distinct expression patterns. Specifically, OsNAM was expressed throughout the SAM, whereas OsCUC3 was expressed only at the base of the SAM, with its expression gradually decreasing as seedlings develop. RNA sequencing analysis showed that the expression of genes related to leaf epidermal cell development, cell wall components, and hormonal signal transduction was altered in response to the loss of function of OsNAM and OsCUC3. Therefore, the boundary genes OsNAM and OsCUC3 not only inhibit the growth of axillary shoots but also regulate the development of aboveground organs, including leaf morphology and number, by maintaining the SAM activity in the main shoot.

Meristem↗
PubMed2026

Defective EV-mediated transport of SHH alters neural fate specification in EPM1 epilepsy.

The extracellular milieu, including extracellular vesicles (EVs), plays a pivotal role in brain development. In this study, we sought to elucidate the pathogenesis of progressive myoclonus epilepsy type 1 (EPM1), a disease caused by mutations in the CSTB gene, using cerebral organoids (COs) derived from patient cells. The results demonstrate that EPM1 COs display increased electrophysiological activity and disrupted excitatory/inhibitory (E/I) balance. Single-cell RNA sequencing analysis of ventral EPM1-COs revealed an abnormal specification of progenitor fate, with a shift toward dorsal neuron identities. We demonstrated that this misspecification is driven by a functional alteration of the ventral signaling niche, resulting from impaired EV dynamics and altered protein cargo. Mechanistically, we identified Sonic Hedgehog (SHH) as a direct physical interactor of CSTB and demonstrated that CSTB deficiency leads to reduced SHH content and secretion. Our findings establish CSTB as a safeguard of ventral patterning and identify the CSTB-SHH-EV axis as a potential therapeutic target for mitigating the E/I imbalance associated with EPM1.

Hedgehog Proteins↗
PubMed2025

The skeletal muscle of aged male mice exhibits sustained growth regulatory transcriptional profile following glucocorticoid exposure compared with young males.

Excess glucocorticoids induce skeletal muscle myopathy by changing gene expression. Advanced age augments glucocorticoid-mediated muscle phenotypes, yet the transcriptional responses underlying those augmented phenotypes are unclear. The purpose of this study was to define the glucocorticoid-responsive transcriptome in young and aged muscle following both acute and more prolonged glucocorticoid treatment. Young (4-mo-old) or aged (24-mo-old) male mice were administered either an acute injection of dexamethasone (DEX) or vehicle or daily DEX or vehicle injections for 7 days. Muscles were harvested 6.5 h after the final or only injection. The tibialis anterior (TA) was selected for RNA sequencing analysis as DEX treatment lowered TA mass specifically in aged males. In silico analyses identified enriched pathways and transcription factors predicted to regulate DEX-sensitive genes. Acute DEX altered similar numbers of genes in young (950) versus aged males (913), although aged males had greater magnitudes of fold change. After 7 days of DEX treatment, aged muscle exhibited more DEGs compared with acute exposure (1,196 vs. 913), whereas young muscle exhibited fewer DEGs than after acute exposure (599 vs. 950). In aged males, glucocorticoid-sensitive genes were consistently enriched for growth regulatory processes across both time points, a pattern that was not evident in young males. Despite those age-associated transcriptional differences, the transcription factors predicted to regulate the glucocorticoid-sensitive genes were similar in young and aged males. These data expand our understanding into how aging modifies the transcriptional response to excess glucocorticoids in skeletal muscle.NEW & NOTEWORTHY Glucocorticoids promote mass loss in certain muscles with advanced age but not at younger ages. In a muscle whose mass is lost in response to elevated glucocorticoids only in advanced age in males, we show that glucocorticoids initiate a unique and exaggerated transcriptional profile after both acute exposure to the hormone and after prolonged treatment that is consistent with muscle atrophy. These findings expand our understanding of the effect primary aging has on glucocorticoid-induced atrophy in males.

Animals↗
PubMed2026

PRDM16 Regulates Prostate Cancer Cell Dormancy and Prevents Bone Metastatic Outgrowth.

UNLABELLED: Understanding dormancy in prostate cancer is challenging because of model availability. In this study, using murine and human prostate cancer cell lines, we generated a stress-induced model of dormancy in vitro and demonstrated that the phenotype could be sustained upon intrailiac artery delivery into the bone marrow microenvironment. RNA sequencing analysis revealed that the transcription factor positive regulatory domain-containing 16 (PRDM16) was commonly upregulated in dormant prostate cancer cells compared with controls. Furthermore, bone marrow-disseminated prostate cancer cells from primary orthotopic tumors were largely positive for PRDM16. Genetic ablation and forced ectopic expression supported a role for PRDM16 in maintaining prostate cancer dormancy in vitro and in vivo. Clinically, PRDM16 negatively correlated with disease recurrence and with the E2F cell-cycle program in disseminated tumor cells derived from the bone marrow of patients with prostate cancer. Gene enrichment and characterization studies implicated PRDM16 as a regulator of metabolic and cell-cycle pathways. Chromatin immunoprecipitation-qPCR further revealed that PRDM16 binds upstream of the promoter of RB1, a potent repressor of E2F activity. Overall, this study developed a straightforward method for inducing cancer cell dormancy and applied this approach to find that PRDM16 governs an intrinsic dormancy program in prostate cancer. SIGNIFICANCE: PRDM16 initiates a dormancy program in prostate cancer cells that is sustained in the bone marrow microenvironment, highlighting PRDM16 as a potential biomarker for relapse and target for eliminating dormant cancer cells.

Male↗
PubMed2026

BET inhibition blunts antibody production and macrophage-mediated fibrosis to restore lung function in murine cGVHD.

Chronic graft-versus-host disease (cGVHD) significantly contributes to late mortality after allogeneic stem cell transplantation, with bronchiolitis obliterans syndrome (BOS) being a particularly lethal and treatment-resistant complication despite available therapies. Bromodomain and extraterminal (BET) proteins are epigenetic readers driving inflammatory transcriptional programs across multiple cell types. We hypothesized that BET inhibition would suppress inflammatory T and B cells and decrease macrophage polarization to a profibrotic phenotype, alleviating disease. In an established BOS cGVHD model, BET inhibition reduced germinal center (GC) formation and responses through a reduction of the CXCL13:CXCR5 axis and inflammatory T follicular helper/GC B cells in the spleen, along with a reduction in plasma cell infiltration within the lung. Mice with cGVHD had elevated pathogenic immunoglobulin G1 (IgG1) and IgM levels, both in circulation and deposited on lung tissue, which were attenuated under BET inhibition. Single-cell RNA-sequencing analysis revealed distinct cell states in the BOS lung vs control. In cGVHD mice, gene set enrichment analysis revealed the upregulation of profibrotic Arginase1 and Tgfb1 expression in alveolar macrophages (AM) and interstitial macrophages (IM), which was significantly reduced with BET inhibition. Furthermore, BET inhibition targeted lung-infiltrating M2 macrophages through the selective depletion of CD206+FcγR+ IM and AM, ultimately resulting in reduced collagen deposition and improved lung function. Our findings reveal a previously unrecognized mechanistic axis of BET regulation during cGVHD fibrosis and highlight BET inhibition as a promising therapeutic strategy.

Animals↗
PubMed2026

A token-pruning framework enables efficient representation of the human genome for RNA modification analysis.

MOTIVATION: Modelling long genomic sequences remains challenging due to extreme sequence length, high redundancy, and the need for biological interpretability. Although Transformer-based architectures have achieved strong performance across genomic tasks, their high computational cost and reliance on fixed tokenization strategies limit their scalability and ability to focus on biologically informative regions. RESULTS: We propose ATSFormer, a token-pruning Transformer framework for efficient and biologically informed genomic sequence modelling. ATSFormer incorporates an attention-guided and parameter-free Adaptive Token Sampling (ATS) module into Transformer layers. Guided by attention-derived importance scores, ATS dynamically retains informative tokens while probabilistically discarding redundant ones, thereby reducing sequence length, FLOPs, and memory usage without introducing additional learnable parameters or extra training procedures. Importantly, the retained tokens correspond to key contributors to model predictions, enabling ATSFormer to highlight biologically meaningful sites and sequence motifs. We evaluated ATSFormer on four benchmark RNA modification datasets derived from RMVar 2.0, covering A-to-I, m1A, m5C, and m7G. Experimental results show that ATSFormer consistently outperforms existing state-of-the-art methods while achieving substantial computational savings. Furthermore, structural analysis using AlphaFold3 supports the biological relevance of the motifs identified by ATSFormer. AVAILABILITY AND IMPLEMENTATION: The source data and code are freely available at GitHub (https://github.com/1gao2/ATSFormer) and Zenodo (https://doi.org/10.5281/zenodo.21813541).

Humans↗
PubMed2026

PETScan: score-based genome-wide association analysis of RNA-Seq and ATAC-Seq data.

MOTIVATION: High-dimensional sequencing data, such as RNA-Seq for gene expression and ATAC-Seq for chromatin accessibility, are widely used in studying systems biology. Accessible chromatin allows transcription factors and regulatory elements to bind to DNA, thereby regulating transcription through the activation or repression of target genes. The association analysis of RNA-Seq and ATAC-Seq data provides insights into gene regulatory mechanisms. Most existing analytic tools exclusively focus on cis-associations, despite regulatory elements being able to physically interact with distant target genes. Furthermore, conventional approaches often utilize Pearson or Spearman correlations, which ignore the count-based nature of RNA-Seq data. RESULTS: To address these limitations, we introduce PETScan, a computationally efficient genome-wide PEak-Transcript Score-based association analysis, utilizing negative binomial models to better accommodate RNA-Seq data. We leverage score tests and matrix calculations for improved computational efficiency, and combine an empirical permutation method with genomic control to ensure valid p-value calculations in studies with limited sample sizes. In real-world datasets, PETScan achieved three orders of magnitude faster than Wald tests, while identifying similar significant gene-peak pairs. AVAILABILITY: The PETScan R package is available on GitHub at https://github.com/yajing-hao/PETScan.

Chromatin Immunoprecipitation Sequencing↗
PubMed2026

ModiCal: A Targeted Calibration Workflow for Site-Specific m5C Validation by Nanopore Direct RNA Sequencing.

Accurate identification of RNA 5-methylcytidine (m5C) at the single-nucleotide resolution remains a central challenge in nanopore direct RNA sequencing (DRS). Current global scanning and modification-aware basecalling methods enable transcriptome-wide profiling but often yield high false-positive rates and lack site-specific accuracy. To address this, we repurposed ModiDeC, originally a de novo multimodification classifier, into a targeted, high-precision validation tool for RNA modification sites with prior biochemical knowledge. This was implemented through a three-step calibration workflow that alternates between biochemical and computational modules using the well-characterized m5C2278 site in 25S rRNA as a starting point. Baseline training uses short synthetic RNAs carrying either a methylated or unmodified C2278 as ground truth, followed by IVT-derived calibration and validation in methyltransferase knockout yeast. The baseline model accurately detected the bona fide m5C2278 site but initially produced off-target predictions. Iterative retraining with unmodified IVT signals progressively reduced and ultimately eliminated false positives while maintaining a strong signal at the bona fide site. The final model retained enzyme-dependent detection in wild-type versus knockout yeast and, when explicitly targeted, was also able to detect the second rRNA site, C2870, which remained invisible in the initial analysis. Application to native human prerRNA processing intermediates further resolved two distinct m5C deposition regimes on 28S rRNA, while generalization to dengue virus genomic RNA confirmed that the same calibration logic transfers across diverse RNA contexts. Together, this study establishes a reproducible and transferable framework that integrates biochemical validation with iterative neural network refinement, providing a route toward reliable site-specific m5C confirmation by nanopore direct RNA sequencing.

RNA Methylation↗
PubMed1977

The use of nuclease P1 in sequence analysis of end group labeled RNA.

A method is described for the direct sequence analysis of 20-25 nucleotides from the termini of 5'- or 3'-end-group [32P] labeled RNA. The method involves partial endonucleolytic digestion of the labeled RNA with nuclease P1 (from Penicillium citrinum) followed by separation of the partial digestion products by two-dimensional homochromatography, the nucleotide sequence being determined by mobility shift analysis. This procedure has been applied to the sequence analysis of the terminal regions of tRNAs and of high molecular weight RNA, such as messenger RNA or viral RNA. A further application involves its use in conjunction with snake venom phosphodiesterase to determine the sequence of 5'-end group labeled oligonucleotides, containing modified bases, derived from T1 or pancreatic RNase digestion of tRNA.

Base Sequence↗
PubMed2025

Transcriptome-Wide Analysis of the 5' Cap Status of RNA Using 5' Monophosphate-Dependent Exonuclease Digestion and RNA Sequencing.

Eukaryotic mRNAs carry an N7-methylguanosine (m7G) cap structure at their 5' extremity, which protects them from the degradation by 5'-3' exoribonucleases and plays a pivotal role in mRNA metabolism, promoting splicing, nuclear export, and translation. Decapping, the enzymatic process that removes this structure, is a key event during cytoplasmic mRNA 5'-3' decay, leading to the degradation of the transcript body by Xrn1. In this chapter, we describe a procedure to assess the cap status of RNA at the transcriptome level. It is based on a treatment of total RNA extracts with a 5' monophosphate-dependent exonuclease, which like Xrn1 specifically degrades decapped RNAs harboring 5' monophosphate extremities, but not RNAs with intact m7G cap. The digested RNAs are then analyzed by RNA sequencing.

Exoribonucleases↗
PubMed1979

Approaches to sequence analysis of 125I-labeled RNA.

A method is described for the initial steps of sequence analysis of RNase T1-and pancreatic RN-ase-resistant oligonucleotides of RNA containing cytidylate residues labeled in vitro with 125I. In many cases an oligonucleotide sequence can be deduced from a consideration of (i) its relative position in the two-dimensional fingerprint (with DEAE thin layer homochromatographic second dimension), (ii) its electrophoretic mobility on DEAE paper at pH 1.9, and (iii) identification of its products of further enzymatic digestion by comparison with a set of marker oligonucleotides. Additional methods including analysis of oligonucleotides following chemical blocking of uridylate residues with CMCT and analysis of products of incomplete enzymatic digestion are also discussed.

Base Sequence↗
PubMed2026

Characterization of METTL3/14-mediated m6A modification in human transcriptome using Nanopore direct RNA sequencing.

Post-transcriptional RNA modifications modulate diverse aspects of RNA metabolism. N6-methyladenosine (m6A), one of the most abundant internal RNA modifications, is deposited by the core methyltransferase complex, METTL3 and METTL14. Oxford Nanopore Technologies (ONT) platform permits direct, single RNA molecule sequencing while preserving native modifications. However, without rigorous benchmarking, the accuracy and reproducibility of modification detection remain uncertain. Here, we leveraged ONT to comprehensively profile bona fide m6A modifications in cellular RNAs at single-nucleotide resolution by integrating two direct RNA sequencing chemistries (RNA002 and RNA004) with the m6Anet and Dorado modification-detection models. We independently depleted METTL3 and METTL14 in human cells and rigorously validated modification calls through several assays and independent orthogonal methods (GLORI and miCLIP). We find that Dorado detected a higher number of m6A events and enabled simultaneous detection of other RNA modifications (5-methylcytosine, pseudouridine, and inosine). Pairing Dorado with an in vitro transcribed, unmodified control under stringent filtering, we provide compelling evidence supporting a global reduction in m6A sites and stoichiometry within coding sequences and across genes, particularly in highly modified genes and sites, and at consensus DRACH motifs. We report a differential and complex regulation of modified transcripts, accompanied by a global reduction in poly(A) tail length. Notably, METTL3 and METTL14 depletion produced distinct transcript-specific effects, supporting non-redundant roles within the m6A writer complex. Together, our study illustrates a notable advancement of ONT capabilities and establishes a robust transcriptome-wide framework for RNA modification detection, thereby laying the groundwork for exploring the contribution of METTL3/METTL14 to cellular functions and disease.

Humans↗
PubMed2024

Identification of RBP binding sites using RNA deaminases.

RNA-binding proteins (RBPs) are critical regulators of gene expression and RNA processing. Identification of their binding sites has important implications for their physiological and disease-related functions. Crosslinking and immunoprecipitation, followed by sequencing (CLIP-seq) and its derivatives, are the most commonly used methods to identify RBP binding sites, but are laborious and require a large amount of starting material. Recent advancements harnessing RNA deaminases in fusion to any RBP of interest, allow for the profiling of RBP binding sites from low-input samples in simpler procedures. Among these efforts, we developed STAMP (Surveying Targets by APOBEC-Mediated Profiling), which efficiently detects RBP-RNA interactions. This chapter describes the detailed protocol for the STAMP method, including plasmid construction, delivery and sorting, library preparation and bioinformatic data analysis.

RNA-Binding Proteins↗
PubMed2019

Genome-wide analysis of polymerase III-transcribed Alu elements suggests cell-type-specific enhancer function.

Alu elements are one of the most successful families of transposons in the human genome. A portion of Alu elements is transcribed by RNA Pol III, whereas the remaining ones are part of Pol II transcripts. Because Alu elements are highly repetitive, it has been difficult to identify the Pol III-transcribed elements and quantify their expression levels. In this study, we generated high-resolution, long-genomic-span RAMPAGE data in 155 biosamples all with matching RNA-seq data and built an atlas of 17,249 Pol III-transcribed Alu elements. We further performed an integrative analysis on the ChIP-seq data of 10 histone marks and hundreds of transcription factors, whole-genome bisulfite sequencing data, ChIA-PET data, and functional data in several biosamples, and our results revealed that although the human-specific Alu elements are transcriptionally repressed, the older, expressed Alu elements may be exapted by the human host to function as cell-type-specific enhancers for their nearby protein-coding genes.

Alu Elements↗
PubMed2025

Analysis of Leishbuviridae from Trypanosomatids.

Over the last decade, considerable progress has been made in unraveling RNA virus diversity. This has contributed to our understanding of the evolution of these viruses, which include emerging zoonotic human pathogens. Current success has been greatly facilitated by the development of next-generation sequencing platforms instrumental for meta-transcriptomic studies. However, due to the rapid evolution of RNA viruses, there are numerous "blind spots" waiting to be explored; one of those is the RNA virome of unicellular eukaryotes. Here, we present the pipeline, which has been successfully used to characterize various types of RNA viruses, including Leishbuviridae (Bunyaviricetes, Hareavirales) in the parasitic flagellates of the family Trypanosomatidae. The pipeline relies on axenic in vitro cell culture and double-stranded RNA enrichment, followed by direct RNA-sequencing. A detailed procedure description starting from the initial total RNA preparation to the final assembly of the viral segments is provided.

High-Throughput Nucleotide Sequencing↗