SEARCH · Search PubMed

Results for “single-cell RNA sequencing”

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.

374 recordsLinked to original sources

PubMed2026

Single-cell transcriptome revealed the aberrant keratinocytes activation in antigen presentation in atopic dermatitis.

BACKGROUND: Atopic dermatitis (AD), a common chronic inflammatory skin disease, has been extensively studied using single-cell genomics. However, keratinocytes, as key effector cells in AD, have underlying mechanisms remain incompletely understood and require further investigation. METHODS: We integrated single-cell transcriptomic data from skin tissues of healthy controls, chronic active AD patients, spontaneously healed AD (SHAD) patients, and an ovalbumin-induced AD mouse model. The study particularly emphasized the gene expression and cellular dynamics of keratinocytes across the different groups, as well as their interactions with immune cells. RESULTS: Compared to healthy controls, we observed significant changes in the keratinocyte transcriptome, cellular state, and keratinocyte-immune cell ligand-receptor interactions in AD skin, particularly the marked activation of genes involved in antigen processing and presentation. Interestingly, such gene activation was not observed in keratinocytes from the ovalbumin-induced AD mouse model, despite its phenotype closely resembling human AD. Furthermore, in SHAD, we identified a recovery of both the ligand-receptor interaction patterns and antigen processing and presentation genes, accompanied by a notable shift in the transcriptome. This involved a significant downregulation of genes related to cytoplasmic transcription and oxidative phosphorylation. Notably, this pattern was not observed in the self-healing mouse model following the removal of ovalbumin stimulation. CONCLUSION: Our results suggest that the persistent activation of antigen processing and presentation pathways in keratinocytes may be a key driver of chronic inflammation in AD. Therefore, redirecting anti-allergic therapeutic strategies from solely targeting immune cells to targeting of keratinocyte-mediated antigen presentation may offer a more effective approach. Furthermore, we raise concerns about the use of ovalbumin-induced mouse models to recapitulate human chronic AD, as the underlying mechanisms may differ significantly.

Dermatitis, Atopic↗

PubMed2026

ScRNA-seq analysis reveals the effects of nitrite stress on the endocrine system of the eyestalk in Litopenaeus vannamei.

Nitrite is a harmful substance generated in Litopenaeus vannamei farming systems, largely originating from the inadequate breakdown of surplus feed and shrimp feces. Its accumulation in the water can affect the growth and physiological functions of shrimp, damage the immune system, and even cause mass mortality, thus becoming a key environmental factor restricting the green development of the industry. Under nitrite stress, the eyestalk, as an important neuroendocrine regulatory center in crustaceans, participates in the stress adaptation of the organism and exerts a protective effect by regulating energy metabolism and immune function. However, the molecular regulatory mechanism of the eyestalk in response to nitrite stress remains unclear. In this study, single-cell RNA sequencing (scRNA-seq) technology was used to analyze the heterogeneity of eyestalk cells in L. vannamei under nitrite stress. A total of 18, 394 high-quality cells were obtained, and six major cell subpopulations, including Neurosecretory cell, Motor neuron, Sensory neuron, Interneuron, Neurogliocyte, and Support cell, were identified. Differential expression analysis identified 839 differentially expressed genes, and different cell types showed distinct specific responses to nitrite stress. Functional enrichment analysis indicated that pathways such as glycolysis, oxidative phosphorylation, ribosome function, and endoplasmic reticulum protein processing were significantly activated, while signal transduction and DNA repair-related pathways were inhibited. Further analysis revealed that nitrite stress could induce mitochondrial function changes and trigger oxidative stress, thereby affecting the neuroendocrine system function of the eyestalk. This study provided insights into transcriptomic responses of the eyestalk to nitrite stress at the single-cell level, laying a theoretical foundation for the management of aquaculture environments.

Animals↗

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

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

Diagnostic value of plasma cell-free DNA metagenomic next-generation sequencing in patients with suspected infections and exploration of clinical scenarios-a retrospective study from a single center.

BACKGROUND: Plasma cell-free DNA metagenomic next-generation sequencing (mNGS) is a non-invasive comprehensive method for the etiological diagnosis of various infectious diseases. However, research on the early diagnosis and real-world clinical impact of plasma mNGS in patients with suspected infection are still limited. MATERIALS AND METHODS: This study retrospectively included 140 patients with suspected infections who underwent early plasma mNGS and conventional culture testing. Referring to the clinical diagnosis of infectious diseases, the diagnostic performance of plasma mNGS and culture tests was compared, and the application scenarios and clinical effects of plasma mNGS were evaluated. RESULTS: The positive rate of plasma mNGS was significantly higher than that of culture methods (55.71% vs 25.10%, p < 0.001) and blood cultures (55.71% vs 12.86%, p < 0.001). Regarding clinical diagnosis, the sensitivity of plasma mNGS was significantly higher than that of culture (58.27% vs 37.80%, p = 0.002). The combination of mNGS and culture achieved a higher detection sensitivity (69.29%), especially in patients with multi-site co-infections (73.68%) and blood infections (73.17%). Plasma mNGS demonstrated higher sensitivity in patients with procalcitonin (PCT) index > 5 ng/ml or human neutrophil lipocalin (HNL) index > 200 ng/ml. In terms of treatment, a total of 69 patients (54.33%) benefited from plasma mNGS. CONCLUSION: This study highlights the significant improvement in pathogen detection performance by combining conventional culture with plasma mNGS detection, especially in patients with multi-site co-infections and blood infections. Early use of plasma mNGS as an adjunct to culture can better guide clinicians to initiate appropriate anti-infective therapy.

Humans↗

PubMed2026

Peripheral pain threshold, glycaemic status, and LAMP3 genetic variation: A community-based analysis.

Diabetic polyneuropathy is a common complication of diabetes, yet substantial inter-individual variation in peripheral pain perception suggests underlying genetic influences. This population-based study investigated clinical, metabolic, and genetic determinants of pain threshold using intraepidermal electrical stimulation in 906 participants from the Iwaki Health Promotion Project 2017. Genome-wide association analysis identified 12 loci showing suggestive associations, among which a missense variant in LAMP3 (rs482912) was prioritized as a biologically plausible candidate. Phenotype-stratified analyses showed that individuals carrying the CT or CC genotypes had lower PINT indices than those with the TT genotype, indicating reduced pain thresholds. Notably, the CC genotype retained an association with lower pain threshold using intraepidermal electrical stimulation under conditions of metabolic stress, including impaired glucose tolerance, elevated HbA1c, and obesity, whereas this association was attenuated in the presence of hypertension. Single-cell RNA sequencing analysis of human skin revealed that LAMP3-positive mature dendritic cells, enriched in immunoregulatory molecules, exhibited transcriptional enrichment of inflammatory, antigen-presenting, and nociception-related pathways, including NF-κB, JAK-STAT, cytokine signaling, and neuroimmune sensitization cascades. Autopsy-based skin analysis further demonstrated genotype-associated differences in dermal LAMP3-positive cell infiltration and CD8-positive T-cell abundance, while CD4-positive T-cell abundance and intraepidermal nerve fiber density remained unchanged across genotypes. Taken together, these findings suggest a potential association between LAMP3 variation and individual differences in peripheral pain threshold and provide biological context supporting a role for neuroimmune interactions in early sensory modulation under metabolic stress. Given the suggestive genetic evidence and indirect mechanistic data, these observations should be interpreted as exploratory and hypothesis-generating.

Humans↗

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

ER proteostasis failure in HYOU1 deficiency alters B cells, neutrophils, and interferon signalling.

Hypoxia upregulated 1 (HYOU1) is a stress-inducible ER chaperone. We investigated 2 unrelated patients carrying biallelic HYOU1 variants and presenting with primary immunodeficiency. Patient 1, homozygous for p.Pro444His, displayed failure to thrive, hypoglycemia, B cell lymphopenia, and neutropenia. Patient 2, compound heterozygous for p.Arg262Gln and p.Pro757_Glu758insAla, exhibited recurrent infections, enteropathy, and hypogammaglobulinemia. In Patient 1, while HYOU1 transcription was preserved, the protein was severely reduced. Tunicamycin treatment of dermal fibroblasts showed a blunted unfolded protein response and defective induction of ER stress-responsive genes. Immunophenotyping showed near-absence of circulating B cells, and single-cell RNA sequencing of bone marrow identified an arrest at the pro-B cell stage. Neutrophils displayed hypogranulation and dysregulated IFN- and apoptosis-associated transcriptional signatures, unresponsive to G-CSF. HYOU1 deficiency hence results in ER stress-induced proteostasis failure that simultaneously impairs adaptive immunity through B cell developmental arrest and innate immunity through neutrophil dysfunction and IFN pathway imbalance. This work expands the spectrum of HYOU1 deficiency and further identifies ER proteostasis as a central determinant of immune homeostasis.

Journal Article↗

PubMed2026

Integrated multi-omics analyses identify an RAS-SLC11A2-associated molecular framework linking iron metabolism with PCOS-related cardiometabolic risk.

INTRODUCTION: PCOS is a common endocrine disorder with elevated cardiometabolic risk, yet the role of the renin-angiotensin system (RAS)-iron metabolism axis in this comorbidity remains unclear. We explored its underlying mechanisms and evaluated the therapeutic potential of gentiopicroside. METHODS: Integrated multi-omics analyses combining transcriptomics, single-cell RNA sequencing, Mendelian randomization, machine learning, molecular docking, and in vitro functional assays were performed to identify shared molecular pathways and therapeutic targets across PCOS, hypertension, NAFLD, and T2DM. RESULTS: SLC11A2 was consistently dysregulated in PCOS transcriptomic datasets, and associated with iron metabolism, inflammatory response and oxidative stress pathways. Genetic analyses validated RAS-related regulation in hypertension susceptibility and revealed shared genetic architecture between PCOS and cardiometabolic traits. Network and single-cell analyses characterized SLC11A2-associated molecular patterns in disease-relevant cell types; machine learning identified disease-classifying molecular signatures. Gentiopicroside alleviated inflammatory and oxidative stress phenotypes, including reduced IL-6 expression and reactive oxygen species accumulation. CONCLUSION: This study defines an RAS-SLC11A2 molecular framework linking iron metabolism dysregulation to PCOS-related cardiometabolic risk, elucidating the mechanisms connecting ovarian dysfunction, inflammation, oxidative stress and hypertension, and supports gentiopicroside as a promising therapeutic candidate.

Humans↗

PubMed2026

Integrated bioinformatics analysis reveals cross-talking hub genes and therapeutic agents between sepsis and acute myocardial infarction.

BACKGROUND: Sepsis and acute myocardial infarction (AMI) are two significant diseases that may share overlapping etiological mechanisms. This study aims to systematically identify core genes common to both conditions and to explore their potential as therapeutic targets and drug candidates through an integrative analysis of clinical data and bioinformatics. METHODS: The AMI dataset was obtained from the GEO database, and RNA sequencing data were collected from blood samples of patients with sepsis at our hospital. Common genes were identified using differential expression gene analysis (DEG) and weighted gene co-expression network analysis (WGCNA). Functional enrichment analyses, including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, were performed. A protein-protein interaction (PPI) network was constructed, and hub genes were identified using the MCC/Degree algorithm. Diagnostic value was assessed via receiver operating characteristic curve analysis. Immune infiltration patterns, single-cell sequencing data, and molecular docking simulations were employed to evaluate immune relevance and identify potential therapeutic compounds. RESULTS: A total of 417 genes were identified between sepsis and AMI, with enrichment analysis revealing significant involvement in inflammatory responses. Three hub genes-JAK2, MYD88, and TIMP1-were selected for further investigation. ROC curves confirmed their strong diagnostic performance for both diseases. Immune infiltration analysis showed that these core genes were significantly correlated with the infiltration levels of various immune cell types. Molecular docking indicated that quercetin exhibited stable binding affinity with the proteins encoded by these genes. qPCR validation further confirmed the upregulation of these three genes, supporting the anti-inflammatory effects of quercetin as a potential targeted therapy. CONCLUSION: JAK2, MYD88, and TIMP1 were identified as shared core genes in sepsis and AMI. These genes not only serve as potential diagnostic biomarkers but also offer novel targets for developing common therapeutic strategies for both conditions. Furthermore, quercetin emerges as a promising candidate for targeted treatment.

Humans↗

PubMed2026

Molecular characterization of 16 MAPK genes in silver carp (Hypophthalmichthys molitrix) and the differences of their mRNA expression between Qiandao Lake and Taihu Lake.

Mitogen-activated protein kinase (MAPK), a serine-threonine protein kinase, is involved in a variety of stress-induced responses and also plays an important regulatory role in cell metabolism. In the study the open reading frames (ORFs) of 16 MAPK genes in silver carp (Hypophthalmichthys molitrix) were obtained and verified, with the evaluations of their taxonomy, structures, conserved motifs, and evolutionary linkages. And the expression patterns of these genes in the silver carp from Qiandao Lake and Taihu Lake were explored for better understanding the response of MAPK genes to different water environment. MAPK genes of silver carp were divided into three subfamilies, including extracellular signal-regulated kinase (ERK) subfamily, p38 subfamily and C-Jun N-terminal kinase (JNK) subfamily. All these genes possessed similar structures and conserved motifs of MAPK family. Realtime qPCR revealed that the expression patterns of 10 MAPK genes (ScMAPK1, ScMAPK3, ScMAPK4, ScMAPK7, ScMAPK15, ScMAPK8a, ScMAPK8b, ScMAPK9, ScMAPK10 and ScMAPK11) in head kidney, spleen and gill of silver carp in Taihu Lake and Qiandao Lake were different. These findings provide a basis for further research on the function of MAPK in silver carp.

Animals↗

PubMed2026

A novel peptide encoded by circTLL1 drives osimertinib resistance in lung cancer by modulating the NT5C2/Ras/PI3K axis.

BACKGROUND: Acquired resistance to osimertinib, a third-generation EGFR tyrosine kinase inhibitor, remains a major clinical challenge in the treatment of non-small cell lung cancer (NSCLC). Although circular RNAs (circRNAs) have been increasingly implicated in drug resistance, most studies have focused on their canonical role as microRNA sponges, while their capacity to encode functional micropeptides remains largely unexplored. This study aimed to identify novel circRNAs involved in osimertinib resistance and to characterize their regulatory functions at the protein level. METHODS: Osimertinib-resistant (OR) NSCLC cell lines were established and validated. High-throughput RNA sequencing was performed to compare the circRNA expression profiles between parental and OR cells. The function of the candidate circRNA was assessed through a series of in vitro and in vivo experiments, including cell viability assays, apoptosis analysis, and xenograft mouse models. Mechanistic investigations involved mass spectrometry, co-immunoprecipitation and western blotting to explore its protein-coding potential and downstream signaling pathways. RESULTS: We identified a novel circRNA, termed circTLL1, that was stably and significantly upregulated in OR-NSCLC cells. Functionally, overexpression of circTLL1 promoted osimertinib resistance, whereas its knockdown restored drug sensitivity both in vitro and in vivo. Mechanistically, we discovered that circTLL1 harbors an open reading frame (ORF) that is translated into a novel 90-amino-acid protein, which we designated circTLL1-90aa. Further investigation revealed that circTLL1-90aa directly interacts with and promotes the degradation of 5'-nucleotidase, cytosolic II (NT5C2), thereby uncoupling nucleotide metabolism from its normal regulatory constraints. The consequent downregulation of NT5C2 leads to elevated GTP levels and leading to the sustained activation of the downstream Ras/PI3K/AKT signaling pathway. CONCLUSION: Our findings unveil a previously unrecognized circRNA/micropeptide/metabolism cascade underlying osimertinib resistance. The identification of the circTLL1-90aa/NT5C2/Ras/PI3K axis not only expands the functional repertoire of the non-coding genome but also provides new insights into the complexity of drug resistance. Given its selective upregulation in resistant cells, circTLL1-90aa holds promise both as a predictive biomarker for treatment stratification and as an actionable therapeutic target, offering a novel strategy to overcome osimertinib resistance in NSCLC patients.

Pyrimidines↗

PubMed2026

Beyond species trees: pervasive gene flow limits phylogenomic resolution in the diversification of Juniperus from the Qinghai-Tibet Plateau.

Understanding how lineages diversify despite persistent ancestral polymorphism and recurrent gene flow remains a central challenge in evolutionary biology. Juniperus distributed across the Qinghai-Tibet Plateau provide an ideal system for addressing this question because repeated geological uplift and climatic oscillations have likely promoted cycles of lineage divergence, range shifts, and secondary contact. Here, we combined approximately 1.08 million genome-wide SNPs from 164 individuals representing thirteen Juniperus lineages with phylogenomic datasets comprising 3,381 nuclear single-copy genes and nearly complete plastomes. We detected extensive phylogenomic discordance and cytonuclear incongruence across genomic datasets. Topology weighting, coalescent simulations, quartet-based tests, and analyses of gene flow and reticulation collectively support the interpretation that these patterns were shaped by the combined effects of prolonged incomplete lineage sorting and gene flow during lineage diversification. Ecological niche analyses further provide a spatial and climatic context in which environmentally similar lineages may have had greater opportunities for secondary contact during historical range shifts. Collectively, our results reveal that the evolutionary history of Qinghai-Tibet Plateau Juniperus is characterized by reticulate diversification rather than strictly bifurcating evolution, and demonstrate how genome-wide discordance can provide biological insights into the evolutionary processes underlying lineage diversification.

Gene Flow↗

PubMed2026

Pervasive hybridization and introgression in Diervilleae (Caprifoliaceae).

Diervilleae (Caprifoliaceae) is a horticulturally important lineage with striking floral diversity and a long history of interspecific crossing, suggesting reticulate evolution. We integrated nuclear SNPs and whole plastome data to reconstruct a phylogenomic backbone for the tribe and to identify hybrids, cultivated accessions, and introgression among lineages. Nuclear and plastid phylogenies consistently recover Weigela and Diervilla as reciprocally monophyletic and resolve four major lineages within Weigela, providing a reproducible framework for revising sectional limits and species boundaries. Cultivated accessions form a well supported clade sister to W. florida and show predominantly W. florida ancestry while retaining contributions from multiple wild lineages, consistent with recurrent crossing, backcrossing, and selection. Analyses of wild populations reveal recurrent hybrids and enable plausible parental combinations to be inferred. Tests across the genome further indicate strong evidence for historical introgression across Diervilleae, with the strongest signals involving W. middendorffiana, W. maximowiczii, and Diervilla. Fossil evidence, divergence time estimation, and paleodistribution modelling together suggest range expansion during the Miocene and Pliocene followed by climate driven contraction, providing a spatiotemporal context for episodic contact, introgression, and the East Asia-North America disjunction.

Hybridization, Genetic↗

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

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

Introgression shapes the genomic conflict landscape of Malus, providing evidence for a reticulate backbone in a woody crop lineage.

Phylogenomic discordance is widespread across plants, but its evolutionary significance is often obscured when conflict is treated primarily as analytical noise rather than as evidence of underlying processes. In woody lineages in particular, incomplete lineage sorting, introgression, and genome duplication can interact over long timescales to produce complex genomic histories that are not adequately summarized by a strictly bifurcating tree. Here, we use Malus as a model woody genus to investigate how these processes structure conflict across a genus-scale, accession-based phylogenomic framework. Using broad taxon sampling, hundreds of nuclear loci, plastid genomes, and genome-wide SNP summaries, we reconstruct a robust nuclear backbone for sampled Malus lineages and evaluate where discordance is concentrated and which processes best explain it. Nuclear analyses resolve eight major clades, whereas conflict is non-random and localized to recurrent hotspots rather than evenly distributed across the tree. Cytonuclear discordance is similarly concentrated, especially around Clade H, represented by sampled accessions of M. tschonoskii, where localized plastid-nuclear disagreement is consistent with candidate plastid capture or organellar introgression. Multiple complementary analyses further indicate that the strongest conflict is not explained by ILS alone, but instead reflects lineage-structured introgression, while polyploid complexes represent additional localized sources of evolutionary complexity. Together, these results provide evidence for a reticulate genomic backbone in Malus and show how integrating nuclear, plastid, and genome-wide conflict analyses can help distinguish background discordance from process-specific signals in woody plant radiations. Several lineage-level reticulation hypotheses identified here should now be tested with broader population-level sampling and curated reference accessions.

Malus↗

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↗