Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “methylation profile”

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.

467 records · Page 9Linked to original sources
PubMed2026

Phenotypic, physiological and transcriptomic analysis of graded salt stress responses in Pyrus betulifolia Bunge and functional characterization of the hub gene PbSTY46.

Pyrus betulifolia Bunge is a salt‑tolerant rootstock for pear, but its salt‑tolerance mechanisms remain largely unknown. In this study, P. betulifolia seedlings were subjected to graded NaCl stress at concentrations of 0 (CK), 50 (T1), 100 (T2), and 200 (T3) mM. We integrated phenotypic observation, physiological assessment, transcriptomic profiling, and functional gene validation to systematically elucidate its salt tolerance mechanisms. Salt stress inhibited seedling growth and root traits in a concentration-dependent manner, and T3 caused the most severe damage. Osmotic solutes responded differentially: soluble sugars peaked under T2, while proline peaked under T3. Antioxidant enzymes showed tissue-specific biphasic responses and declined after prolonged T3 stress. Meanwhile, chlorophyll and photosynthesis decreased, whereas anthocyanin increased, indicating a metabolic shift from photosynthesis to photoprotection. Transcriptome analysis revealed distinct responses depending on stress intensity: mild stress induced membrane lipid remodeling, moderate stress activated circadian rhythm and hormone signaling, and severe stress enhanced phenylpropanoid biosynthesis and thiamine metabolism. Gene Set Enrichment Analysis (GSEA) further highlighted progressive enrichment of phenylpropanoid biosynthesis, heme binding, and oxidoreductase activity. Weighted Gene Co‑expression Network Analysis (WGCNA) identified a blue module significantly positively correlated with root traits, from which the hub gene PbSTY46 was identified. Functional validation via overexpression, loss‑of‑function mutants, and pharmacological interventions (MeJA/DIECA) confirmed that PbSTY46 acts through JA signaling to enhance antioxidant enzyme activities and thereby confer salt tolerance. Collectively, P. betulifolia adopts a "survival‑first" strategy that coordinates growth arrest, osmotic homeostasis, and ROS scavenging. These findings establish PbSTY46 as a key regulator that links JA signaling to antioxidant defense. Thus, PbSTY46 represents a promising candidate for marker‑assisted breeding of salt‑tolerant pear cultivars.

Salt Stress↗
PubMed2026

Uncovering hidden complexity in the Apis mellifera mitotranscriptome: a polyadenylation-centered perspective.

Mitochondrial transcription is gaining increasing attention as researchers seek to better understand the full coding potential of mitochondrial DNA (mtDNA). Emerging evidence suggests that mtDNA may encode additional elements beyond classical oxidative phosphorylation genes, pointing to a more complex transcriptional architecture than previously recognized. In this study, we explored the mitochondrial transcriptome of Apis mellifera (Insecta: Hymenoptera), with a particular focus on polyadenylation-associated features. Our analysis revealed that both sense and antisense transcripts undergo polyadenylation, although transcript abundance and poly(A) tail lengths varied markedly across mitochondrial genes. Several transcripts exhibited alternative isoforms, either extended or truncated, frequently including intergenic regions. These regions may represent functional non-coding elements or structural variants rather than conventional untranslated regions (UTRs). Interestingly, some transcripts also contained non-templated nucleotide additions particularly cytosine residues immediately upstream of the poly(A) tails. Monocistronic units that included portions of downstream intergenic regions were among the most abundantly represented, suggesting a possible regulatory role for these sequences. To experimentally validate our in silico findings, we performed RT-qPCR to assess relative gene expression and applied 3' RACE-PCR to define transcript boundaries. These approaches confirmed the presence of multiple transcript isoforms and supported the involvement of polyadenylation in shaping mitochondrial RNA diversity. Together, our findings reveal a previously underappreciated level of complexity in the A. mellifera mitochondrial transcriptome and highlight the potential regulatory significance of polyadenylation dynamics and intergenic region transcription.

Animals↗
PubMed2026

Integrated analysis uncovers exogenous induction and molecular regulation of erinacine A accumulation in Hericium erinaceus.

Erinacine A, a cyathane-type diterpenoid mainly from Hericium erinaceus mycelia, exhibits prominent neurotrophic and neuroprotective activities, making it a promising candidate for managing neurodegenerative diseases. However, its low abundance and unclear genetic regulatory mechanisms hinder its application as a nutraceutical. This study aimed to decipher its regulatory mechanisms and enhance production. Four exogenous inducers were screened, with salicylic acid (SA) and ergosterol (ERG) significantly increasing erinacine A content by 62.21% and 146.70% at 20 days, respectively. Transcriptome and WGCNA of inducer-treated sample identified darkorange and magenta modules associated with erinacine A biosynthesis, with the eri gene cluster enriched in the darkorange module and eriG and eriF as hub genes. Forward genetic analysis via QTL mapping of the HeD127 dikaryon population revealed significant phenotypic variation in erinacine A content (0.341-13.085 mg/g) and identified two loci (erA-1 and erA-2) explaining 18.63% of phenotypic variation. Integrating these forward and reverse genetic analyses revealed that salicylic acid and ergosterol synergistically regulate core carbon metabolic pathways to augment acetyl-CoA supply for the mevalonate pathway, suppressed competitive metabolism, enhanced diterpene skeleton construction and structural modification. These results deepen our understanding of the genetic and molecular basis governing accumulation of erinacine A, and facilitate its application in neuroprotective pharmaceuticals.

Diterpenes↗
PubMed2026

Morphological changes and transcriptomic insights into skeletal development of embryos and larvae of the sea urchin Strongylocentrotus intermedius.

To explore morphological features and molecular dynamics underlying skeletogenesis in the sea urchin Strongylocentrotus intermedius, we conducted combined morphological observation and comparative transcriptome analyses across representative embryonic and larval developmental stages. Morphological results showed that triradiate spicules first emerged at the gastrula stage. The 8-arm pluteus stage was identified as a key phase for skeletal remodeling, during which new three-radiate crystals transformed into complex stereoscopic ossicles including tube feet, spines and test plates. Transcriptomic data indicated that most differentially expressed genes (DEGs) were downregulated from the blastula to gastrula. The altered expression of basal metabolic genes and extracellular matrix genes including Colp2α and calm may be correlated with the linear mineralization of early spicules, which potentially reflects an energy adjustment pattern in developing larvae. During the transition from 6-arm to 8-arm pluteus, expression changes of calmodulin-like, Colp2α and SISin18G001660 suggest potential associations with regional calcium deposition and modifications of skeletal matrix properties. This work systematically characterizes morphological traits and transcriptional dynamics of skeletogenesis in S. intermedius. Its early spiculogenesis follows the conserved developmental pattern of echinoderms, while massive formation of stereoscopic ossicles occurs at the 8-arm pluteus stage. Stage-specific transcriptional changes across key larval skeletogenic stages are uncovered, offering transcriptomic resources for functional verification of skeletal regulatory genes.

Animals↗
PubMed2026

Transcriptomic analysis provides molecular insights into the innate immune defense of Mactra veneriformis against Vibrio alginolyticus infection.

Mactra veneriformis is an economically important bivalve mollusc in China, but its aquaculture is frequently threatened by Vibrio infections, particularly Vibrio alginolyticus. To investigate the molecular immune response of M. veneriformis to V. alginolyticus, we performed RNA-seq analysis of hepatopancreatic tissues collected at 48 h post-infection, the peak mortality time point, with PBS-injected individuals used as controls. Infection with V. alginolyticus caused severe histopathological damage in the hepatopancreas and resulted in a cumulative mortality of 53.3% over 14 d, compared with 3.3% in the control group. Transcriptomic analysis identified 2623 differentially expressed genes (DEGs), including 1585 significantly up-regulated genes and 1038 down-regulated genes. KEGG enrichment analysis demonstrated that DEGs were significantly enriched in immune related and metabolism pathways, including the JAK-STAT signaling pathway, RIG-I-like receptor (RLR) signaling pathway, and cytochrome P450 (CYP450) signaling pathway. Collectively, these findings revealed candidate immune related genes (tlr3, tlr5, myd88, nfkb1, il-17d, and ifi44l), a putative TLR-MyD88-NF-κB signaling axis, and KEGG signaling pathways, including JAK-STAT, RLR and CYP450, that may be involved in the innate immune response of M. veneriformis to V. alginolyticus infection. These results provide a transcriptomic basis for understanding host-pathogen interactions in this species and highlight candidate genes and pathways for future functional validation and potential application in disease-resistance breeding.

Animals↗
PubMed2026

Transcriptome analysis of the pectoral fin degeneration in half-smooth tongue sole (Cynoglossus semilaevis).

Appendage degeneration is a notable morphological feature of some teleosts with specialized benthic lifestyles. The half-smooth tongue sole (Cynoglossus semilaevis) undergoes severe pectoral fin regression during metamorphosis. However, the molecular basis underlying rapid pectoral fin degeneration remains unclear. Here, we performed time-series transcriptome sequencing on pectoral fins at pre-metamorphosis, metamorphosis peak and post-metamorphosis to characterize the molecular changes associated with pectoral fin degeneration. Transcriptional dynamics and functional enrichment showed that no significant enrichment of classical apoptosis-related transcriptional pathways was detected during pectoral fin degeneration. Instead, sustained downregulation of twist1b, identified as a transcriptomic candidate, together with significant upregulation of ssh1, coupled with enrichment of lysosome and ubiquitin-proteasome system (UPS) pathways, suggested enhanced tissue remodeling during pectoral fin degeneration. Temporal expression clustering revealed heterochronic misalignment in the developmental gene expression: upstream initiator tbx5 was upregulated at early metamorphosis, while downstream maintenance signal fgf10 decreased synchronously. Distal patterning gene hoxd12a exhibited premature expression and rapid decay, losing sustained late-phase expression. Moreover, transient elevation of gli3 during metamorphosis may contribute to restricted distal fin growth. We conclude that pectoral fin degeneration in C. semilaevis is associated with heterochronic disruption of developmental signaling and extensive tissue remodeling. This study provides transcriptomic insights into pectoral fin degeneration in tongue soles and establishes a basis for future functional studies of appendage reduction in teleosts.

Animals↗
PubMed2026

Metatranscriptomic analysis of viral sequences associated with Culex nigripalpus at an Alabama aquaculture site.

Mosquitoes associated with aquaculture habitats can harbor diverse viruses, yet the viromes of many locally abundant species remain poorly characterized. At an aquaculture-associated site in Auburn, Alabama, we surveyed mosquito populations and found Culex nigripalpus to be the dominant species collected. To characterize viruses associated with this mosquito, we performed RNA-seq on pooled female Cx. nigripalpus and compared complementary bioinformatic workflows for viral detection and genome recovery. One workflow removed host-associated reads by mapping to the closest available mosquito reference genome prior to assembly, whereas a second workflow used fully de novo assembly and viral database annotation. Additional protein-level filtering, cross-workflow comparison, and comparison of Trinity and rnaSPAdes assemblies were used to prioritize well-supported viral candidates. Across the original analyses, 16 submitted accessions corresponding to 12 collapsed virus/name groups were recovered, including Merida virus, Hubei mosquito virus 5, Zhejiang mosquito virus, Hubei virga-like virus 3, Rinkaby virus, Elemess virus, Qingnian mosquito virus, Serbia narna-like virus 2, XiangYun narna-levi-like virus 8, Ecclesville picorna-like virus, and baculovirus-like fragments. Several candidates were supported across multiple workflows, while others were recovered only under specific analytical conditions, indicating that candidate recovery was influenced by assembly and filtering choices. Selected viral contigs were independently supported by RT-PCR amplification. Overall, these results provide a first characterization of viral sequences associated with Cx. nigripalpus from an Alabama aquaculture-associated site and show that comparison across assembly and filtering strategies helped prioritize the most consistently supported viral candidates.

Animals↗
PubMed2026

Comparative transcriptomic analysis of the gills and hepatopancreas of freshwater-cultured Litopenaeus vannamei under chronic nitrite stress.

To investigate the differences in molecular responses between the gills and hepatopancreas of freshwater-cultured Litopenaeus vannamei under chronic nitrite stress, a 30-day chronic stress experiment was conducted with a control group and a stress group. Transcriptomic analysis of the gills and hepatopancreas was performed using Illumina sequencing; differentially expressed genes (DEGs) were identified, and GO, KEGG, GSEA, PPI, and RT-qPCR validation were carried out. The results showed that 196 DEGs (161 up-regulated and 35 down-regulated) were identified in the gills, and 287 DEGs (199 up-regulated and 88 down-regulated) in the hepatopancreas, with only 18 DEGs shared between the two tissues. DEGs in the gills were enriched in oxidoreductase activity, glycerophospholipid metabolism, and tyrosine metabolism; DEGs in the hepatopancreas were enriched in lipid transporter activity, phagosome, ECM-receptor interaction, and riboflavin metabolism. GSEA revealed significant suppression of the mTOR pathway in the gills and the Polycomb complex pathway in the hepatopancreas. PPI network analysis identified hub genes P5CS and eEF2 in the gills, and PER, TUBB1, SHMT, and TUBB4B in the hepatopancreas. RT-qPCR validation was consistent with the RNA-seq results (R2 = 0.764). This study indicates that, under chronic nitrite stress, the gill response is centered on redox regulation and inhibition of growth metabolism, whereas the hepatopancreas response primarily involves lipid transport, cytoskeletal remodeling, and phagosome activation. The two tissues synergistically adapt through fundamental biosynthetic and motor protein pathways. This research provides molecular evidence for deciphering the nitrite tolerance mechanisms in freshwater-cultured shrimp.

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

Comparative transcriptome analysis provides insights into dorso-ventral color pattern formation of Holothuria edulis.

Animal body color patterns are highly diverse and play critical roles in camouflage, intraspecific communication, and environmental adaptation. Holothuria edulis, an important echinoderm inhabiting tropical waters, exhibits a typical dorsoventral dichromatism. This unique body color difference represents a key phenotypic trait for its habitat adaptation; however, the core differential genes regulating this trait remain to be elucidated. In this study, comparative transcriptome sequencing was performed on the dorsal and ventral body wall tissues of H. edulis, leading to the identification of a number of differentially expressed genes (DEGs), followed by GO functional annotation and KEGG pathway enrichment analysis. GO enrichment analysis indicated that the DEGs were significantly enriched in functional categories such as extracellular region, peptidase inhibitor activity, and tetrapyrrole binding. KEGG pathway analysis further revealed significant enrichment of protein digestion and absorption, the TNF signaling pathway, and cholesterol metabolism. Notably, the pigmentation-related gene FMO2 was highly expressed in the dorsal body wall tissue, whereas cyp1a1, ZIC1, Slc7a11, WNT-1, and ADAMTS20 were highly expressed in the ventral body wall tissue. This study identified DEGs and enriched pathways associated with dorsoventral body color differences in H. edulis, providing new insights into the molecular regulatory mechanisms underlying body color pattern formation. From the perspective of aquaculture applications, body color is one of the important traits affecting the quality and market value of sea cucumber products. Elucidating the molecular mechanisms of body color variation can provide a scientific basis for molecular marker-assisted breeding of superior sea cucumber variety.

Animals↗
PubMed2026

Unveiling the molecular basis of gonadal development: Multi-omics uncovers sex-related genes and steroid pathways in Sinonovacula constricta.

The razor clam Sinonovacula constricta is an economically important cultured mollusk in China, but the molecular mechanism of its gonadal development and sexual differentiation remains unclear. This study integrated gonadal transcriptomic, proteomic, and metabolomic analysis to identify key sex-related molecules. Transcriptome analysis identified 2795 DELs and 6497 DEGs between sexes, including the sex-related genes Fem-1b, Fem-1c, GUCY1B2 and FAT4, as well as a regulatory network of 39 lncRNA-mRNA pairs involving Tektin-4, Ropporin-1, Histone H1, and FoxN4. Proteomic analysis revealed 3217 DEPs: Tektin family members, Ropporin-1 and Tssk proteins were upregulated in the testis, while histone H1 and FAT4 were upregulated in the ovary. Metabolomic analysis detected 409 DEMs, with uridine identified as a potential sex differential marker (upregulated in the ovary), and 23 gonadal development-related DEMs showed sex-specific upregulation. Integrative transcriptome-proteome analysis identified 1543 co-expressed DEGs/DEPs enriched in nucleosome assembly, oxidative phosphorylation, and carbon metabolism, including key sex-related genes AKAP14, Tektin/Tssk families, Histone H1, and FAT4. Transcriptome-metabolome integration identified 32 shared KEGG pathways (e.g., biosynthesis of unsaturated fatty acids, pyrimidine metabolism), while proteome-metabolome integration revealed 5 (positive ion) and 6 (negative ion) co-enriched pathways, with alanine, aspartate and glutamate metabolism and oxidative phosphorylation being functionally relevant to gonadal development. Collectively, these results reveal the molecular basis of gonadal development, highlight critical sex-related genes and steroid metabolic pathways, and provide valuable resources for future reproduction and breeding in S. constricta.

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

Transpulmonary proteomic gradient analysis in women with pulmonary arterial hypertension associated with systemic sclerosis.

This study investigated proteomic alterations in the pulmonary circulation of patients with pulmonary arterial hypertension associated with systemic sclerosis (PAH-SSc) by analyzing the transpulmonary protein gradient and comparing the proteomic profiles with systemic sclerosis (SSc) without PAH. Twenty women were included (10 PAH-SSc, 64.6 ± 10.8 years; 10 SSc, 62.8 ± 11.5 years). The transpulmonary gradient was defined as the difference in biomarker concentrations between wedge-position and pulmonary artery blood samples. Peptides were analysed using liquid chromatography-mass spectrometry, and differentially abundant proteins were identified with Proteome Discoverer. Protein-protein interaction networks were generated with STRING and visualized in Cytoscape. A total of 270 proteins were detected, with no significant transpulmonary gradient alterations. However, patients with PAH-SSc showed distinct proteomic profiles compared to SSc. Multivariate analysis identified 48 differentially abundant proteins in pulmonary artery plasma, with 15 overrepresented and 33 downregulated in PAH-SSc. Among these, the downregulation of transforming growth factor-beta-induced protein ig-h3 (TGFβI/ig-h3) points to a potential involvement of the TGF-β-related extracellular matrix remodelling pathway in PAH-SSc. However, further validation in larger and independent cohorts is required before its relevance as a biomarker or therapeutic target can be established. In conclusion, while no transpulmonary proteomic gradient was observed, the proteomic profiles of PAH-SSc and SSc were different. The profile in PAH-SSc was characterized by differences in immune response, lipid metabolism, and hemostatic proteins. SIGNIFICANCE: This study offers the first proteomic characterization of the transpulmonary gradient in PAH-SSc and SSc. Although no differences in the gradient were found, the pulmonary artery plasma proteome of PAH-SSc patients showed a distinct pattern compared to SSc. Several proteins associated with immune function, haemostasis, and cellular processes were altered, which may indicate specific pathophysiological features of PAH-SSc or suggest how lung dysfunction develops in SSc. Targeting dysregulated proteins like TGFβI/ig-h3 or addressing immune-coagulation imbalances may support future research studies. Overall, these findings refine the molecular profile of PAH-SSc and provide a basis for future large-scale studies aimed at clarifying disease mechanisms and identifying clinically relevant molecular signatures.

Humans↗
PubMed2026

Ribosome stalling position, spacing, and A-site occupancy impact translation and cotranslational mRNA decay in plants.

Ribosomes can pause during mRNA translation, but what causes pausing, how pauses affect protein production, and whether they trigger cotranslational mRNA decay are poorly understood in plants. Here, we investigate the causes and consequences of ribosome pausing in Arabidopsis and maize. This is accomplished by sizing, mapping, and quantifying footprints of individual ribosomes (monosomes) and closely spaced ribosome pairs (disomes) at single-codon resolution on open reading frames (ORFs). Ribosome footprinting was combined with 5'P-degradome-seq to examine the coincidence of pausing with cotranslational decay under control conditions and brief hypoxia in Arabidopsis. The data resolve two monosome conformations and three disome configurations. These include monosomes with a vacant or occupied A-site and disomes that have collided or are separated by one or two codons. Pausing is prevalent at initiation, termination, and di-Proline codons. Di-Proline pauses do not trigger cotranslational decay but appear important in cotranslational protein processing. Brief hypoxia induces stalling of A-site vacant ribosomes at Aspartate codons, often coinciding with 5'P peaks, indicating that rate-limiting decoding can trigger cotranslational mRNA decay. Notably, actively transcribed and translated hypoxia-response mRNAs accumulate 1- to 2-codon-separated disomes and are actively degraded. Comparative analysis of footprints in the two species reveals ribosome conformations and codon-specific pausing can be conserved or lineage-specific, as exemplified by pausing at di-Prolines and on Conserved Peptide upstream ORFs. In sum, the stalling of ribosomes at specific codons, coupled with ribosome A-site occupancy and disome spacing, modulates protein production and cotranslational mRNA decay in plants.

Ribosomes↗
PubMed2026

Mineral-driven molecular signatures of energy metabolism underpin sperm motility in buffalo.

The success of spermatogenesis depends on the interplay of various biomolecules that ultimately determine sperm quality. In this study, RNA-seq analysis of frozen-thawed buffalo sperm (n = 8) revealed the presence of 263 mineral-associated genes (>1 FPKM) in high (n = 4) and 181 in low motile (n = 4) sperm groups. Among these, 177 mineral-associated genes were commonly expressed between them, and the majority were upregulated (>1 fold), LOC102391588 (ncRNA; 37-fold), ZNF699 (26.5-fold), MYZAP (13-fold), etc., in the high motile group. The expression of selected mineral-associated genes was validated. The top enriched functions in commonly expressed genes were regulation of transcription by RNA polymerase II (FDR: 4.7 × 10⁻2; ZNF331, ZNF692, ZNF180, etc.), followed by spermatogenesis (FDR: 2.9 × 10⁻2; CALR3, ADAM18, ADAM29, etc.), proton transmembrane transport (FDR: 4.0 × 10⁻2; ATP6V0E1, ATP1A4, ATP6V1B2, etc.) and flagellated sperm motility (FDR: 2.3 × 10⁻1; CATSPERD, EFCAB6, CABS1, etc.). Additionally, the chromatin remodeling pathway (FDR: 3 × 10⁻2; PTP4A1, PPM1A, DUSP1, etc.) emerged as the most significant and may suggest that these minerals influence genome packaging and sperm functionality. Mineral-associated genes were predominantly associated with zinc (49%), followed by calcium (20%), phosphorus (10%), iron (5%), sodium (2%), potassium (2%), copper (1%) and other trace elements (11%). Although the current study uses frozen-thawed sperm, the findings indicate that mineral-associated genes are crucial for promoting membrane stability, energy production, motility and chromatin integrity, which may contribute to the superior fertilizing ability of sperm.

Animals↗
PubMed2026

An oxidative stress - and immunotherapy-related six-gene signature defines immune subtypes and predicts prognosis and immunotherapy response in hepatocellular carcinoma.

BACKGROUND: Oxidative stress and the tumor immune microenvironment jointly shape hepatocellular carcinoma (HCC) progression and response to immunotherapy, yet integrated biomarkers linking these processes are lacking. METHODS: Transcriptomic and clinical data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets were used to identify oxidative stress- and immunotherapyrelated differentially expressed genes (OSIRDEGs). Functional enrichment, weighted gene co-expression network analysis (WGCNA) and LASSO-Cox regression were used to construct a prognostic signature. Consensus clustering, TIDE, CIBERSORT and ssGSEA characterized immune phenotypes. Somatic mutation, copy-number and drug-response data were integrated to assess genomic alterations and drug sensitivity. Expression of model genes was validated by qRT-PCR and western blotting in HCC cell lines. RESULTS: We identified 24 OSIRDEGs enriched in cell-cycle and mitotic pathways. WGCNA intersection yielded 18 module genes, from which a six-gene signature (BUB1B, CDKN2A, CENPE, HMMR, PTTG1, SPP1) was derived. The signature robustly stratified patients into high- and low-risk groups with significantly different progression-free and disease-free survival in both TCGA-LIHC and GSE14520. Based on signature expression, two molecular subtypes were defined, exhibiting distinct survival, immune landscapes and predicted immunotherapy responsiveness. Model genes harbored recurrent alterations and showed significant correlations with anticancer agents. All six genes were upregulated at mRNA and protein levels in metastatic HCC cell lines versus normal hepatocytes. CONCLUSIONS: We systematically explored the landscape of OSIRDEGs in HCC, and proposed a validated six-gene signature that refines prognostic stratification, delineates immunerelevant HCC subtypes and highlights candidate biomarkers for therapeutic selection and mechanistic investigation.

Humans↗
PubMed2026

Transcriptomic insights into the coordinated regulation of signaling, apoptosis, immunity, and metabolism during Sinonovacula constricta larval metamorphosis.

Metamorphosis is a critical ontogenetic transition for marine bivalves, marking the shift from planktonic to benthic lifestyles, where successful transformation dictates survival. The razor clam Sinonovacula constricta is economically important; however, low larval metamorphosis rates remain a major bottleneck in seedling production. To elucidate the mechanisms governing this process, we performed a comparative transcriptome analysis of S. constricta larvae at pre- and post-metamorphosis stages using Illumina sequencing. A total of 3701 differentially expressed genes (DEGs) were identified, including 3254 up-regulated and 447 down-regulated genes. Functional annotation of the respective top 20 significantly up-regulated and down-regulated DEGs indicated their potential pivotal roles in signal transduction (e.g., up-regulated: CAV1, CHRNA2; down-regulated: APP, NOTCH1), cellular proliferation and differentiation (e.g., up-regulated: TUBA, EGF1; down-regulated: KIF23, TTC25), transcriptional and epigenetic regulation (e.g., up-regulated: NFIL3; down-regulated: OVO, HMX1), substance transport (e.g., up-regulated: LRP2, LRP1B; down-regulated: SLC51A, Slc33a1), substance metabolism (e.g., up-regulated: CPK3, CYP26A1; down-regulated: RDMT1, ADAC), immunomodulation (e.g., up-regulated: CPN2, CRISP2), and protein homeostasis (e.g., up-regulated: HSP27, NAS-27). Functional enrichment analysis further revealed that DEGs were significantly enriched in pathways related to signal transduction and developmental regulation (e.g., Ras, TNF), cell death and homeostasis (e.g., apoptosis), immune responses (e.g., Toll-like receptor), energy metabolism (e.g., lipid), cardiovascular related (e.g., Fluid shear stress), cell junction and architecture (e.g., Tight junction), and infectious disease (e.g., measles). These results suggest a synergistic interplay between signaling, apoptosis, immunity, and metabolism during S. constricta metamorphosis. This study advances our understanding of marine bivalve metamorphosis and offers candidate genes for further mechanistic studies.

Animals↗
PubMed2026

Integrative machine learning and transcriptomic analysis reveals molecular mechanisms underlying low survival rate in larval Chinese Bahaba (Bahaba taipingensis).

Chinese Bahaba (Bahaba taipingensis) is a Class I protected marine fish endemic to China. Low larvae survival during artificial breeding severely hinder population recovery. To investigate the molecular mechanism of high mortality in larval fish, this study performed RNA-seq on liver from naturally deceased (ND) and mass-dead (MD) individuals, combined with least absolute shrinkage and selection operator (LASSO) regression and random forest (RF) algorithms to screen for core signature genes. A total of 873 differentially expressed genes (DEGs) were identified, including 112 upregulated and 761 downregulated genes. GO and KEGG enrichment analyses revealed significant enrichment in amino acid metabolism disorders, one‑carbon folate pool impairment, PPAR signaling abnormalities, ECM-receptor interaction, focal adhesion pathway, indicating widespread metabolic suppression accompanied by extracellular matrix remodeling and signaling disturbances in the livers of MD fish. MAD pre-filtering combined with dual machine learning algorithms yielded 18 robust core signature genes, among which SLC38A4, MMP1, FADD, FKBP5, and APOB were consistently identified as high-frequency core genes by both algorithms. SLC38A4 exhibited the highest importance score in the RF model and was significantly downregulated, making it the primary molecule distinguishing ND from MD phenotypes. ROC curve analysis showed that both models achieved an AUC of 1.000 (95% CI lower bound: 0.610), confirming the precise discriminatory ability of the core genes. GSEA further demonstrated significant enrichment of this core gene set in ND samples. This study provides the first systematic elucidation of the molecular mechanisms underlying liver dysfunction in low survival rate B. taipingensis, characterized by amino acid transport impairment, metabolic reprogramming, and structural remodeling, offering theoretical foundations for health assessment, early mortality risk warning, and artificial breeding conservation of this species.

Animals↗