Search PubMed⌕ Search

SEARCH · Search PubMed

Search Search PubMed

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.

66 records · Page 3Linked to original sources
PubMed2026

Transcriptomic responses to developmental temperature in two field-collected Spodoptera exigua populations from Korea.

The beet armyworm, Spodoptera exigua, is a polyphagous insect whose development and seasonal occurrence are strongly influenced by temperature. However, transcriptomic responses to developmental thermal regimes remain insufficiently characterized in field-collected populations. In this study, we compared two Korean field-collected populations of S. exigua: a Haenam population collected in May and initially maintained at 15 ± 1 °C (HN), and a Jeju population collected in July and initially maintained at 27 ± 1 °C (JJ). F1 larvae from each population were reared under three fluctuating developmental temperature regimes: low (15-21 °C), middle (21-27 °C), and high (27-33 °C), followed by RNA-seq analysis. Differential expression analysis revealed population-associated variation in transcriptomic responses across developmental temperatures. HN exhibited a larger number of differentially expressed genes under the high-temperature regime, suggesting stronger transcriptomic sensitivity to elevated developmental temperature. Functional enrichment analyses identified population-associated differences in pathways related to heat response, oxidative metabolism, cytoskeletal organization, cuticle-associated processes, lipid metabolism, and immune-related functions. In JJ, heat-response and cuticle-related expression patterns were more prominent under warmer developmental conditions, whereas HN showed broader changes in stress- and metabolism-associated pathways under high temperature. Overall, this study provides a comparative transcriptomic analysis of two field-collected S. exigua populations under different developmental temperature regimes and identifies RNA-seq-based molecular response patterns associated with population-specific thermal response profiles.

Animals↗
PubMed2026

Genome-guided stage- and tissue-resolved transcriptome analysis of Serrodes campana identifies sex-biased antennal expression and candidate chemosensory-related genes.

Serrodes campana is an erebid moth of ecological and forestry relevance; its larvae are mainly associated with the soapberry tree, Sapindus mukorossi, whereas adults exhibit fruit-piercing behavior. However, stage- and tissue-resolved transcriptomic resources for this species remain limited. Here, using a chromosome-level reference genome, we performed a genome-guided transcriptome analysis of S. campana based on 12 RNA-seq libraries representing major developmental stages and key adult tissues. Global transcriptomic analyses revealed pronounced transcriptional differentiation across developmental stages and tissue types. Tissue-enriched gene sets and functional enrichment analyses identified distinct molecular signatures associated with developmental, sensory, and pheromone-associated tissues. Comparative analysis of female and male antennae further revealed sex-biased expression of several candidate chemosensory-related genes. Among 153 curated chemosensory-related candidate genes, most odorant receptor genes showed strong antennal enrichment, whereas other major chemosensory gene families displayed broader but still tissue-preferential expression patterns. In addition, an exploratory comparison of female terminal abdominal gland tissue and male terminal abdominal coremata revealed divergent expression profiles and highlighted candidate genes potentially associated with pheromone-related physiology, reproduction, and tissue-specific signaling. Together, this study provides the first genome-guided stage- and tissue-resolved transcriptomic resource for S. campana and offers a useful foundation for future studies of chemosensory detection, sex-biased gene expression, and pheromone-associated biology in this species.

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

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

Transcriptomic insights into thermal stress reveal physiological trade-off between thermal stress adaptation and reproductive investment in Spodoptera litura.

Spodoptera litura, a highly polyphagous lepidopteran pest, poses a major threat to agricultural productivity due to its remarkable adaptability to diverse environmental conditions. Although heat stress is known to trigger transcriptional reprogramming in insects, the molecular mechanisms underlying thermal stress responses in S. litura remain poorly understood. In the present study, fourth-instar larvae were exposed to acute heat stress (44 °C) and compared with control conditions (27 ± 1 °C) to investigate heat-induced transcriptional alterations affecting physiology and reproduction. High-quality RNA-Seq data achieved more than 80% mapping efficiency, with a total of 15,782 transcripts were identified. Transcriptome analysis of S. litura larvae showed 323 differentially expressed genes (DEGs), of which 262 genes were significantly upregulated and 61 were downregulated in heat-stressed larvae compared to the control group. The DEGs were associated with stress response, reproduction, signalling, proteostasis, detoxification, oxidative stress, metabolism, development, and chromatin regulation. Heat shock proteins genes, including HSP70, HSP90, and HSP27, together with co-chaperones such as TRET-1, STIP1, and Starvin, were strongly upregulated, indicating enhanced cellular protection against protein damage and oxidative stress under heat stress. Conversely, key reproductive and cell cycle-related genes, including BARR, CAPD2, FEO, CDK2 and MORULA, were significantly downregulated, suggesting reproductive impairment and developmental arrest. RT-qPCR validation corroborated the RNA-Seq findings, demonstrating a heat-induced physiological trade-off that prioritizes survival over reproduction. Consistent with these molecular responses, heat-stressed insects exhibited marked reproductive impairment, including significant reductions in gonadosomatic index, eupyrene sperm bundle count, mating frequency, mating success, female calling behaviour, copulation duration, fecundity, and egg fertility. Collectively, these findings provide comprehensive insights into the molecular basis of thermal adaptation in S. litura and demonstrate that acute heat stress compromises reproductive fitness while activating conserved stress-response pathways that promote short-term survival.

Animals↗
PubMed2026

A contextual activity score (CAS) for inferring ADAR-associated transcriptional activity across RNA-seq, single-cell, and spatial transcriptomics.

BACKGROUND AND OBJECTIVE: Adenosine-to-inosine RNA editing, catalyzed by Adenosine Deaminases Acting on RNA (ADARs), is a widespread modification involved in neural function, immune regulation, and cancer. The Alu Editing Index (AEI) is the standard metric to estimate ADAR activity but requires raw sequencing reads and is poorly suited for single-cell and spatial transcriptomic data. This study aimed to develop an alternative framework for inferring ADAR-associated transcriptional activity from gene expression data across diverse transcriptomic technologies. METHODS: We developed the Contextual Activity Score (CAS), a framework based on transcriptional signatures from ADAR perturbation experiments. Context-specific signatures were generated for human neurons, mouse neurons, and cancer models to infer ADAR1 and ADAR2 activity. CAS was computed from normalized gene expression matrices using regulon-based enrichment analysis. Performance was evaluated by comparing with the Alu Editing Index across bulk RNA sequencing datasets, simulated sequencing depths, and library preparation protocols. RESULTS: CAS showed strong concordance with the Alu Editing Index across multiple datasets, while remaining robust to reduced sequencing depth and different library protocols. Unlike the Alu Editing Index, CAS can be applied to single-cell and spatial transcriptomic data and enables the independent assessment of ADAR2 activity. In cancer and neuronal contexts, CAS captured biologically meaningful variations in ADAR-associated transcriptional activity at sample, cell-type, and spatial levels. CONCLUSION: CAS provides a scalable approach applicable across multiple RNA-seq protocols for estimating ADAR-associated transcriptional activity using gene expression data. This method, implemented in an open-source R package for broad adoption, expands the ability to study ADAR-associated transcriptional activity across transcriptomic modalities where direct editing quantification is challenging, such as single-cell and spatial transcriptomics.

Adenosine Deaminase↗
PubMed2026

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

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

Expression profiles of miRNAs in ruminant intermediate hosts with cystic echinococcosis.

Cystic echinococcosis (CE), caused by the larval stage of Echinococcus granulosus sensu lato (s.l.), is a parasitic zoonotic disease recognized by the World Health Organization as a neglected tropical disease of significant public health concern. Despite ongoing control programs, CE remains endemic, underlining the need for integrated control strategies that involve new diagnostic and therapeutic tools. Recent investigations have spotlighted microRNAs (miRNAs) as key regulators in parasite development, immunomodulation, and as potential diagnostic and therapeutic targets. In the present research, a molecular study was conducted to investigate hydatid cyst samples (protoscoleces and germinal membranes) collected in southern Italy from different ruminant species (sheep, cattle, and water buffaloes), naturally infected with CE, with the ultimate goal of establishing a more comprehensive picture of miRNA expression patterns in these intermediate hosts. The bioinformatic analysis of hydatid cyst samples revealed 168 mature miRNAs. Among these, egr-miR-10-5p, egr-let-7-5p, and egr-miR-71-5p were the most abundant, with egr-miR-10-5p showing particularly high expression levels. No significant differences in miRNA abundance between host species were found. In contrast, when focusing on the comparison between protoscoleces and sterile germinal membranes, 24 miRNAs were found to be differentially expressed. Targeted qPCR of four selected miRNAs (egr-miR-71-5p, egr-let-7-5p, egr-miR-125-5p, and egr-miR-10-5p) showed clear overexpression in protoscoleces and in fertile germinal membranes compared with sterile ones. The differential miRNA expression patterns provide insight into the molecular mechanisms controlling the parasite's lifecycle and may guide the development of novel intervention methods to enhance CE control in endemic areas.

Animals↗
PubMed2026

Biological characterization of Candida parapsilosis haploids induced by voriconazole.

OBJECTIVES: Candida parapsilosis is an important opportunistic fungal pathogen causing serious human infections in nosocomial settings. It has long been thought that C. parapsilosis has a diploid genome with a high homozygosity between chromosome homologs. METHODS: In this study, we report the discovery of C. parapsilosis haploids induced by voriconazole, a triazole with broad antifungal activity against fungal pathogens, in an experimental evolutionary assay. RESULTS: The haploid strains were able to undergo auto-diploidization under in vitro culture conditions or during systemic infection at a low frequency. Compared to the progenitor diploid strain, C. parapsilosis haploid and auto-diploid strains exhibited a reduced ability of invasive growth and biofilm formation. Global transcriptional expression analysis indicated that haploid and auto-diploid strains had a similar transcriptomic profile, which showed a remarkable difference from the progenitor diploid strain perhaps due to the loss of chromosome heterozygosity. Moreover, the haploid and diploid strains had distinct fungal burdens in different animal tissues, suggesting the haploid state could have a colonization advantage over the diploids in certain tissues such as the brain and spleen. CONCLUSIONS: The discovery of C. parapsilosis haploids not only sheds lights on the biology of this important fungal pathogen, but also provides a tool for genetic modifications for the field.

Voriconazole↗
PubMed2026

Ecr positively regulates activity of the PhoQ/PhoP signalling system in Klebsiella pneumoniae.

BACKGROUND: The rising prevalence of polymyxin resistance in multidrug-resistant Klebsiella pneumoniae presents a critical situation with limited therapeutic options. METHODS: Methods Genomic sequencing of 15 clinical polymyxin-resistant K. pneumoniae strains with multidrug resistance revealed that MgrB inactivation, predominantly disrupted by insertion sequences (ISs) in the IS1, IS4, and IS5 families, was the leading cause of polymyxin resistance. Comparative transcriptomics of wild-type, ΔmgrB, and ΔmgrBΔphoP were performed to elucidate the MgrB-PhoPQ regulatory network. RESULTS: This study conducted a system-wide analysis of the regulatory network and identified a species-specific PhoPQ regulon in K. pneumoniae.Beyond the classical MgrB-PhoPQ-ArnBCADTEF pathway, we identified a previously unannotated PhoPQ-regulated gene, 144 bp LN739_RS09850, encoding an Ecr homologue from Enterobacter cloacae. This protein has been reported to confer colistin heteroresistance, with the underlying mechanism not yet functionally validated. This study revealed that overexpression of Ecr homologues decreased colistin susceptibility in both K. pneumoniae and E. cloacae, but this phenotype was abolished upon phoP deletion, confirming PhoP's essential role. Consistent with this dependency, comparative transcriptomics of Ecr-overexpressing K. pneumoniae vs. control revealed significant upregulation of mgrB, phoPQ, arnBCADTE, and pmrD. Two-hybrid bacterial assays further demonstrated direct Ecr-PhoQ interaction. Electrophoretic mobility shift assay confirmed that PhoP directly binds to the ecr promoter in vitro, and a β-galactosidase reporter assay demonstrated that PhoP enhanced ecr promoter activity, indicating that PhoP regulates ecr expression by directly controlling its transcription. CONCLUSION: Collectively, these findings suggest that PhoP may directly activate the transcription of Ecr, with Ecr feedback activating the PhoPQ system via interaction with PhoQ, leading to induction of the arn operon and consequent polymyxin resistance.

Klebsiella pneumoniae↗
PubMed2026

Effects of rumen fluid transplantation on longissimus dorsi muscle development in Xizang sheep: An association analysis based on transcriptomic and serum metabolomic profiles.

This study aimed to investigate the effects of rumen fluid transplantation (RFT) on the growth and development of the longissimus dorsi muscle in female Xizang sheep. After RFT, muscle lightness differed significantly between the two groups, with the LDC group showing significantly higher lightness than the LDT group. In contrast, no significant differences were observed between groups in other muscle phenotypic traits, including drip loss, pH, cooking loss, shear force, redness, and yellowness. Antioxidant-related indices (SOD, GSH-PX, MDA, CAT, and T-AOC) also showed no significant differences between groups. Histological analysis revealed that muscle fiber length, width, and density were significantly greater in the experimental group than in the control group. Transcriptomic analysis identified 515 differentially expressed genes (DEGs), of which 419 were downregulated. KEGG analysis indicated that genes involved in muscle development-related pathways, such as cell adhesion and the PI3K-Akt signaling pathway, were predominantly downregulated. Key serum metabolites (L-kynurenine, IPA, allantoin, and propionylcarnitine) showed highly significant positive correlations with muscle fiber growth indices. In contrast, metabolites such as l-carnitine, acetylcarnitine, and citrulline were negatively correlated with muscle fiber growth, but positively correlated with the expression of muscle structure-related genes (COL11A1 and EFNA5) and with meat lightness. Overall, this study provides new insights into the potential molecular basis by which RFT influences muscle growth and development. However, the mechanisms by which RFT affects muscle development and meat quality-related traits remain unclear and warrant further investigation.

Animals↗
PubMed2026

Transcriptomic insights into the molecular mechanism of antifouling agent-induced settlement inhibition in the Pacific oyster Crassostrea gigas.

Marine biofouling remains a persistent challenge to maritime industries and marine ecosystems worldwide. In this study, we systematically evaluated the acute toxicity, settlement inhibitory efficacy, and underlying molecular mechanisms of an N-oleyl-1,3-propanediamine-based antifouling agent using pediveliger larvae of the Pacific oyster Crassostrea gigas. The 96 h-LC50 of the agent was determined to be 0.81 mg/L, and exposure to 1.68 mg/L achieved complete larval settlement inhibition without inducing significant acute toxicity. Transcriptomic analysis identified 791 differentially expressed genes, dominated by downregulated genes associated with ribosomal function, translation, cell adhesion, and cytoskeletal organization. The agent exerts its inhibitory effect primarily through the global suppression of protein synthesis, disruption of cell-substrate adhesion and cytoskeletal integrity, and induction of proteotoxic stress responses. These findings reveal a multi-pathway molecular mechanism underlying antifouling agent-induced settlement inhibition in oyster larvae and provide key molecular biomarkers to support the development of eco-friendly antifouling technologies.

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

Genome-wide identification and expression analysis of the CREB/ATF family and its potential role in melanogenesis in the Manila clam (Ruditapes philippinarum).

Ruditapes philippinarum is an economically important bivalve species in China, and shell color is a trait of ecological and commercial significance. Melanin is a key determinant of shell color, and members of the CREB/ATF family have been reported to participate in melanogenesis in other organisms. In this study, members of the CREB/ATF family were systematically identified at the whole-genome level based on genomic and transcriptomic datasets, followed by analyses of their phylogenetic relationships, gene structures, and expression patterns. A total of six CREB/ATF family members were identified and classified into five subfamilies. Expression profiling and RT-qPCR validation revealed that most CREB/ATF genes were highly expressed in the mantle and displayed clear differences among shell-color phenotypes. Except for RpATF4 and RpCREBZF, most members exhibited relatively high expression levels in dark-colored shell strains, particularly in black and zebra-striped clams. Moreover, most genes showed low expression during early embryonic and larval stages but increased expression at the single-siphon spat and juvenile stages. These results suggest that the CREB/ATF family may be involved in melanin-associated shell-color regulation in R. philippinarum, providing important candidate genes and a theoretical basis for further elucidating the molecular mechanisms of shell-color formation in mollusks.

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

Muscular fiber properties and multi-omics investigation of larval and adult locomotor muscle in Microhyla fissipes.

During metamorphosis, Microhyla fissipes undergoes a critical transition from an aquatic to a terrestrial lifestyle, accompanied by significant remodeling of skeletal muscle. Notably, larval tail muscle degenerates, while adult hindlimb muscle develops. However, the molecular mechanisms that orchestrate these muscle type-specific adaptations to the changing environment remain unclear. In this study, histological observation, transcriptomics, and metabolomics were integrated to compare locomotor muscles from two stages: larval muscle from tail versus adult muscle from hindlimb. Our results revealed that adult muscle fibers exhibited reduced diameter and shorter sarcomere length compared to those of tadpoles. Transcriptomic analysis identified 4103 differentially expressed genes (DEGs), including 2182 up-regulated and 1921 down-regulated genes. Up-regulated genes were mainly involved in energy metabolism and cellular homeostasis pathways, including PPAR signaling and oxidative phosphorylation, whereas down-regulated genes were associated with carbohydrate metabolism and cell proliferation pathways, such as glycolysis/gluconeogenesis and PI3K-Akt signaling. Metabolic profiling indicated a metabolic shift from anaerobic to aerobic energy production, with 57 differential metabolites identified, mainly involved in protein metabolism and insulin-related pathways. Integrated multi-omics analysis further highlighted the AMPK and FoxO signaling pathways play key roles in this process. In conclusion, our findings demonstrate that the metabolic and structural differences between larval and adult skeletal muscles are mediated by AMPK- and FoxO-dependent signaling pathways, providing novel insights into the molecular mechanisms underlying adaptive development and locomotor transition in anuran amphibians.

Animals↗