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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

PubMed2026

PGR expression as a pharmacogenomic companion biomarker to GENE70-derived genomic risk in ER-positive/HER2-negative breast cancer.

BACKGROUND: The biology of the estrogen receptor-positive (ER+) and human epidermal growth factor receptor 2-negative (HER2-) breast cancers is heterogeneous even when they are categorized by their risk via genomics. Transcriptomic PGR expression reflects endocrine pathway activity and may provide complementary biological information within established GENE70-derived genomic-risk categories. Whether this molecular marker improves the biological interpretation of genomic-risk stratification beyond conventional clinicopathological assessment remains uncertain. OBJECTIVES: The aim of this study was to determine whether transcriptomic PGR expression provides complementary biological and prognostic information within reconstructed GENE70-derived genomic-risk categories and refines the characterization of endocrine-related tumour biology in ER-positive/HER2-negative breast cancer. METHODS: This study analysed publicly available transcriptomic and clinical data from three cohorts: METABRIC (discovery cohort), GSE96058/SCAN-B cohort (validation cohort) and TCGA-BRCA cohort (molecular validation cohort). The GENE70-derived genomic-risk score was reconstructed for each cohort using matched genes. Cox regression, Kaplan-Meier analysis and subgroup comparisons were used to assess relationships between PGR expression, clinicopathologic variables, molecular features and survival outcomes. RESULTS: Across the three independent cohorts, low transcriptomic PGR expression was consistently associated with higher GENE70-derived genomic risk, increased MKI67 expression, reduced ESR1 expression and enrichment of the Luminal B subtype. Survival findings differed between cohorts. In the discovery METABRIC cohort, transcriptomic PGR expression showed heterogeneous associations with survival, particularly within GENE70-derived high-risk subgroups, whereas the external GSE96058/SCAN-B validation cohort demonstrated consistent associations between low PGR expression and poorer overall survival in both the overall ER-positive/HER2-negative population and GENE70-derived high-risk subgroups. CONCLUSION: These findings suggest that transcriptomic PGR provides complementary biological and prognostic information within GENE70-derived genomic-risk categories. However, because treatment response was not evaluated in the present study, the findings should not be interpreted as evidence of predictive or pharmacogenomic utility and prospective studies incorporating treatment-response analyses are required before such applications can be established.

Humans↗

PubMed2026

Transcriptomic and Metabolomic Profiling Identifies a Core Gene-Metabolite Axis Driving African Swine Fever Virus Replication in the Soft Tick Ornithodoros lahorensis.

African swine fever virus (ASFV) causes an incurable swine disease with nearly 100% mortality, posing a catastrophic threat to global pig production. The soft tick Ornithodoros lahorensis acts as a critical biological vector that sustains persistent ASFV replication and mediates long-distance viral transmission, yet the molecular mechanisms governing ASFV-tick interplay remain poorly understood. Here, we integrated transcriptomics and metabolomics to systematically dissect molecular changes in O. lahorensis across three infection stages: Uninfected control, early infection (7 days post-infection, dpi), and late persistent infection (21 dpi). Multi-omics integration revealed that ASFV extensively remodels tick host metabolism, predominantly activating purine/pyrimidine metabolism, lipid biosynthesis, and energy metabolism. We further characterized a conserved regulatory module consisting of 12 core genes and 8 signature metabolites that collectively support ASFV genome replication and virion assembly. Three hub metabolic genes (TK1, ATP5F1B, and IMPDH) were selected for functional validation via siRNA silencing in ticks; individual gene silencing suppressed ASFV loads by 89.2%, 91.5%, and 87.8%, respectively (p < 0.001***). This work represents the first comprehensive multi-omics investigation of ASFV infection in O. lahorensis. We identified tick-specific molecular targets to block vector-mediated ASFV spread and established a standardized multi-omics analytical pipeline for tick-virus interaction research. Our findings elucidate the mechanistic basis of long-term ASFV persistence in soft ticks and deliver novel actionable clues for developing vector-targeted ASF intervention strategies.

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

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

Transcriptomic and RNAi analyses reveal chloride channel 3-associated osmoregulation in Litopenaeus vannamei under low-salinity stress.

Chloride channels and transporters are important for cellular volume regulation and salinity adaptation in euryhaline crustaceans, yet the intestinal transcriptional relationship between plasma-membrane and intracellular chloride pathways remains unclear in Litopenaeus vannamei. In this study, RNA interference of anoctamin 1 (ANO1) was combined with intestinal transcriptome sequencing under the production-relevant low-salinity condition of salinity 3. ANO1 silencing produced a focused transcriptional response, with 16 differentially expressed genes (DEGs) identified (11 upregulated and 5 downregulated). Functional enrichment indicated that these genes were associated with transporter activity, cytoskeletal organization, extracellular matrix-receptor interaction, membrane lipid metabolism, and vesicular processes. Notably, a transcript encoding chloride channel protein 3 (CLC-3) was significantly upregulated following ANO1 knockdown, suggesting a potential transcriptional relationship between ANO1 and CLC-3 in chloride homeostasis. Based on this finding, CLC-3 was selected for full-length cDNA cloning, sequence characterization, salinity-gradient expression analysis, and RNAi-based functional assessment. The cloned CLC-3 cDNA was 2883 bp in length and encoded an 850 amino acid protein containing a conserved voltage-gated chloride channel (Voltage-CLC) domain and two cystathionine β-synthase domains. Phylogenetic analysis placed LvCLC-3 within the intracellular CLC-c clade, and tissue distribution analysis showed the highest CLC-3 expression in the intestine. Intestinal CLC-3 expression responded nonlinearly to salinity variation, peaking at salinity 20. Under salinity 3, CLC-3 knockdown reduced ANO1, Na+/K+-ATPase alpha subunit, and Na+-K+-2Cl- cotransporter transcript levels, whereas glutamate-gated chloride channel expression increased. Mild hepatopancreatic structural alterations were also observed after CLC-3 knockdown. These findings suggest that CLC-3 is a salinity-responsive intracellular chloride-transporter candidate associated with intestinal ion-transport-related transcriptional responses after ANO1 suppression in L. vannamei, although the underlying physiological mechanism requires further validation.

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

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

Integrated transcriptomic and metabolomic analysis of fluoride tolerance-related pathways and differentially expressed genes in silkworm strain XSKD.

XueSong KD (XSKD) silkworm strain exhibits prominent fluoride tolerance, yet the underlying molecular mechanisms of fluoride tolerance remains unclear. In the present study, fourth-instar pre-molting XSKD silkworms were used as experimental materials for integrated transcriptomic and untargeted metabolomic analyses. In total, 572 differentially expressed genes and 90 differential metabolites were screened. GO enrichment and KEGG enrichment based on the hypergeometric distribution model revealed that 13-Hydroxy-9Z,11E-octadecadienoic acid (13-(S)-HODE) acts as the core differential metabolite, which is significantly enriched in the linoleic acid metabolism pathway. Within this pathway, LOC101737302 and CYP338A1 display opposite expression trends and show correlations with pathway metabolites. Based on multi-omics data, this study preliminarily characterizes the lipid metabolic response under fluoride stress, providing omics dataset support for further in-depth exploration of the molecular mechanism of fluoride tolerance in silkworms.

Animals↗

PubMed2026

Integrated transcriptomic and metabolomic analyses provide new insights into the response of black rockfish (Sebastes schlegelii) larvae to temperature fluctuations.

Sebastes schlegelii usually encounter elevated and fluctuating water temperatures near its upper thermal limit in summer, yet the hepatic responses of larvae to repeated temperature fluctuation regimes remain unclear. To address this question, S. schlegelii larvae were exposed for 8 days to four thermal regimes: constant 18 °C (CT), constant 28 °C (HT), intermittent cooling from 18 to 8 °C followed by recovery to 18 °C (FL), and intermittent warming from 18 to 28 °C followed by recovery to 18 °C (FH). Survival rate was evaluated, and integrated liver transcriptomic and metabolomic analyses were performed. Final survival rates were 96.67% in the CT group, 97.78% in the FL group, and 77.78% in the FH group. Survival rate in the HT group (38.89%) was significantly lower than that in the other three groups (P < 0.05). HTvsCT, FLvsCT, FHvsCT, and FHvsHT comparisons identified 2598, 1207, 622, and 2404 differentially expressed genes and 627, 606, 690, and 610 differential metabolites, respectively. KEGG enrichment analyses of DEGs and SDMs in HTvsCT highlighted HSP-mediated proteostasis, endoplasmic-reticulum protein processing, branched-chain and sulfur amino acid metabolism, glutathione metabolism, and central carbon metabolism, with upregulated hsp90aa1, bckdha, gclc, and pfkp and reduced levels of branched-chain amino acids and methionine. Compared with HT, FH showed attenuated disturbances in proteostasis, amino acid and redox regulation, and central carbon metabolism, together with recovery-associated glycerophospholipid turnover. FL primarily induced polyunsaturated fatty acid (PUFA)-related membrane lipid remodeling. These findings indicate that hepatic responses differed between continuous high-temperature exposure and temperature fluctuations and between fluctuation regimes.

Animals↗

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

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

Distinct cell morphotypes of Aureobasidium melanogenum ZN exhibit differential functional profiles in promoting maize growth.

Black yeast-like fungi of the genus Aureobasidium exhibit morphological plasticity, but whether distinct cellular states within the same genetic background are associated with different plant growth-promoting functions remains unclear. Here, yeast-like cells (YL), swollen cells (SC), and chlamydospores (CH) of Aureobasidium melanogenum ZN were characterized. YL was associated mainly with siderophore production and laccase activity, SC with extracellular polysaccharide accumulation, and CH with phosphate mobilization and higher ammonia and IAA production. Whole-genome and comparative genomic analyses revealed a shared repertoire related to nutrient acquisition, auxin-associated metabolism, extracellular oxidation, and carbohydrate remodeling, with expansions in nutrient- and cell-surface-related gene families. Transcriptomic and metabolomic analyses showed distinct deployment of these capacities, with CH exhibiting broad reprogramming of tryptophan-associated, nitrogen, phosphate, central-carbon, and amino-acid metabolism. In maize, CH at the optimal inoculation concentration of 105 CFU·mL-1 produced the strongest growth promotion, increasing plant height, dry biomass, root length, root surface area, and root volume by 58.6%, 365.1%, 191.0%, 194.3%, and 222.4%, respectively. Consistent with this pronounced growth phenotype, maize root transcriptomics showed coordinated CH-induced responses involving root development, nutrient transport, redox regulation, and root-interface remodeling. Root-zone tracking showed greater short-term stability and persistence of CH. These findings identify cellular state as an important functional dimension of Aureobasidium-plant interactions and provide a basis for developing fungal inoculants with defined beneficial cellular states.

Zea mays↗

PubMed2026

Integrated physiological and transcriptomic analyses reveal coordinated gill responses to heat stress in pikeperch (Sander lucioperca).

Climate change-driven warming of aquatic environments has made thermal stress an increasingly important factor influencing fish physiological homeostasis. Given their central roles in respiration and osmoregulation, gills are particularly responsive to variations in ambient temperature. Histological examination, physiological measurements, and transcriptome profiling were integrated to investigate the mechanisms associated with heat stress-induced gill injury in pikeperch (Sander lucioperca). Histological analysis revealed that exposure to 29 °C directly caused structural damage to the gills of pikeperch. Oxidative status was evaluated by measuring malondialdehyde (MDA) levels and the activities of antioxidant enzymes, including superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT). MDA accumulation was significantly enhanced under heat stress, while antioxidant enzyme activities (SOD, POD, and CAT) displayed a transient increase followed by a subsequent decline. Transcriptome profiling showed marked enrichment of the protein processing in endoplasmic reticulum pathway after heat stress, suggesting activation of endoplasmic reticulum (ER) stress in pikeperch gills. With increasing stress duration, the unfolded protein response (UPR) appeared unable to re-establish ER homeostasis, shifting ire1 and atf6 toward a pro-apoptotic state. Protein-protein interaction (PPI) analysis further highlighted hub genes potentially involved in heat stress-induced ER stress and apoptosis. TUNEL staining and western blotting collectively confirmed that heat stress triggered apoptosis in pikeperch gill tissue. Overall, this study provides new insights into the physiological and molecular responses of pikeperch gills to heat stress and enhances our understanding of thermal stress adaptation in cold-water aquaculture species under climate change.

Animals↗

PubMed2026

Distinct molecular responses to acute cold exposure revealed by comparative transcriptomic and metabolomic profiling in the bay scallop Argopecten irradians.

Acute cold stress can elicit distinct molecular responses even when bay scallop populations show similar phenotypic outcomes. We compared a seventh-generation fast-growing bay scallop line (BS) with a commercial control population (CC) during a 72-h acute cold exposure at -1 ± 0.3 °C. RNA-seq was used as the discovery layer, representative BS cold-responsive genes were evaluated by qRT-PCR, and paired LC-MS profiles provided a comparative metabolic layer. At baseline, 138 genes differed between BS and CC; after cold exposure, 134 of these baseline differences disappeared and 61 of 65 cold-state differences newly emerged. BS showed a larger transcriptomic response magnitude than CC, with 1129 cold-responsive genes compared with 28 genes in CC, and this ordering remained robust across multiple sensitivity analyses. Survival after 72 h was identical in BS and CC (83/90, 92.2% in each population). Biochemical responses were time-dependent and marker-specific: CAT, LZM, T-SOD and T-AOC showed population-by-time interactions, whereas GSH-Px and MDA did not, and the 72-h differences were not consistently favourable to BS. Metabolomic cold effects were strongly concordant between populations, and no feature showed a significant population-by-cold interaction. Features putatively assigned to arachidonic acid metabolism were enriched, but this provider-annotated pathway signal remains exploratory because authentic-standard confirmation was not performed. These findings indicate population-specific differences in molecular responsiveness but do not establish superior cold tolerance in BS.

Animals↗

PubMed2026

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

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

Animals↗

PubMed2026

Antennal transcriptome analysis of chemosensory proteins in the raspberry weevil, Aegorhinus superciliosus (Coleoptera: Curculionidae).

Aegorhinus superciliosus (Coleoptera: Curculionidae) is a polyphagous pest of economic importance in southern Chile, the chemical ecology of which remains poorly characterized. Across insect species, chemosensory proteins, including odorant receptors (ORs), gustatory receptors (GRs), ionotropic receptors (IRs), odorant-binding proteins (OBPs), chemosensory proteins (CSPs), and sensory neuron membrane proteins (SNMPs), mediate the detection of chemical cues involved in host selection, reproduction, and other ecologically relevant behaviors. In this study, the antennal transcriptome of adult A. superciliosus was sequenced and analyzed using a de novo RNA-seq approach. Three independent biological replicates per sex were used for RNA-seq, and the same number of independent biological replicates was used for RT-qPCR validation; sequencing yielded 147,409,936 high-quality reads after quality filtering. A total of 112 candidate chemosensory genes were identified, comprising 43 ORs, 34 OBPs, 10 CSPs, 18 IRs, 5 GRs, and 2 SNMPs. Phylogenetic analyses assigned these candidate proteins to established clades, providing a comparative framework for functional inference for ORs and OBPs. Sex- and tissue-biased expression analyses revealed that several ORs, including AsupOR4, AsupOR19, and AsupOBP13, exhibit antennal enrichment and sex-specific expression patterns. Notably, AsupOR19 and AsupOBP13 displayed strong female-biased expression. In addition, transcripts of selected ORs and OBPs were detected in non-antennal tissues, such as the rostrum and legs, suggesting potential functional versatility beyond canonical olfaction. Together, these findings represent the first molecular identification of the chemosensory repertoire of A. superciliosus. This study establishes a foundation for reverse chemical ecology approaches aimed at identifying behaviorally active volatile organic compounds (VOCs) toward environmentally sustainable strategies for integrated pest management.

Animals↗

PubMed2026

Transcriptomics reveals species-specific adaptive strategies to calorie restriction in two Argopecten scallops with distinct lifespans.

Calorie restriction (CR) is a well-established non-genetic intervention for lifespan extension in multiple model organisms. Seasonal food shortage in cold and temperate seas may mimic CR, inducing in bivalves a response similar to that in vertebrates and thereby prolonging life expectancy. However, the relationship and the mechanism underlying the food availability and lifespan in bivalves remain largely unexplored. Two closely related scallop species the short-lived warm-water Argopecten irradians (lifespan <2 years) and the longer-lived cold-water Argopecten purpuratus (7-10 years) provide an ideal comparative system to investigate species-specific adaptive strategies. In this study, we subjected both species to CR for 30 and 56 days and performed comparative transcriptomic profiling, weighted gene co-expression network analysis (WGCNA), and physiological assays to elucidate their distinct molecular responses. Transcriptomic analysis revealed that A. purpuratus exhibited substantially more DEGs than A. irradians at both time points under CR, with both species showing downregulation of metabolic pathways but to different extents. A. irradians mounted an early nutrient-sensing response at 30 days (IGF1R, PIK3R3, INSR suppression), indicating acute sensitivity to limitation; by contrast, A. purpuratus displayed delayed FoxO activation at 56 days, along with its downstream effectors NFKBIA, CREB3L4, and SMAD4, suggesting a gradual adaptive program may link to its extended lifespan. WGCNA identified three negatively correlated modules in each species, with coral2 being the most prominent in A. irradians and darkolivegreen in A. purpuratus. The former was dominated by ciliary motility genes, whereas the latter featured coordinated repression of oxidative phosphorylation. Additionally, both species exhibited conserved suppression of mTOR/S6K growth signaling and activation of cellular maintenance programs. Collectively, these findings expand the understanding of CR-mediated longevity regulation in bivalves and provide candidate gene resources for future functional studies and breeding programs.

Pectinidae↗

PubMed2026

Sex- and development-specific transcriptomic profiling of venom and silk genes in the wolf spider Pardosa astrigera provides insights into ecological adaptation and predatory strategies.

Spider venom and silk glands are two major secretory systems that contribute to prey capture, defense, and reproduction, but their sex- and development-specific molecular regulation in wandering wolf spiders remains poorly understood. Here, the transcriptome of Pardosa astrigera, an important agricultural natural enemy in China, revealed significant sex- and development-associated molecular differentiation among adult females, adult males, and spiderlings. A total of 100,025 unigenes were obtained, of which 23,852 were functionally annotated, providing a comprehensive transcriptomic resource for this species. Differential expression patterns showed marked variation among groups, with 531, 1792, and 832 DEGs detected in PAF vs PAS, PAM vs PAS, and PAF vs PAM, respectively. These genes were mainly associated with metabolic, oxidation-reduction, cuticle development, MAPK signaling, and lysosome pathways. Fifteen co-expression modules revealed distinct expression patterns. The turquoise, pink, yellow, and red modules were development-related, whereas the blue module was male-biased. Venom- and spidroin-related genes were distributed across multiple modules, suggesting coordinated regulation. Overall, 42 venom peptides, 21 venom proteins, and 11 spidroins were identified. Representative genes showed strongly biased expression, including spiderling-biased U3_Pp1a and U5_Pp1e, female-biased U4_Pp1a, and male-biased SMase D_108750 and PaTuSp_108466. These findings reveal sex- and development-biased expression patterns of venom- and silk-related candidate genes in P. astrigera and may provide molecular insights into ecological adaptation and predatory strategies in wandering wolf spiders.

Animals↗