SEARCH · Search PubMed

Results for “Differentially expressed genes”

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.

680 records · Page 5Linked to original sources

PubMed2026

Targeting ncRNA control networks with engineered exosomes to overcome therapy resistance in thyroid cancer.

Papillary thyroid cancer (PTC) is the most prevalent endocrine malignancy, accounting for over 90% of thyroid cancers. While differentiated thyroid cancers (DTCs) typically have favorable outcomes, a significant subset progresses to radioactive iodine-refractory (RAIR) disease, characterized by impaired iodine uptake and a 10-year survival rate below 10%. Genetic alterations and dysregulated signaling pathways underlie this transition. Non-coding RNAs (ncRNAs), including microRNAs (miRNAs), circular RNAs (circRNAs), and long non-coding RNAs (lncRNAs), play critical regulatory roles in tumor biology and may be transported via exosomes, facilitating intercellular communication and contributing to RAIR-PTC. This systematic review, conducted according to PRISMA 2020 guidelines, evaluated the role of exosomal ncRNAs in RAIR-PTC. A comprehensive search of PubMed, PubMed Central, and Google Scholar identified studies published within the past 15 years in English. Following stringent quality appraisal, studies with a non-bias score above 40% were included. Of 961 identified publications, 96 high-quality studies met inclusion criteria. Evidence indicates that therapy resistance in RAIR-PTC is driven by convergent ncRNA regulatory networks that suppress sodium-iodide symporter (NIS) expression and activate oncogenic pathways, most notably MAPK, PI3K/AKT/mTOR, and Wnt/β-catenin signaling. Multiple ncRNAs converge on key regulatory nodes, forming redundant circuits that sustain dedifferentiation, metabolic adaptation, and impaired iodide transport. Several consistently dysregulated ncRNAs directly or indirectly regulate NIS expression and trafficking, highlighting actionable targets. Exosomes emerge as biologically compatible, programmable delivery vehicles capable of transporting therapeutic ncRNA payloads independent of endogenous packaging mechanisms. These findings support a precision therapeutic paradigm in which engineered exosomes reprogram ncRNA networks to restore iodine-handling pathways and overcome therapy resistance in RAIR-PTC.

Humans↗

PubMed2026

Addressing lignin composition and content via Arabidopsis arogenate dehydratase knockout and over-expression genotypes.

Following the down-selection of 14 Arabidopsis thaliana arogenate dehydratase (ADT) knockout and over-expression (OE) genotypes, the most highly contrasting quadruple knockout adt3/4/5/6 and ADT OE genotypes were subjected to proteomics, metabolomics, and scanning electron microscopy (SEM) analyses as needed, with results compared to Columbia wild-type (WT). The basal adt3/4/5/6 stem cross-sections, ∼70% lignin content reduced, exhibited buckled vessel cell walls and partially detached xylary fibers, in contrast to WT and ADT4m/5 m OE genotypes that did not. Anatomical defects primarily resulted from guaiacyl lignin level reductions in vessels with concomitant increased stem syringyl:guaiacyl (S/G) ratios. Phenylpropanoid and various upstream shikimate-chorismate pathway enzyme abundances, as well as specific monolignol oxidases (laccases/peroxidases), generally increased in adt3/4/5/6 at different stem and rosette leaf growth/development stages, relative to WT. Opposite effects were largely observed with the ADT5m OE genotype. By contrast, flavonoid and glucosinolate pathway enzyme amounts varied. Such enzyme abundance increases were overall unproductive as adt3/4/5/6 was unable to restore WT, ADT4 OE, ADT5 OE, ADT5m OE, and ADT4m/5 m OE secondary metabolite (lignin, phenylpropanoid, lignan, flavonoid, phenolic acid, and glucosinolate) levels. Conversely, ADT OE genotypes did not significantly increase programmed lignin levels or alter S/G compositions. In sum, proteomics analyses of adt3/4/5/6 and adt5 'perceived' that lignin and low molecular weight secondary metabolite amounts were not at 'programmed' levels as for WT and ADT OE genotypes but observed increases in relevant pathway protein abundances were futile. Notably though, proteomics analyses did not lead to predicting that lignin and associated biochemical pathways would have reduced metabolite levels, relative to WT and ADT OE genotypes. Genotype adt3/4/5/6, possibly the highest lignin level reduced genotype reported, did not utilize other phenolics to compensate. By contrast, the differential temporal and spatial deposition of cell wall oxidases again indicate the exquisite control over lignin deposition, and our lack of knowledge of precise lignin structure and assembly in subcellular regions of the lignified cell walls.

Lignin↗

PubMed2026

Comparative Transcriptomics Reveals Shared Downstream Pathways in Craniofacial Pathology.

Treacher Collins syndrome and Nager syndrome are craniofacial developmental disorders caused by defects in ribosome biogenesis and RNA splicing, respectively, yet they exhibit overlapping abnormalities affecting neural crest cell-derived craniofacial structures. To investigate shared downstream pathogenic mechanisms, we performed a comparative transcriptomic analysis of zebrafish polr1c and sf3b4 mutant models from our previous studies. Comparative analysis identified 17 shared differentially expressed genes (DEGs) between polr1c and sf3b4 mutants, with the majority of shared genes dysregulated in the same direction, indicating a coordinated rather than random transcriptional response. Gene ontology analysis identified ATP-dependent protein folding chaperone activity as the only shared molecular function, driven in part by upregulation of hsp90aa1.2, indicating a common proteostasis response. Because chaperone activity is linked to extracellular matrix (ECM) protein processing, we cross-referenced DEGs from both mutants against the curated zebrafish matrisome. Three of the 17 shared DEGs (serpinh1b, il11a, and lepa) were matrisome-associated and upregulated in both mutants. Serpinh1b, a collagen-specific chaperone, was strongly expressed in craniofacial cartilage and mesenchymal populations during pharyngeal arch development and exhibited nearly identical fold changes in both mutants. Il11a is of particular interest because its receptor, IL11RA, is known to be associated with human craniosynostosis, suggesting potential relevance to craniofacial development. Together, it is possible to hypothesize that shared chaperone-associated transcriptional changes, together with altered ECM-related gene expression, may contribute to polr1c- and sf3b4-associated craniofacial disorders, warranting further functional validation.

Extracellular Matrix↗

PubMed2026

Integrative analysis of transcriptome and DNA methylome dynamics during caudal fin regeneration in silver pomfret (Pampus argenteus).

Caudal fin regeneration in teleost fish is a complex, multi-stage process involving coordinated molecular and cellular changes. While the role of epigenetic regulation particularly DNA methylation has been studied in model freshwater species such as zebrafish, its contribution to regeneration in marine teleosts remains largely unexplored. In this study, we integrated transcriptomic and DNA methylomic data to characterize the temporal dynamics of gene expression and methylation during caudal fin regeneration in the silver pomfret (Pampus argenteus). Using RNA-sequencing and reduced representation bisulfite sequencing (RRBS) at three biologically critical time points 1, 3, and 7 days post-amputation (dpa), we characterized the spatiotemporal molecular landscape of caudal fin regeneration. These time points capture the key transitional phases of wound healing and inflammation (1 dpa), blastema formation and progenitor proliferation (3 dpa), and regenerative outgrowth with tissue remodeling (7 dpa), enabling robust detection of the major molecular programs underlying epimorphic regeneration. Concurrently, CG-methylome analysis identified thousands of dynamically changing differentially methylated regions (DMRs). A strong global inverse correlation was observed between promoter methylation and gene expression. Integrative analysis pinpointed key regeneration genes (fgf20a, msxb, sox9b) whose expression was associated with dynamic methylation changes in their promoters or gene bodies. We conclude that DNA methylation is a dynamic and key regulatory layer that acts in concert with transcriptional reprogramming to coordinate tissue regeneration, providing new insights into the epigenetic mechanisms underlying complex regenerative processes in teleosts.

Animals↗

PubMed2026

Comparative analyses of olfactory receptor repertoires in Schizothorax fish based on the chromosome-level genomes: Implications for regulatory roles of dietary differentiation and ploidy variation.

The olfactory receptor (OR) genes constitute the molecular basis of fish olfaction, mediating survival behaviors and environmental adaptation while coevolving with habitat-driven evolution. Schizothorax, a cyprinid genus endemic to the Qinghai-Tibetan Plateau, exhibits remarkable dietary divergence and ploidy variation in response to plateau environmental changes, which presumably facilitates the adaptive evolution of OR genes. However, the evolutionary patterns of OR genes associated with trophic divergence and ploidy variation in this genus remain unclear. In this study, three species were selected: the herbivorous diploid S. macropogon, the carnivorous diploid S. lantsangensis, and the herbivorous tetraploid S. curvilabiatus. S. macropogon possessed 142 OR genes (92.25% functional), primarily located on chromosomes 14 and 24, with the fewest sequence clusters. Such compact gene repertoire and highly overlapping chromosomal clusters indicated specialization for a herbivorous olfactory niche. S. lantsangensis contained 127 OR genes (93.70% functional), concentrated on chromosomes 4 and 5, with fewer sequence clusters and a scattered distribution, reflecting evolution of OR genes under carnivorous feeding habits. The herbivorous tetraploid S. curvilabiatus exhibited striking features: 316 OR genes (94.30% functional), the most subfamilies, unique ε and κ OR subfamilies, and species-specific motifs. These characteristics revealed that ploidy, rather than herbivory, dominated OR gene evolution. In conclusion, dietary differentiation and ploidy variation together drove olfactory adaptive evolution in Schizothorax, providing new insights into vertebrate OR gene ecological adaptation.

Animals↗

PubMed2026

Integrative genomic and transcriptomic analyses identify key regulators of skin pigmentation in Larimichthys crocea.

The yellow body coloration of large yellow croaker (Larimichthys crocea) constitutes a crucial economic trait, yet its underlying genetic regulatory mechanisms remain poorly understood. This study systematically elucidated the molecular basis of body color variation by integrating genome resequencing and skin transcriptome analyses, combined with the contextual analysis of key pigmentation-related genes and phenotypic histological validation. 200 phenotyped individuals (including yellow-selected lines, F1 progeny, and normal control groups, all derived from a well-characterized aquaculture stock) identified 39 significantly associated SNPs (-log₁₀(P) ≥ 6), mapping to multiple candidate genes. These genes were significantly enriched in pathways related to pigment deposition (GO:0033059), melanosome organization (GO:0032438), melanogenesis, and tyrosine metabolism. Cross-developmental stage transcriptome analysis revealed 2395 differentially expressed genes (DEGs). Multi-omics integration identified eight overlapping candidate genes, including tyrp1, slc45a2, oca2, and dgat2, among which tyrp1 was prioritized for in-depth validation based on its core regulatory role in eumelanin synthesis, significant SNP association signal, and consistent downregulation in transcriptomic data. Experimental validation demonstrated that the g.895C > T mutation in exon 2 of tyrp1b was strongly significantly associated with the yellow phenotype: the frequency of mutant genotypes (TT/CT) reached 92.86%in the yellow-selected group, whereas the control group exclusively exhibited the wild-type genotype (CC). qPCR confirmed significantly downregulated tyrp1b expression in the skin of yellow individuals, consistent with the transcriptome trend. Histological and stereomicroscopic observations of skin tissues further validated the physiological basis of the yellow phenotype, revealing a significant reduction in melanophore number and abnormal melanosome morphology in yellow-phenotype individuals, accompanied by increased xanthophore density. These results suggest that tyrp1b mutation is strongly associated with the yellow phenotype. However, the presence of a wild-type CC individual in the yellow group indicates that this mutation is not strictly required for yellow coloration, suggesting that other genetic or environmental factors may also contribute to the phenotype, Additionally, downregulation of the carotenoid metabolism gene bco2 coupled with upregulation of xdh, together with the functional changes of slc45a2 and oca2, may synergistically promote xanthophore pigment deposition, contributing to the yellow phenotype. As melanin synthesis in large yellow croaker relies on the conserved tyrosinase pathway and transporter proteins, mutations in associated genes (tyrp1b, slc45a2, oca2) represent a primary underlying cause for the loss of melanin-based coloration and transition to a yellow phenotype in L. crocea. These findings provide key molecular targets and a theoretical foundation for molecular breeding of body color in this species, and also enrich the understanding of xanthism regulatory mechanisms in teleosts.

Animals↗

PubMed2026

Genome-wide scans reveal candidate genes associated with wing morph differentiation in Tetrix japonica.

Wing dimorphism is an important dispersal-related trait in insects, but its genomic basis remains poorly understood in pygmy grasshoppers. Here, we integrated genome-wide single-nucleotide polymorphism (SNP) analyses, population structure inference, selection scans, and functional annotation to investigate genomic differentiation between long- and short-winged Tetrix japonica. Principal component analysis (PCA), ADMIXTURE, and phylogenetic analyses revealed weak genome-wide separation between morphs, indicating differentiation on a largely shared genetic background. Genome-wide scans based on the fixation index (FST), nucleotide diversity ratios, and Tajima's D, using 50-kb non-overlapping windows and empirical top-5% outlier thresholds, identified multiple candidate regions across seven chromosomes. The broader long- and short-winged candidate sets spanned 9.35 Mb and 9.37 Mb and directly overlapped 82 and 77 genes, respectively. Candidate genes were associated with signaling/hormone regulation, membrane transport, metabolism, cytoskeletal organization, extracellular matrix structure, and development. Short-winged candidate genes were significantly enriched for ABC-type transporter activity and ATP hydrolysis activity. Because all individuals originated from a single laboratory-maintained population with weak genome-wide structure, these regions should be regarded as candidate loci from a screening-stage analysis that require validation in independent populations and by functional assays, rather than as confirmed targets of selection.

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

Dynamic evolution of chaperone-mediated autophagy is associated with tumor microenvironment remodeling and prognostic stratification in lung adenocarcinoma: insights from single-cell transcriptomics, ensemble machine learning, and experimental validation.

BACKGROUND: Lung adenocarcinoma (LUAD) shows prognostic heterogeneity, and tumor-node-metastasis (TNM) staging is limited for individualized management. Chaperone-mediated autophagy (CMA) maintains proteostasis, but its role during adenocarcinoma in situ (AIS)-minimally invasive adenocarcinoma (MIA)-invasive adenocarcinoma (IAC) progression remains unclear. METHODS: Single-cell RNA sequencing (scRNA-seq) data from GSE189357 and bulk transcriptomes from The Cancer Genome Atlas (TCGA)-LUAD and Gene Expression Omnibus (GEO) cohorts were integrated. CMA activity, cell-cell communication, weighted gene co-expression network analysis (WGCNA), tumor-normal differential expression, machine-learning survival modeling, tumor microenvironment (TME) features, drug sensitivity, and EPC1 function were analyzed. RESULTS: CMA-high tumor epithelial cells increased from AIS (58.1%) to MIA (65.7%) but declined in IAC (44.4%; p < 0.001). CMA-low cells preferentially received fibroblast-derived extracellular matrix cues. A CMA-negatively correlated module identified 69 core genes. Random survival forest (RSF) performed best among 117 machine-learning combinations (mean concordance index > 0.873). High-risk patients had worse survival across cohorts, and the risk score was independently associated with overall survival (hazard ratio = 16.013, 95% confidence interval: 9.579-26.768, p < 0.001). High-risk tumors showed proliferative activation and M0 macrophage enrichment, whereas low-risk tumors showed stronger immune-related signaling. EPC1 overexpression suppressed malignant phenotypes in A549 cells. CONCLUSION: CMA dynamics are associated with stromal and immune remodeling during LUAD progression. A CMA-based model provides robust prognostic stratification and may offer a basis for future TME-guided studies.

Chaperone-mediated autophagy↗

PubMed2026

Occupationally relevant vibrations and the brain: frequency-dependent proteomics signatures in a rat model.

INTRODUCTION: Occupational exposure to whole-body vibration (WBV), particularly in agricultural environments, has been associated with adverse cognitive and physiological effects. This study examined the neurophysiological impact of WBV in a rat model at 4 Hz and 30 Hz, frequencies representative of off-road and on-road vehicle operation. METHODOLOGY: Forty-four Sprague-Dawley rats were assigned to control (0 Hz), low-frequency (4 Hz), or high-frequency (30 Hz) vibration conditions. After three days of exposure, brain tissues were collected and analyzed using mass spectrometry-based proteomics to identify differentially expressed proteins. RESULTS: Proteomic profiling revealed distinct, frequency-dependent alterations in brain protein expression. Compared with controls, 32 cognition-related proteins were differentially regulated at 4 Hz and 29 at 30 Hz, with 13 differing between the two vibration conditions. Principal component analysis showed clear separation among groups, indicating unique proteomic signatures for each exposure frequency. Functional enrichment and protein-protein interaction analyses demonstrated involvement of synaptic plasticity, cytoskeletal organization, calcium regulation, and neurotransmitter release. Exposure to 4 Hz was associated with the upregulation of proteins involved in calcium homeostasis and synaptic integrity, suggesting potential disruption of cognitive processes. In contrast, 30 Hz increased the expression of proteins related to axonal guidance and neuroprotection, indicating a less clearly adverse response that may reflect adaptive or potentially beneficial effects. DISCUSSION: These findings provide new insight into biological mechanisms underlying WBV-induced cognitive changes and underscore the importance of vibration frequency in shaping neurophysiological outcomes. They also establish a foundation for future studies integrating proteomics with behavioural assessments in animals and humans.

Animals↗

PubMed2026

Blood Bile Acids for Inflammatory Bowel Disease Diagnosis and Disease Activity Assessment: A Metabolomics Meta-Analysis.

Alterations in circulating bile acids (BAs) have been reported in inflammatory bowel disease (IBD), but the consistency of these changes across clinically relevant comparisons remains unclear. Our goal was to investigate systemic BA alterations in IBD using a metabolomics meta-analysis with an exploratory analysis of BA-related gene expression as a supporting context. A systematic review and meta-analysis of 28 metabolomics studies examined blood BA profiles associated with IBD, IBD diagnosis, and disease activity assessment. Univariate analysis and logistic regression modeling of two independent IBD cohorts explored the blood BA-related genes and IBD. Across 28 studies that comprised 5056 IBD patients, 1721 healthy controls, and 314 non-IBD patients, 131 BAs were reported. Eight predefined clinical comparisons were eligible for the meta-analysis. Lower secondary BA levels were consistently observed in IBD patients compared with controls, between UC and CD, and in active versus remission patients. Deoxycholic acid, glycodeoxycholic acid, and taurodeoxycholic acid were frequently decreased, whereas glycocholic acid was increased in certain comparisons. Transcriptomics analyses revealed differential expression of several BA-related genes in blood, including SLC51A, ABCB4, and ACOT8, across the comparisons. Our findings identify consistent circulating BA alterations in IBD and highlight the relevance of blood BA for future biomarker research in the diagnosis and disease activity assessment.

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

Evolutionary diversification of invertase paralogs couples carbon metabolism and sexual reproduction in fission yeasts.

Dynamic patterns of gene gain and loss play a major role in the diversification of eukaryotes, reflecting adaptation to a broad range of ecological contexts. Reconstructing the evolutionary history of genes provides a powerful framework for understanding how functional innovation shapes life-history traits. Here we report a comprehensive analysis of gene gain and loss across the fission yeast clade, whose evolutionary trajectory remains elusive. Reductive evolution of metabolic genes is a major contributor to species diversification, as observed in other fungal taxa. Notably, we uncovered an evolutionary scenario in which an ancestral gene duplication was followed by lineage-specific loss of one or the other paralog, except in S. pombe, which retained both. We demonstrate that these paralogs encode catalytically-active invertases, named Inv1 and Inv2, with distinct enzymatic properties, localization, regulation, and physiological roles. Inv1 is a secreted enzyme subject to glucose catabolite repression and is the sole invertase required for sucrose assimilation, resembling canonical yeast invertases. In contrast, Inv2 is intracellular, constitutively expressed, and required for inducing sexual differentiation in response to nutrient availability. Overall, these findings reveal an unexpected role for carbon metabolism in modulating the haploid-diploid cycle of fission yeasts, suggesting that diversification of core metabolic functions may contribute to adaptation to environments with distinct sugar compositions.

Evolution↗

PubMed2026

Gene-environment interaction between perinatal oxytocin exposure and Pten mutation shapes epigenetic reprogramming of oxytocin signaling and behavior in mice.

Synthetic oxytocin (Pitocin) is the most commonly used pharmacologic agent for induction and augmentation of labor. Beyond its uterotonic effects, oxytocin plays a critical role in neurodevelopment and social behavior. Dysregulated oxytocin signaling has been implicated in autism spectrum disorder (ASD), raising concern that perinatal exposure to exogenous oxytocin may have lasting neurodevelopmental consequences. This study aimed to determine whether offspring harboring a genetic predisposition for ASD are differentially impacted by perinatal oxytocin exposures, with a focus on long-term oxytocin signaling and autism-like behavior. Pregnant mice carrying offspring with heterozygous mutations in phosphatase and tensin homolog deleted on chromosome ten (Pten), a well-established monogenic risk factor for ASD, received continuous oxytocin versus phosphate-buffered saline (PBS) control via micro-osmotic pumps during late gestation. Wild-type (WT) offspring exposed to each treatment served as a secondary control. Adult offspring were assessed for oxytocin receptor (Oxtr) methylation in the frontal cortex and hippocampus, oxytocin expression in the hypothalamus, serum oxytocin levels, and were subject to a battery of social and anxiety-related behavior tests. Perinatal oxytocin exposure produced genotype-dependent effects in offspring. Epigenetic analyses revealed bidirectional remodeling of Oxtr methylation in the frontal cortex and hippocampus, with increased exon 1 methylation in WT mice and decreased methylation in Pten-mutant mice, resulting in significant genotype-treatment interactions. Hypothalamic oxytocin expression increased following treatment regardless of genotype, though baseline levels were higher in Pten-mutant mice. Neither oxytocin treatment nor genotype impacted long-term serum oxytocin levels. Behavioral outcomes were modest but context-specific: repetitive behaviors and cognition performance were unchanged, but oxytocin-treated Pten-mutant mice exhibited increased anxiety-like behavior alongside improved social memory. In contrast, oxytocin-treated WT mice showed reduced social novelty preference. Exploratory analyses suggested potential sex-dependent trends. Our findings support a model in which genetic susceptibility shapes the epigenetic encoding of early-life hormonal signals, thereby recalibrating oxytocin system function and downstream behavioral outcomes. Together, these data highlight the context-dependent effects of perinatal oxytocin exposure and argue against uniformly beneficial or detrimental effects, emphasizing the importance of gene-environment interactions in neurodevelopmental trajectories.

Animals↗

PubMed2026

The Role of Small Segmental Duplications in Generating Identical Isoforms Through Alternative Splicing Sites.

Alternative splicing plays a crucial role in expanding proteomic diversity but can also generate identical isoforms under certain conditions. While mutually exclusive splicing of tandem exons has occasionally been reported to produce identical isoforms, the extent to which other splicing events contribute to this phenomenon remains unclear. In this study, we demonstrate that alternative 5' and 3' splice site selection can also lead to the formation of identical isoforms, providing an additional type of splicing event for functional redundancy in transcriptomes. To address this, we analyzed reference genome annotations from 15 plant species, including Arabidopsis thaliana and wheat (Triticum aestivum), obtained from the RefSeq database. Identical isoforms were computationally defined as transcripts with distinct exon-intron structures but identical coding sequences. Our analysis reveals that the majority of alternative 5' and 3' fragments originate from small segmental duplications, suggesting that sequence repetition within gene regions facilitates the emergence of such splicing patterns. We also observed differences in the annotated 5' UTRs of some identical isoforms. However, since the alternative splicing sites themselves were not located within UTRs, these differences may reflect annotation uncertainty rather than genuine AS-derived variation. Given that UTR predictions in reference databases are not always precise, such observations should be interpreted cautiously. Expression analysis using an isoform-specific k-mer approach confirmed that identical isoforms can be differentially regulated. These findings suggest that, beyond expanding protein diversity, alternative splicing can also generate redundant isoforms that are differentially expressed at the RNA level, indicating potential regulatory roles. By elucidating the structural and regulatory factors contributing to the formation and retention of identical isoforms, our study provides new insights into the evolutionary and functional significance of alternative splicing in plants.

Alternative Splicing↗

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↗

PubMed2026

Genomic characterization of a hypervirulent Aeromonas veronii NN0115 from Nile tilapia and head kidney transcriptome of infected fish reveals B-cell-dominated immune response with specific immunoglobulin downregulation.

Aeromonas veronii is a pathogen of multiple fish species, yet systematic understanding of its infection in Nile tilapia (Oreochromis niloticus) remains limited. A dominant strain, NN0115, was isolated from a natural outbreak and identified as A. veronii by 16S rRNA and whole-genome average nucleotide identity (ANI, 96.33%). Experimental infection revealed high virulence (LD50 = 3.41 × 106 CFU/mL, equivalent to 8.53 × 104 CFU/fish). The genome is 4.58 Mb (58.57% GC) and encodes 4216 proteins. Virulence factor analysis identified 1253 genes, dominated by motility-related (264) and immune modulation (208) factors. Genomic island GI2 harbors 7 virulence genes and two dual-function resistance-virulence genes. The strain is resistant to 9 of 25 agents tested but carries three RND efflux pump genes whose predicted resistance was not phenotypically observed. The head kidney transcriptome of tilapia at 24 h post-bacterial infection identified 773 differentially expressed genes; among them, 57 were immunoglobulin (Ig) genes, and 56 were down-regulated. Integration of published single-cell transcriptomic data showed that non-Ig B-cell marker genes were down-regulated by 32%, whereas Ig genes were reduced by 63%, indicating selective transcriptional suppression of Ig genes rather than a general decrease in B-cell transcriptional activity. Together, this study provides a comprehensive characterization of a highly virulent A. veronii from Nile tilapia and reveals that selective downregulation of B-cell Ig genes is the dominant transcriptional feature of the host head kidney response.

Animals↗

PubMed2026

Transcriptomic insights into exogenous fatty acid-enhanced halotolerance in Zygosaccharomyces rouxii.

BACKGROUND: High salinity restricts microbial growth during brine-based food fermentation. Although exogenous unsaturated fatty acids improve the salt tolerance of Zygosaccharomyces rouxii, the associated transcriptional mechanisms remain unclear. This study investigated the transcriptomic response of Z. rouxii CGMCC 3791 to palmitoleic acid (C16:1) under high salt conditions. RESULTS: Cells were cultured in yeast extract peptone dextrose (YPD) containing 120 g L-1 NaCl, with or without 20 μg mL-1 C16:1. They were analyzed by RNA sequencing. Principal component analysis clearly separated the two treatments. Using q < 0.05 and |log2 fold change| > 1, 23 differentially expressed genes were identified - three upregulated and 20 downregulated. INO1, MLS1, POX1, MEP2, and SOD5 were among the major responsive genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses indicated that peroxisome-associated functions, lipid metabolism, oxidative stress responses, nitrogen utilization, and mitogen-activated protein kinase (MAPK) signaling were the principal C16:1-responsive processes. CONCLUSION: Exogenous C16:1 elicited a focused transcriptional adjustment rather than broad transcriptome-wide reprogramming in salt-stressed Z. rouxii. The results indicated that peroxisome-linked lipid processes and redox regulation were candidate mechanisms underlying fatty-acid-associated halotolerance and provided targets for improving the robustness of high-salt food fermentation. © 2026 Society of Chemical Industry.

Zygosaccharomyces rouxii↗