Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Target 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.

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

898 records · Page 7Linked to original sources
PubMed2026

Simultaneously PYCR-1 and ALH-6 inhibition exacerbates 6-PPD quinone toxicity via disrupting proline and glutamate metabolisms and activating insulin signals in Caenorhabditis elegans.

Glutamate synthesized from the proline can serve as a precursor for key intermediate metabolites of citric acid cycle. Recently, we observed reduced glutamate content and expression of alh-6 controlling glutamate synthesis by 6-PPD quinone (6-PPDQ) in Caenorhabditis elegans. However, possible effect of 6-PPDQ on proline synthesis and the association with 6-PPDQ toxicity induction remain unclear. After 0.1-10 μg/L 6-PPDQ exposure, proline content was further reduced, and expression of pycr-1 governing proline biosynthesis was decreased. In 6-PPDQ exposed nematodes, RNA interference (RNAi) of pycr-1 decreased α-ketoglutarate content, enhanced mitochondrial dysfunction, reduced nicotinamide adenine dinucleotide (NADH) and reduced flavine adenine dinucleotide (FADH₂) contents, inhibited mitochondrial complex I/II activities, and decreased expressions of gas-1 and mev-1. Moreover, compared to single RNAi, double RNAi of pycr-1 and alh-6 exacerbated the 6-PPDQ toxicity in reducing α-ketoglutarate, NADH, and FADH₂ contents, and suppressing mitochondrial complex I/II activities and gas-1 and mev-1 expressions. Additionally, double RNAi of pycr-1 and alh-6 intensified toxicity of 6-PPDQ on longevity and caused upregulation of insulin ligand and receptor genes and downregulation of daf-16 and its targeted genes in 6-PPDQ exposed nematodes. Furthermore, after 6-PPDQ exposure, daf-16 RNAi suppressed pycr-1 and alh-6 expressions, suggesting formation of a regulatory feedback loop between pycr-1/alh-6 and daf-16. Our findings highlight involvement of disrupted proline and glutamate metabolisms in 6-PPDQ-induced mitochondrial dysfunction and reduced longevity.

Animals↗
PubMed2026

Molecular diagnostic yield and barriers in inherited retinal diseases: a retrospective cohort study.

OBJECTIVE: To evaluate the diagnostic yield of panel-based genetic testing for inherited retinal diseases (IRDs) and identify barriers to molecular resolution. DESIGN: Retrospective cohort. PARTICIPANTS: A total of 404 patients with clinically confirmed IRDs who were evaluated at the Adult Inherited Retinal Dystrophy Service, Ontario, Canada (October 2021-September 2024). METHODS: Patients underwent targeted massive parallel sequencing panel testing. Diagnostic yield was calculated, and unresolved cases were reviewed. Associations between yield, phenotype, ethnicity, and sex were assessed using χ² analysis. RESULTS: Of 685 referrals, 570 had confirmed IRDs. After we excluded 140 pending results and 26 patients who declined testing, 404 patients were analyzed. At referral, 94 patients (23.2%) had a previous molecular diagnosis, and 138 (34.0%) were diagnosed through clinic-initiated testing, giving an overall yield of 57.4%. Yield varied significantly by phenotype (χ², P = 1.4 × 10⁻⁶), from 94.4% in vitelliform macular dystrophies to 25.0% in vitreoretinopathies, with no sex association (P = 1.0). Disease-causing variants were identified in 83 IRD-associated genes, most frequently ABCA4, USH2A, and BEST1. Of 172 unresolved cases, 62 (36.0%) had negative panels, and 110 (63.9%) were inconclusive, including 30 with unphased pathogenic variants in recessive genes and 10 with high-suspicion variants of uncertain significance. Key barriers included limited family availability for phasing, restricted access to functional assays, and lack of public coverage for whole-exome or whole-genome sequencing. CONCLUSIONS: Massive parallel sequencing-based panel testing achieved a 57% diagnostic yield in this IRD population. Success was strongly phenotype-dependent with substantial heterogeneity. Whole-exome sequencing, whole-genome sequencing, family segregation, and functional genomics could improve diagnostic outcomes and management.

Humans↗
PubMed2026

TNFα-dependent modulation of WT1-MMP9 regulatory axis links developmental and inflammatory pathways in glaucoma.

Glaucomas are heterogeneous optic neuropathies associated with extracellular matrix dysregulation, abnormal ocular morphogenesis, and inflammatory signaling. Targeted deep sequencing of 586 primary congenital glaucoma (PCG) cases and 1,757 controls identified rare pathogenic variants in multiple genes, including WT1 and MMP9. Notably, WT1 variants clustered within the nuclear export sequence. Further, functional analyses showed that combined wt1-pax6 suppression in zebrafish disrupted ocular morphogenesis, highlighting developmental interdependence. In human trabecular meshwork cells, WT1 acted as a transcriptional repressor of MMP9, while TNF-α signaling triggered nitric oxide-dependent nuclear export of WT1, resulting in delayed MMP9 upregulation. This effect was reversible by inhibiting nuclear export or nitric oxide synthase. A patient-derived mutation in the nuclear-export region of WT1, disrupted this regulatory switch, causing abnormal MMP9 expression. These findings position WT1 as an important regulator linking developmental and inflammatory mechanisms in glaucoma pathogenesis.

anterior segment dysgenesis↗
PubMed2026

Mechanisms of Hematopoietic Stem Cell Aging and Emerging Rejuvenation Strategies.

Hematopoietic stem cell (HSCs) aging is a complex biological process driven by both cell-intrinsic alterations and extrinsic cues from the bone marrow niche. Understanding these mechanisms is critical for developing therapies against aging-related hematopoietic disorders. This review synthesizes recent advances in the molecular mechanisms underlying HSCs aging, including microenvironmental aging, genomic instability, epigenetic dysregulation, mitochondrial dysfunction, and aberrant nuclear mechanotransduction. We summarize that the functional decline of HSCs during aging drives a compensatory expansion of the phenotypically defined stem cell pool, leading to an aberrant increase in cell number. We also highlight aging-associated HSCs heterogeneity, including CD150high and P-selectin-positive subsets that enrich for myeloid-biased or functionally compromised HSCs states while emphasizing that surface phenotype alone may not fully indicate functional rejuvenation. Finally, we discuss emerging rejuvenation strategies-including targeting myeloid-biased HSCs, modulating inflammatory pathways, and implementing epigenetic or metabolic interventions-supported by cutting-edge technologies such as single-cell multi-omics, gene editing, and computational modeling. These approaches hold promise for counteracting age-related hematopoietic decline and restoring immune competence.

Humans↗
PubMed2026

Flavoromics-based profiling reveals taste and aroma differences between infant formula and breast milk.

Flavor differences between infant formula (IF) and breast milk (BM) are considered a potential factor affecting infants' acceptance of IF. This experiment employs flavoromics combined with multivariate statistical analysis to systematically compare the flavor profiles of IF and BM. Electronic tongue analysis and amino acid correlation revealed that IF was characterised by pronounced saltiness and umami richness, whereas BM exhibited greater bitterness and astringency. Volatile compound profiling identified five key flavor constituents in IF, predominantly aldehydes such as hexanal and pentanal. In contrast, BM contained a broader array of compounds-including acids, aldehydes, and esters-resulting in a more complex flavor profile. Kyoto Encyclopedia of Genes and Genomes (KEGG)-based metabolic pathway annotation, together with fatty acid profiling, suggested that some volatiles may be associated with lipid oxidation, Maillard reaction and sulfur-containing amino acid degradation pathways, offering a theoretical basis for the targeted optimisation of IF flavor.

Humans↗
PubMed2026

Importin-8 silencing reduces Ataxin-3 aggregation and alleviates Machado-Joseph disease/spinocerebellar ataxia type 3.

Machado-Joseph disease (MJD) is a polyglutamine neurodegenerative disorder characterized by the formation of neuronal intranuclear Ataxin-3 inclusions. The nucleocytoplasmic transport is known to be profoundly dysregulated from early stages of neurodegeneration, further aggravating its progressive nature. To understand the contribution of these proteins in MJD, we screened the effect of 34 transport proteins on Ataxin-3 aggregation. Importin-8 (Ipo8 gene) was selected for further neuropathology and behavior analysis to evaluate its modulatory effect in vitro and in vivo. Our results showed that Ipo8 mRNA expression is increased in MJD and that, upon Ipo8 downregulation, Ataxin-3 aggregation and neuronal dysfunction are reduced in vitro and in vivo. Furthermore, Ipo8 silencing was demonstrated to alleviate ataxic traits in vivo. Collectively, our findings indicate that Importin-8 may exert its protective effect by modulating NF-kb/p65 and Argonaute-2 proteins, also previously linked to MJD, standing as a promising therapeutic target to delay disease progression.

Argonaute-2↗
PubMed2026

In silico identification of DNMT1 inhibitors from the PlantCyc database through computational approach to assess the anti-cancer potential of nutraceutical compounds in breast cancer.

Breast cancer accounts for a disproportionate share of global cancer-related deaths, with 670,000 fatalities and 2.3 million new diagnoses recorded in women during 2022 alone. Existing treatment modalities carry considerable toxicity burdens, and resistance to available agents remains an unresolved clinical problem. DNA methyltransferase 1 (DNMT1), the enzyme chiefly responsible for maintaining genome-wide methylation patterns during DNA replication, has been mapped out as a high-value target in breast cancer because its dysregulation silences tumour suppressor genes through promoter hypermethylation. The present work involves hierarchical in silico workflow to screen 4549 plant-derived compounds from the PlantCyc database (v16.0.3) against the human DNMT1 catalytic domain (PDB ID: 4WXX). Ten top-scoring compounds were taken forward for molecular docking via AutoDock Vina; Quercetin and Kaempferol both recorded the highest binding affinities at -9.5 kcal/mol, Wogonin (-9.3 kcal/mol) and Xanthohumol (-8.1 kcal/mol) also emerged as strong binders. Pharmacokinetic evaluation using ADMET-AI confirmed that all 10 compounds met Lipinski's rule of five, with human intestinal absorption values at or above 0.98. Wogonin and Xanthohumol were selected for a 100 ns all-atom molecular dynamics (MD) simulation in GROMACS due to their well-rounded ADMET profiles and limited existing data on their specific interactions with DNMT1 in breast cancer. Across all measured trajectory metrics, backbone RMSD, residue fluctuation, radius of gyration, solvent-accessible surface area, and intermolecular hydrogen bond count, Wogonin formed a more stable, compact complex. These findings suggest that Wogonin and Xanthohumol are non-toxic nutraceutical candidates suitable for DNMT1 targeted epigenetic therapy, with computational foundation strong enough to facilitate future in vitro and in vivo validation work.

Humans↗
PubMed2026

In vivo genome-wide CRISPR screens identify FOXR1 as a suppressor of CD8+ T cell antitumor immunity.

T cell dysfunction critically limits the efficacy of T cell-based immunotherapies in solid tumors, yet the intrinsic regulators of T cell dysfunction remain incompletely understood. Through an in vivo genome-wide CRISPR screen in tumor-infiltrating CD8+ T cells, we identified Forkhead Box R1 (FOXR1) as a potent transcriptional suppressor of CD8+ T cell effector functions. Genetic ablation of FOXR1 significantly enhanced cytokine production and cytotoxic capacity in both murine and human CD8+ T cells, whereas its overexpression impaired T cell activation and effector molecule expression. Mechanistically, multiomics integration of RNA-seq, CUT&Tag-seq, and ATAC-seq revealed that FOXR1 binds directly to promoter regions of key effector genes, including IL2, GZMB, and PRF1, and represses their expression. Importantly, FOXR1 deletion in human anti-CD19 CAR T cells improved their efficacy against solid tumors, demonstrating that FOXR1 is a checkpoint of T cell effector function and targeting FOXR1 is a promising strategy to enhance CAR T cell efficacy against solid tumors.

Animals↗
PubMed2026

Genomic insights into end-use grain quality and nutritional traits of an ancient Indian dwarf wheat ( Triticum sphaerococcum Percival) population using a multi-locus genome-wide association study.

BACKGROUND: Triticum sphaerococcum, an ancient hexaploid wheat species, is renowned for its stress resilience and superior nutritional quality. A panel of 116 T. sphaerococcum accessions (the largest known collection at a single site globally), with six bread wheat released varieties, was evaluated for its potential for genetic quality improvement. Field experiments were conducted under standard, heat and moisture-deficit conditions across two cropping seasons for ten grain end-use quality and nutritional traits. RESULTS: Genotypes showed highly significant differences (P ≤ 0.001) for measured traits, with high broad-sense heritability resulting from substantial genotypic variance contributions. Triticum sphaerococcum consistently outperformed T. aestivum across environments, with moisture-deficit stress proving more detrimental to quality parameters than heat stress, while micronutrient content increased under stressed conditions. Trait correlations revealed that the gluten index (GI) correlated negatively with the grain hardness index (GHI), wet gluten (WG), and water-binding capacity (WB), while positively correlating with dry gluten (DG) and protein content (PRO), whereas grain iron (GFE), zinc (GZN), and protein showed consistent positive interrelationships. Two superior accessions, PAUTS10 (WG 35.13%, DG 13.71%, PRO 16.42%, GZN 50.89 ppm) and Sonamoti (WG 33.33%, DG 12.92%, PRO 16.27%, GZN 56.03 ppm), were identified, surpassing the best check variety HD3226 for quality and nutritional parameters. Multi-locus genome-wide association studies identified 30 stable quantitative trait nucleotides across environments, with candidate gene analysis revealing genes involved in transcription regulation, biosynthetic processes, metal ion homeostasis, and transport. CONCLUSIONS: Triticum sphaerococcum demonstrated superior grain quality and micronutrient potential compared with modern wheat, highlighting its value as a genetic resource for biofortification. The identification of elite accessions and stable quantitative trait nucleotides (QTNs) provides useful targets for breeding programs aimed at improving protein and micronutrient content. Integrating ancient germplasm with modern genomic tools can accelerate the development of nutritionally enhanced wheat varieties. © 2026 Society of Chemical Industry.

Triticum↗
PubMed2026

Modelling the effects of biological intervention in a dynamical gene network.

Cellular response to environmental and internal signals can be modeled by dynamical gene regulatory networks (GRN). In the literature, three main classes of gene network models can be distinguished: (1) non-quantitative (or data-based) models which do not describe the probability distribution of gene expressions; (2) quantitative models which fully describe the probability distribution of all genes co-expression; and (3) mechanistic models which allow for a causal interpretation of gene interactions. We propose two rigorous frameworks to model gene alteration in a dynamical GRN, depending on whether the network model is quantitative or mechanistic. We explain how these models can be used for design of experiment, or, if additional alteration data are available, for validation purposes or to improve the parameter estimation of the original model. We apply these methods to the Gaussian graphical model, which is quantitative but non-mechanistic, and to mechanistic models of Bayesian networks and penalized linear regression.

Gene Regulatory Networks↗
PubMed2026

Genome-wide identification and functional validation of asparagine synthetase genes (NtASNs) in Nicotiana tabacum.

Asparagine (Asn) is pivotal for plant nitrogen (N) metabolism and plays indispensable roles in plant growth, development, and stress tolerance. However, the systematic characteristics and core functions of asparagine synthetase genes (NtASNs) in tobacco remain unclear. Through a comprehensive genome-wide investigation, nine members of the NtASN gene family were identified. Subsequent CRISPR/Cas9-mediated knockout and overexpression assays of these NtASN genes revealed that NtASN1e, NtASN2a, and NtASN2b are the core genes responsible for Asn biosynthesis in tobacco. Their knockout reduced asparagine synthetase activity and Asn content, delayed seed germination by 2-3 days, and displayed elevated oxidative injury when exposed to salinity conditions. In contrast, overexpression of these genes elevated Asn accumulation. Subcellular localization analysis indicated that NtASN1e was localized to both the cytoplasm and chloroplasts, whereas NtASN2a exhibited dual localization in the cytoplasm and endoplasmic reticulum, and NtASN2b was mainly localized in the cytoplasm. This study systematically clarifies the evolutionary characteristics and core functions of the NtASN gene family and provides candidate genes for optimizing nitrogen metabolism and improving salt-stress adaptation in tobacco. These findings hold important practical significance for molecular breeding and product quality improvement in industrial crops.

Nicotiana↗
PubMed2026

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

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

Humans↗
PubMed2026

Algae-to-host horizontal gene transfer in Paramecium bursaria is associated with host adaptation during endosymbiosis.

Paramecium bursaria maintains a stable endosymbiosis with green algae, yet the evolutionary consequences of this association remain unclear. Here, we screened the host genome for algal-derived horizontally transferred genes (HTGs) using a lineage-aware workflow designed to detect horizontal gene transfer (HGT) between two defined lineages. We identified 16 candidate HTGs, including four putative newly transferred genes and 12 homologous transferred genes, most of which were functionally associated with redox homeostasis and metabolism. Five HTGs showed symbiosis-dependent expression. RNAi knockdown of GH32s and SATs reduced host proliferation, total cell area, and motility, while GH32s knockdown also reduced endosymbiont load. Duplication patterns suggest that most transfers may have occurred after the P. bursaria lineage diverged from the sampled Paramecium species but before its lineage-specific whole-genome duplication (WGD). The HTGs also showed host-associated shifts in GC content and gene length, while representative HTGs retained conserved domains and functional motifs. Together, our results support algae-to-host HGT in P. bursaria and suggest that some transferred genes may contribute to metabolic integration during endosymbiosis.

Gene Transfer, Horizontal↗
PubMed2026

Comprehensive identification and evolutionary analysis of the Wnt gene family in bivalves: Insights into the larval development of the noble scallop Chlamys nobilis.

The Wnt gene family regulates fundamental developmental processes in metazoans, but its evolutionary composition and developmental deployment in bivalves remain largely unresolved. Here, we performed a comparative genomic analysis of Wnt genes in 19 bivalve species and examined developmental expression profiles in the noble scallop Chlamys nobilis, with Crassostrea gigas and Chlamys farreri used for cross-species comparison. A total of 235 Wnt genes were identified and assigned to 12 subfamilies. No reliable Wnt3 ortholog was detected in any analyzed bivalve, supporting the view that Wnt3 loss occurred early during lophotrochozoan evolution rather than representing a lineage-specific absence. Most Wnt proteins retained the conserved WNT domain, indicating strong structural conservation, whereas lineage-specific copy-number variation and gene loss were observed among species. C. farreri and C. gigas each retained 12 Wnt genes and lacked Wnt3, whereas C. nobilis lacked Wnt3, Wnt7, and Wnt16. Developmental transcriptome analysis and RT-qPCR revealed clear stage-specific expression patterns. In C. gigas, Wnt2/10/A were highly expressed during earlydevelopment and peaked around the D-shaped larval stage, while Wnt8 and Wnt11 showed distinct stage-specific peaks. By contrast, Wnt1/5/6/9 were more active during later larval development or juvenile formation. These results provide a comparative framework for bivalve Wnt evolution and identify candidate Wnt genes potentially involved in larval development and aquaculture-relevant developmental transitions.

Animals↗
PubMed2026

Genomic and Phenotypic Characterization of Two Novel Enterobacter Phages With EDTA-Enhanced Antibiofilm Activity.

Multidrug-resistant members of the Enterobacter cloacae complex (ECC) are increasingly linked to difficult-to-treat infections and biofilm-mediated antimicrobial tolerance. Here, two lytic phages, vB_EhoIP_HHH and vB_EluM_RZH, displaying podovirus-like and myovirus-like morphology, respectively, were isolated from the River Chelt. HHH has a 39,582 bp genome (51.2% GC, 63 ORFs), while RZH has a 174,197 bp genome (39.4% GC, 314 ORFs), with neither genome carrying antimicrobial resistance, virulence or lysogeny-associated genes. VIRIDIC and VICTOR analyses placed HHH within Kayfunavirus and RZH within Karamvirus, supporting their classification as distinct species. Both phages demonstrated rapid adsorption, short latent periods and stability across physiological pH and temperature ranges. A phage cocktail targeting MDR ECC strain was evaluated with EDTA against established biofilms. Crystal violet assays showed the greatest biomass reduction at MOI 10 with 0.5-0.75 mM EDTA. Bliss independence analysis revealed localized synergy within this window but significant overall antagonism at higher EDTA concentrations. CFU enumeration confirmed greater activity against 24 h than 48 h biofilms. The optimized combination also reduced recoverable bacteria in a fibroblast infection model while maintaining low LDH release. These findings identify two novel lytic Enterobacter phages and support a narrow EDTA concentration window for enhanced phage-mediated antibiofilm activity.

Biofilms↗
PubMed2026

Genome-wide characterisation of the myosin light chain gene family in Chinese perch (Siniperca chuatsi) and its expression patterns in muscle fibre types and injury response.

The Class II myosin light chain (myl) genes in Chinese perch (Siniperca chuatsi) have not yet been systematically characterised, and relationships with muscle fibre specification, development, and injury-associated remodelling remain unclear. In this study, fast and slow muscle fibres were initially distinguished using myofibrillar ATPase histochemistry. Subsequently, genome-wide mining identified 16 Class II myl genes, comprising eight essential and eight regulatory light-chain subunits. Their conserved-domain features, chromosomal distribution, phylogenetic relationships and expression profiles were analysed. Transcriptomic profiling showed that summed myl transcript abundance was higher in fast muscle than in slow muscle, accounting for 67.2% of the pooled myl transcript pool across the two muscle types (paired t-test, raw P = 0.036). mylpfa, myl1 and mylz3 were the major fast-muscle-associated genes, whereas myl10, myl2b and myl13 were preferentially expressed in slow muscle at the transcript level. These patterns support these genes as candidate fibre-type-associated expression markers. Developmental profiling identified stage-associated myl expression patterns, including a possible expression shift between mylpfb and mylpfa. In the descriptive injury-repair time course (d0-d7), FPKM profiles indicated that fast-muscle-associated genes (mylpfa, mylz3 and myl1) were lower at d1 and recovered by d3, whereas several slow-muscle-associated genes showed biphasic transcript-level increases. The slow-muscle-associated RLC gene mylpfb showed a delayed expression peak at d7. Notably, the embryonic isoform myl6l showed a modest increase from approximately 2 FPKM at d0 to 4-5 FPKM after injury, suggesting a possible injury-associated expression pattern that requires further validation. Together, these findings provide a genome-wide description of the Chinese perch myl gene family and identify candidate fibre-type-associated genes and descriptive injury-associated isoform expression patterns.

Animals↗
PubMed2026

Genome-wide characterization of heat shock protein genes reveals thermal stress-responsive candidates in Litopenaeus vannamei.

Heat shock proteins (HSPs) are conserved molecular chaperones involved in protein folding, refolding, aggregation prevention, and degradation of damaged proteins. However, the genomic organization and thermal responsiveness of HSP genes in the Pacific white shrimp (Litopenaeus vannamei) remain incompletely understood. Here, we performed a genome-wide analysis of the HSP gene family and examined its phylogenetic relationships, structural features, duplication patterns, sequence variation, interaction networks, and transcriptional responses to acute heat stress. A total of 34 HSP genes were identified and classified into the HSP90, HSP70, HSP40/DNAJ, HSP60, and small HSP families. Phylogenetic, motif, gene structure, synteny, and subcellular localization analyses revealed evolutionary conservation and structural diversification among family members. Three duplicated gene pairs were identified, comprising two segmental duplications and one tandem duplication. All pairs exhibited Ka/Ks ratios below 1, consistent with purifying selection of varying strength. Sequence analysis identified 295 nonsynonymous single-nucleotide polymorphisms, of which 12 were consistently predicted to be deleterious by multiple algorithms. Protein-protein interaction analysis indicated enrichment of protein-folding and cellular stress-response functions. RT-qPCR analysis showed significant induction of HSPA4, HSP90AA1, TRAP1, BiP, and DNAJA1 after 6, 12, and 24 h of exposure to 34 °C, whereas DNAJC3 was significantly induced only at 12 h. All six genes reached their highest transcript abundance at 12 h. These findings may provide a genomic framework for HSP genes in L. vannamei and identify candidate genes and variants associated with thermal stress responses.

Animals↗