Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “RNA Splicing”

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.

129 recordsLinked to original sources
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

Suppression of HIV-1 replication in CEM-A cell cultures by trans-splicing group I introns targeting PAS/PBS sequences and conditionally expressing ΔN-Bax.

Anti-HIV group I introns containing antisense guide sequences directed against the HIV-1 primer activation signal and primer-binding site (PAS/PBS) were designed and evaluated. Because PAS/PBS sequences are present in the viral RNA species examined, these RNAs can serve as trans-splicing substrates. The introns were active against both artificial target RNAs and viral RNA generated during infection. Cleavage and degradation of targeted viral RNA may have contributed to suppression, whereas inclusion of a 3' exon encoding the proapoptotic protein ΔN-Bax was associated with increased programmed cell death and may have augmented suppression of viral replication. In cultured CEM-A cells, transgene expression of these introns markedly suppressed HIV-1 replication, with p24 levels falling below the assay detection limit in selected clones. RESULTS: RT-PCR and sequence analysis detected splice products containing the expected PAS/PBS junctions. In the dual-luciferase assay, intron expression reduced normalized Gaussia luciferase signal by approximately 70% relative to the negative control. Qualitative Annexin V imaging and caspase-3 assays were consistent with infection-dependent apoptosis after ΔN-Bax splice-product formation. Transient expression of each intron in HEK293T cells followed by infection with VSV-G-pseudotyped HIV-1NL4-3 at an MOI of 2 reduced p24 levels by approximately 50% at 4 days post-infection. Construct 128L produced the strongest RT-PCR band under the tested conditions and was selected for subsequent experiments. A canonical splice product and a low-abundance noncanonical splice product were detected; both involved the intended HIV-derived target RNA, although transcriptome-wide off-target splicing was not assessed. Heterogeneous transformed HEK293T populations showed an approximately 2-log10 reduction in p24. In selected clonal HEK293T and CEM-A lines, p24 was below the assay detection limit at the measured endpoints, including up to 90 days after infection in some CEM-A clones. CONCLUSIONS: PAS/PBS-targeting group I introns suppressed HIV-1-associated p24 production in the tested cell-culture models. Linking the introns to a ΔN-Bax 3' exon was associated with infection-dependent apoptosis and may further limit viral replication and spread. The use of highly conserved, functionally constrained target sequences may reduce the likelihood of escape, but viral evolution and transcriptome-wide off-target effects were not assessed. This conditional death-upon-infection strategy warrants further evaluation in primary-cell and in vivo models.

Humans↗
PubMed2026

Novel splice site variants in GBA1 are associated with Gaucher disease and genotype-phenotype correlations.

BACKGROUND: Variants in GBA1 are associated with neurodegenerative disease. This study aimed to explore pathogenic GBA1 variants. METHODS: Four patients with progressive myoclonic epilepsy (PME) and extremely low β-glucosidase levels were recruited. Whole-exome sequencing and long-range PCR were performed to identify GBA1 variants. Bioinformatic analyses were used to predict the impact of the identified variants. A literature review was performed to explore the genotype-phenotype correlations. GBA1 expression data across different brain regions and developmental stages were analyzed using the BrainSpan database. RT-PCR was performed to verify the splicing effects. RESULTS: Compound heterozygous GBA1 variants were identified in four patients. Five distinct variants were detected, including two novel splice site variants (c.308-2A>G and c.762-2A>C) and three previously reported variants. All identified variants were rare or absent in gnomAD. Splice site variants c.308-2A>G and c.762-2A>C were predicted to cause aberrant splicing. Minigene-based splicing assays coupled with RT-PCR and Sanger sequencing confirmed that both variants cause complete exon skipping (exon 4 and exon 7, respectively). All patients presented with PME onset in childhood/adolescence, intellectual regression, low β-glucosidase, and diffuse brain atrophy and were subsequently diagnosed with Gaucher disease type 3. GBA1 expression in the brain showed two distinct peaks: one in infancy and another after five years of age. The onset age of PME aligned with the second GBA1 expression peak (after five years of age). CONCLUSION: This study identified compound heterozygous GBA1 variants, including two novel candidate pathogenic splice site variants, in Gaucher disease type 3 patients, expanding the known mutational spectrum.

Humans↗
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↗
PubMed2026

Genetic regulation of CPEB3-mediated alternative polyadenylation associated with survival of patients with hepatocellular carcinoma.

BACKGROUND: Alternative polyadenylation (APA) is a key post-transcriptional mechanism that regulates gene expression by modulating 3'UTR length, its dysregulation has been implicated in carcinogenesis. How genetic variants influence APA to affect hepatocellular carcinoma (HCC) prognosis remains unclear. METHODS: Prognosis-APA quantitative trait loci (apaQTL) were performed using genotype and APA profiling from TCGA data. A two-stage survival analysis in 848 Chinese and 369 TCGA LIHC patients and functional validation were used to identify prognostic apaQTL in HCC progression. RESULTS: A total of 2,025 and 817 significant APA events were identified in Chinese and TCGA cohort, respectively. Besides, 859 events were associated with poor prognosis in HCC and enriched in RNA splicing / metabolism pathways. We detected 32,034 significant apaQTLs, predominantly enriched in 3'UTRs and RBP-binding regions. CPEB3 was prioritized as a key APA regulator RBP; its low expression correlated with poor patient survival and promoted proliferation, migration, and invasion in HCC cells. Notably, a functional apaQTL variant rs2037547, located in GSK3B and mediated by CPEB3, demonstrated a poor survival of HCC patients in both cohort (pooled HR=1.29, p=0.016). Mechanistically, rs2037547 promoted aberrant APA at proximal poly(A) sites of GSK3B through CPEB3, leading to increased expression of short 3'UTR isoform. This regulatory alteration enhanced HCC cell proliferation, invasion, and migration, and contributed to HCC progression. CONCLUSION: These findings elucidated the distinct role of apaQTL-mediated APA dysregulation in HCC prognosis, providing insights for prognostic stratification and potential targets for personalized therapy in HCC.

RNA-binding proteins↗
PubMed2026

Frequent mutations in the BIRC3 gene promote metastatic potential of nasopharyngeal carcinoma cells through the TRAF2-NF-κB pathway.

Nasopharyngeal carcinoma (NPC) is a head and neck cancer characterized by highly locoregionally invasive behavior attributable to the latent infection with Epstein-Barr virus (EBV) and genomic instability. It is well established that EBV-encoded oncogenic molecules actively contribute to the malignant behavior of NPC cells. However, the mechanism by which aberrant genomic alterations enable NPC cells to become aggressive remains largely unknown. In the present study, whole-exome sequencing (WES) revealed that the gene encoding the baculoviral IAP repeat-containing 3 (BIRC3) protein was frequently mutated in circulating tumor cells (CTCs) but not in paired primary tumor cells from patients with metastatic NPC. A minigene assay indicated that the c.637 A > G mutation disrupted normal mRNA splicing, resulting in the partial deletion of Exons 2 and 3 and altered stability of BIRC3 mRNA. In vitro experiments demonstrated that ectopic expression of the BIRC3c.637A>G mutant enhanced NPC cell invasive properties, including proliferation, resistance to apoptosis, migration, and invasion. Furthermore, overexpression of wild-type BIRC3 promoted invasive characteristics in NPC cells through the TRAF2-NF-κB signaling axis. In summary, BIRC3 acts as a regulator of the malignant features of NPC cells. Frequent BIRC3 mutations in CTCs, such as the c.637 A > G mutation, further enhance the metastatic potential of disseminated NPC cells by inducing aberrant alternative splicing. These findings suggest the therapeutic feasibility of targeting the BIRC3/TRAF2/NF-κB axis in the treatment of NPC.

Humans↗
PubMed2026

[Analysis of clinical phenotypes and pathogenicity of a c.4476+5G>T variant of SCN1A gene in a Chinese pedigree affected with Genetic epilepsy with febrile seizures plus].

OBJECTIVE: To explore the pathogenicity and characteristics of a heterozygous splicing variant of SCN1A gene in a Chinese pedigree affected with Genetic epilepsy with febrile seizures plus (GEFS+). METHODS: A retrospective analysis was carried out on the clinical data and results of genetic testing of a GEFS+ pedigree consisting of 5 members who had visited the First Affiliated Hospital of Zhengzhou University on July 1, 2024. Pathogenicity of the splicing variant of the SCN1A gene was validated with a minigene splicing assay. This study was approved by the Medical Ethics Committee of the the First Affiliated Hospital of Zhengzhou University (Ethics No.: KS-2018-KY-36). RESULTS: The proband, a 24-year-old female, presented with FS in conjunct with focal seizures, and both of her younger brothers had Dravet syndrome. All of the three patients had carried a c.4476+5G>T variant of the SCN1A gene, which was unreported previously. Minigene experiment verified that the variant could cause loss of the first 7 bps of exon 24 and 138 bps from exon 23 of the SCN1A gene, resulting in alteration p.V1447_1495delfs*6 and affecting splicing. Based on the guidelines from American College of Medical Genetics and Genomics (ACMG), the variant was predicted as likely pathogenic (PVS1+PM2_Supporting). CONCLUSION: The c.4476+5G>T variant at an intronic site of the SCN1A gene probably underlay the pathogenesis of GEFS+ in this pedigree.

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

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

Integrative multi-omics reveals a fibroblast-centered, ZFHX3-prioritized regulatory framework linking sick sinus syndrome and atrial fibrillation.

OBJECTIVE: To define shared genetic and multi-scale mechanisms underlying comorbidity between sick sinus syndrome (SSS) and atrial fibrillation (AF). METHODS: We integrated genome-wide association study (GWAS) summary statistics for SSS and AF with Genotype-Tissue Expression (GTEx) expression and splicing quantitative trait loci (eQTL/sQTL), atrial single-cell and spatial transcriptomics, and epigenomics. We identified trait-relevant tissues and pathways, prioritized shared cell types, quantified genome-wide and local genetic sharing, detected joint loci by cross-trait meta-analysis, and linked loci to regulatory programs via colocalization and cell-prioritized co-expression networks. RESULTS: Both traits showed strongest enrichment in cardiac tissue, especially Heart Atrial Appendage. Fibroblasts from the left atrial appendage were consistently prioritized as the key shared cell population. SSS and AF displayed significant positive genome-wide genetic correlation, with multiple locally shared regions, including six major loci. Cross-trait meta-analysis identified eight joint-phenotype SNPs implicating four susceptibility genes. ZFHX3 was the leading tissue-cell-gene candidate, acting as a hub in fibroblast co-expression modules and colocalizing with cardiac regulatory signals. CONCLUSION: Shared liability for SSS and AF is highly tissue- and cell-specific, converging on regulatory networks in atrial appendage fibroblasts, with ZFHX3 serving as a central mechanistic and biomarker node.

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

Single-cell RNA sequencing provides further insights into the immunostimulatory action of freeze-dried Lactiplantibacillus plantarum on Penaeus vannamei shrimp.

Immunostimulation through dietary interventions opened new avenues in developing disease control and prevention tools for shrimp aquaculture. We have previously shown that feeding with freeze-dried Lactiplantibacillus plantarum (LAB) increased disease resistance of Penaeus vannamei against both Vibrio parahaemolyticus and white spot syndrome virus (WSSV) based on bulk RNA sequencing of shrimp gills. This tissue participates in ion transport and serves as a first line of defense against environmental stressors and pathogenic infections. However, characterization of their cell composition and functions remains limited. Here, we implemented a single-cell RNA sequencing approach to further gather insights into how feeding with freeze-dried LAB modulates host immunity which may not be evident with bulk RNA sequencing approach. A total of five clusters with unique transcriptional signatures were identified, corresponding to pillar cells, septal cells, and sessile hemocytes. Pseudo-bulk analyses at global- and cluster-levels showed differential expression of genes related to host immunity and metabolism. We further revealed how overall transcriptomic changes are not exclusively caused by gene expression changes but may also be driven by cell population dynamics. This study highlighted how single-cell RNA sequencing approach may shed light on the mechanisms of action of immunostimulants which may be masked in bulk transcriptome analyses.

Animals↗