SEARCH · Search PubMed

Results for “Functional genes”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

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

1,017 records · Page 3Linked to original sources

PubMed2026

Deciphering CD8+ T cell exhaustion in human cancers through single-cell and spatial transcriptomics.

Exhausted CD8+ T cells (Tex) within the tumor microenvironment (TME) represents a critical barrier limiting anti-tumor immune responses. Tex cells are characterized by upregulated inhibitory immune checkpoint receptors, reduced cytotoxicity, and functional heterogeneity. Their genomic features and regulatory networks remain poorly defined, and only a minority of patients respond to immune checkpoint blockade (ICB) therapy. Single-cell RNA sequencing (scRNA-seq), through high-resolution transcriptomic profiling, has revealed diverse Tex subpopulations, identified subpopulation-specific marker genes and regulatory pathways. Spatial transcriptomics has further mapped the spatial distribution of Tex and their interaction networks with immune cells, tumor cells, and stromal cells, elucidating the impact of spatial heterogeneity on Tex functionality. Current studies indicate that the exhausted state of Tex is dynamic and modifiable, with functional differences among subpopulations closely associated with tumor progression and therapeutic response. However, the genomic characteristics, epigenetic regulation, and spatial interaction mechanisms of Tex require further exploration. This review summarizes recent advances in high-resolution omics technologies for precisely dissecting Tex heterogeneity, functional features, and interactions with other cells. It emphasizes the central value of optimizing Tex-targeted tumor immunotherapy strategies, providing theoretical foundations and directional guidance for developing more effective anti-tumor immunotherapies.

Humans↗

PubMed2026

A machine learning-derived and functionally validated circadian rhythm signature predicts clinical outcomes and in silico drug sensitivity in colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) displays considerable heterogeneity in clinical outcomes, highlighting the need for reliable prognostic biomarkers. While the aberrant expression of circadian rhythm-related genes has been implicated in cancer pathogenesis, its comprehensive role in CRC progression and predicted therapeutic vulnerabilities remains inadequately characterized. METHODS: Bulk and single-cell RNA-sequencing data were integrated from multiple CRC cohorts. A circadian rhythm signature (CRS) was developed through machine learning algorithms and validated for prognostic value. Comprehensive analyses of tumor microenvironment, genomic alterations, and drug sensitivity were performed. Furthermore, the biological function of the core gene, BHLHE40, was validated in CRC cell lines through CCK-8, EdU, and wound healing assays. RESULTS: Single-cell analysis demonstrated an elevated expression signature of circadian rhythm-related genes in dendritic cells. The optimized CRS, comprising 14 circadian rhythm-related genes, successfully categorized patients into high- and low-risk groups. Patients with a high CRS showed markedly poorer overall survival and computationally inferred immunosuppressive features, including reduced CD8+ T cell infiltration and increased M2 macrophage polarization. Genomic analysis revealed enhanced mutation burden in TP53 and alterations in RTK-RAS/WNT pathways. Notably, in vitro assays confirmed that BHLHE40 is significantly overexpressed in CRC cells. Knockdown of BHLHE40 markedly inhibited tumor cell proliferation and migration. Drug sensitivity profiling identified bexarotene and SMER-3 as potential therapeutic options for high-CRS patients. A nomogram integrating CRS with clinical parameters demonstrated superior predictive accuracy for 1-, 3-, and 5-year survival. CONCLUSIONS: The CRS represents a promising prognostic biomarker that reflects tumor immune status and genomic features, providing valuable insights for personalized treatment strategies in CRC.

Circadian rhythm↗

PubMed2026

Statistical test to compare the linkage model and the admixture model based on central limit results.

In the Admixture Model, the probability that an individual carries a certain allele at a specific marker depends on the allele frequencies in K ancestral populations and the proportion of the individual's genome originating from these populations. The markers are assumed to be independent. The Linkage Model is a Hidden Markov Model that extends the Admixture Model by incorporating linkage between neighboring loci. We prove consistency and asymptotic normality of maximum likelihood estimators for the ancestry of individuals in the Linkage Model, complementing earlier results by (Pfaff et al., 2004; Pfaffelhuber and Rohde, 2022; Heinzel, 2025) for the Admixture Model. These results are used to prove that a statistical test that allows for model selection between the Admixture Model and the Linkage Model is an asymptotic level-α-test. Finally, we demonstrate the practical relevance of our results by applying the test to real-world data from The 1000 Genomes Project Consortium (2015).

Genetic Linkage↗

PubMed2026

Transcriptomic responses of Porphyrophora sophorae larvae during licorice root colonization reveal coordinated remodeling of translation, mitochondrial energy metabolism and defense-related genes.

BACKGROUND: Porphyrophora sophorae is a subterranean piercing-sucking scale insect that damages licorice (Glycyrrhiza uralensis) roots, but the molecular responses associated with larval root colonization remain insufficiently defined. METHODS: We compared non-parasitic larvae (NP) and root-colonizing larvae (RC) using six RNA-seq libraries, de novo transcriptome assembly, DESeq2-based differential expression analysis, GO/KEGG enrichment, annotation-based candidate gene screening, and RT-qPCR validation of selected genes. RESULTS: Sequencing yielded 260.91 million clean reads, and de novo assembly produced 60,794 non-redundant transcripts. DESeq2 identified 703 FDR-significant DEGs, including 49 upregulated and 654 downregulated genes in RC larvae. Upregulated genes were mainly associated with translation- and ribosome-related processes, whereas downregulated genes were enriched in mitochondrial, oxidation-reduction, energy metabolism, and oxidative phosphorylation-related functions. Annotation-based screening identified 75 FDR-significant candidate genes associated with chemosensation, defense-related responses, and energy metabolism, with mitochondrial energy metabolism-related genes forming the largest module. RT-qPCR validation based on the raw Ct data showed concordant expression directions for ten selected transcript targets. CONCLUSIONS: Root colonization in P. sophorae larvae was associated with coordinated transcriptional remodeling involving selective activation of translation-related processes, adjustment of mitochondrial energy metabolism, and changes in defense-related gene expression. These results provide candidate molecular targets for future functional studies of host contact, feeding establishment, and physiological adjustment in this subterranean scale insect.

Animals↗

PubMed2026

Investigation of Fatty Acid Metabolism-Associated Molecular CPOX and the Underlying Mechanism in Follicular Lymphoma.

Dysregulated lipid metabolism is a key driver of follicular lymphoma (FL). This study aimed to explore the lipid metabolism-related genes (LMRGs) and clarify the underlying roles and mechanisms in FL. Bioinformatics methods, including differential analysis, WGCNA, machine learning, and Mendelian randomization, were utilized to select the LMRGs in FL. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were conducted to investigate the function of the key LMRG. Receiver operator characteristic (ROC) was used to evaluate the diagnostic value of the key gene CPOX. A pan-cancer analysis investigated CPOX's expression level and immune correlations. In vitro experiments using FL cell lines (WSU-FSCCL, DOHH2) validated CPOX expression, and CPOX knockdown in DOHH2 cells was used to assess its impact on viability, migration, invasion, and fatty acid metabolism. CPOX was confirmed to be a risk factor, significantly overexpressed in FL, and exhibited effective diagnostic ability in FL (AUC = 0.731). Functional analysis linked CPOX to mitochondrial function, oxidative phosphorylation, and heme metabolic process. Pan-cancer indicated the dysregulated CPOX across multiple cancers and closely correlation with immune characteristics. Experimentally, CPOX was higher in the more invasive DOHH2 cells; and CPOX knockdown suppressed FL progression and reduced lipid droplet formation, triglyceride, total cholesterol, and free fatty acid levels. In conclusion, this study fills the gap in understanding the significance of lipid metabolism-related molecules in FL, and innovatively proposes that CPOX is a risk factor for FL. Knockdown of CPOX inhibits the FL progression, which is regulated by fatty acid metabolism.

Lymphoma, Follicular↗

PubMed2026

ScRNA-seq analysis reveals the effects of nitrite stress on the endocrine system of the eyestalk in Litopenaeus vannamei.

Nitrite is a harmful substance generated in Litopenaeus vannamei farming systems, largely originating from the inadequate breakdown of surplus feed and shrimp feces. Its accumulation in the water can affect the growth and physiological functions of shrimp, damage the immune system, and even cause mass mortality, thus becoming a key environmental factor restricting the green development of the industry. Under nitrite stress, the eyestalk, as an important neuroendocrine regulatory center in crustaceans, participates in the stress adaptation of the organism and exerts a protective effect by regulating energy metabolism and immune function. However, the molecular regulatory mechanism of the eyestalk in response to nitrite stress remains unclear. In this study, single-cell RNA sequencing (scRNA-seq) technology was used to analyze the heterogeneity of eyestalk cells in L. vannamei under nitrite stress. A total of 18, 394 high-quality cells were obtained, and six major cell subpopulations, including Neurosecretory cell, Motor neuron, Sensory neuron, Interneuron, Neurogliocyte, and Support cell, were identified. Differential expression analysis identified 839 differentially expressed genes, and different cell types showed distinct specific responses to nitrite stress. Functional enrichment analysis indicated that pathways such as glycolysis, oxidative phosphorylation, ribosome function, and endoplasmic reticulum protein processing were significantly activated, while signal transduction and DNA repair-related pathways were inhibited. Further analysis revealed that nitrite stress could induce mitochondrial function changes and trigger oxidative stress, thereby affecting the neuroendocrine system function of the eyestalk. This study provided insights into transcriptomic responses of the eyestalk to nitrite stress at the single-cell level, laying a theoretical foundation for the management of aquaculture environments.

Animals↗

PubMed2026

Assembly and Characterization of the First Complete Mitochondrial Genome of Tussilago farfara L.: Insights into Biological Functions and Phylogenetic Relationships within the Asteraceae Family.

Tussilago farfara L., a member of the Asteraceae family, is an economically valuable species due to its edible and medicinal properties. To elucidate the structural characteristics, genetic mechanisms, and evolutionary pathways of the organelle genomes of T. farfara, we sequenced, assembled, and annotated its mitochondrial genome for the first time. The complete mitochondrial genome of T. farfara spans 306,024 bp and contains 33 mitochondrial protein-coding genes (PCGs), 3 rRNAs, and 22 tRNAs. Analysis of the nucleotide substitution rate and genetic diversity revealed that most mitochondrial genome genes may have undergone purifying selection, indicating a slow evolutionary rate and a relatively conserved genomic structure. We further identified 13 fragments of chloroplast-derived DNA integrated into the mitochondrial genome, evidencing intracellular gene transfer. Collinearity analysis showed that Arctium lappa shares the most extensive mitochondrial homologous sequences and the highest sequence similarity with T. farfara. Phylogenetic analysis based on the mitochondrial genome helped to clarify the evolutionary and taxonomic position of T. farfara within the Asteraceae family. The mitochondrial genome sequence of T. farfara provides a valuable genomic resource for species identification and for evolutionary studies within the Asteraceae family.

Genome, Mitochondrial↗

PubMed2026

Whole-exome characterization of host genetic variation in HIV-associated genes across the high-prevalence Mizo population, Northeast India.

BACKGROUND: The Mizoram state of Northeast India has one of the highest HIV prevalence rates in Asia, yet the host genetic factors influencing HIV susceptibility in this Tibeto-Burman population remain uncharacterised. METHODS: We performed whole-exome sequencing using Illumina NovaSeq 6000, mean coverage 100X on 76 HIV-negative Mizo individuals. Variants were called using GATK HaplotypeCaller v4.3 against GRCh38p14, annotated with ANNOVAR, and filtered using hard-quality thresholds (QD ≥ 2, SOR ≤ 3, MQ ≥ 40, DP ≥ 10, GQ ≥ 20). The allele frequencies were compared against gnomAD v2.1.1 population databases. Hardy-Weinberg equilibrium was assessed using the Wigginton exact test with Bonferroni correction. RESULTS: Post-quality filtering resulted in 12,011 sample-variants across 2,821 unique positions from 36 HIV-associated loci (33 protein-coding genes, 2 chemokine ligands, and 3 lncRNA targets). Of these, 784 observations (51 unique positions) were high-impact nonsynonymous or loss-of-function variants. ADAR rs2229857 (p.K384R, NM_015840) was the most frequently observed variant (Mizo carrier frequency = 0.895; 95% CI: 0.806-0.946). CXCR1 rs16858808 (p.R335C) showed the greatest population enrichment (Mizo carrier frequency = 0.197; 95% CI: 0.123-0.300; 7.65-fold carrier-frequency enrichment versus gnomAD South Asian; CADD = 15.60). Sixteen of 20 tested variants deviated from Hardy-Weinberg equilibrium after Bonferroni correction (p < 0.0025), predominantly showing excess homozygosity consistent with the endogamous Mizo population. The protective variant CCR5-Δ32 was absent in all the 76 individuals tested. CONCLUSION: This first whole-exome characterization of HIV host genes in the Mizo population identifies CXCR1 rs16858808 as the most population-enriched functional variant and reveals a pervasive endogamy signature. These findings provide a population-specific genetic framework for future HIV susceptibility studies and ART pharmacogenomics research.

Humans↗

PubMed2026

Correlative analysis of endogenous miRNA expression profiles underlying brown planthopper adaptation to resistant rice.

The brown planthopper (Nilaparvata lugens Stål, BPH) is a major insect pest threatening global rice production. However, the molecular mechanisms underlying the adaptation of BPH populations with different virulence levels to resistant rice cultivars remain poorly understood. MicroRNAs (miRNAs), as key post-transcriptional regulators, play critical roles in host adaptation in herbivorous insects. In this study, we analyzed the miRNA expression profiles of a high-virulent population (IR56p) and a low-virulence population (TN1p) after feeding on susceptible (TN1) and resistant (IR56) rice cultivars. Our findings reveal distinct miRNA-mediated regulatory strategies employed by the two populations. The IR56p population showed downregulation of miRNAs including miR-10, miR-124, and miR-316, showing an inverse correlation with increased expression of predicted target genes involved in detoxification (carboxylesterase, UDP-glycosyltransferase) and effector function (calmodulin). In contrast, several miRNAs highly expressed in IR56p, including miR-307, miR-317, and miR-275, were predicted to target rice genes associated with hormone signaling, cell wall biosynthesis, and oxidative homeostasis, suggesting a possible but unproven inter-species regulatory role that requires functional validation. Collectively, these descriptive and correlative findings provide hypothesis generating insights into insect-plant coevolution and identifies candidate molecular targets for future functional validation and RNA interference-based pest management strategies.

Animals↗

PubMed2026

Genome-wide identification and expression profiling of HSD3B and SDR42E1 genes in the Pacific oyster (Crassostrea gigas): potential associations with gonadal development.

Sex steroids are lipid-soluble signaling molecules that regulate sex differentiation, reproductive development and physiological homeostasis in animals. 3β-Hydroxysteroid dehydrogenase/Δ5-Δ4 isomerase (3β-HSD) is a key steroidogenic enzyme, whereas SDR42E1, an extended short-chain dehydrogenase/reductase, has been implicated in sterol- and steroid-related metabolism. However, the composition, evolutionary relationships and expression patterns of the HSD3B- and SDR42E1-related genes in bivalve gonadal development remain poorly characterized. In this study, five PF01073-containing genes, comprising three CgHsd3b and two CgSdr42e1 genes, were identified in the Pacific oyster Crassostrea gigas. Phylogenetic analysis separated the proteins into HSD3B-related and SDR42E1-related groups, and gene-structure and motif analyses indicated subfamily-level divergence. All five proteins retained the SDR domain but differed in exon-intron structure and motif composition. Each contained the extended-SDR TGxxGxxG motif, whereas exact classical [ST]GxxxGxG and NNAG motifs were absent. Tyr- and Lys-equivalent residues were conserved, while the HSD3B1 Ser-equivalent position contained Thr in two C. gigas proteins and Ser in one. These features support their classification as extended-SDR proteins but do not establish enzymatic activity or substrate specificity. The three CgHsd3b genes were dispersed on one chromosome, whereas CgSdr42e1-1 and CgSdr42e1-2 were adjacent on another chromosome, suggesting a possible local duplication event for the CgSdr42e1 pair. Public RNA-seq data showed distinct tissue- and gonadal-stage expression patterns, with several genes displaying gonad-biased or female-stage-associated expression. Independent RT-qPCR profiling of the representative genes CgHsd3b-3 and CgSdr42e1-1 detected stage-dependent expression, although tissue rankings differed from those in the public RNA-seq datasets. These differences may reflect the use of independent biological samples, tissue composition, normalization procedures, and platform-specific measurements. Because enzymatic assays, metabolite measurements, cellular localization, and functional perturbation were not performed, the results identify candidate genes whose expression is associated with gonadal development rather than demonstrating regulatory roles. This study provides a comparative framework for future functional investigation of sterol- and steroid-related metabolism in bivalves.

Animals↗

PubMed2026

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

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

Animals↗

PubMed2026

Novel environmental contaminant 6PPD-quinone promotes malignant phenotypes in colorectal cancer cells and identifies candidate response-associated genes.

6PPD-quinone (6PPDQ), an oxidative transformation product of the widely used tire antioxidant 6PPD, is a ubiquitous environmental contaminant with bioaccumulation potential and widespread human exposure. Recent epidemiological evidence indicates a positive association between urinary 6PPDQ levels and colorectal cancer (CRC) risk; however, its biological effects on CRC-related phenotypes and associated molecular responses remain unclear. We integrated bioinformatics analysis, prognostic modeling, molecular docking and dynamics simulations, and in vitro experiments to investigate cellular and molecular responses to 6PPDQ in CRC models. Predicted 6PPDQ targets were intersected with CRC prognosis-related genes from The Cancer Genome Atlas, followed by functional enrichment and LASSO regression to construct a prognostic risk model, with 1-, 3-, and 5-year AUC values of 0.727, 0.754, and 0.778, respectively. Molecular docking and 100-ns molecular dynamics simulations suggested interactions between 6PPDQ and candidate proteins, including CPT2, SHC2, SRMS, and STK35. Functional assays showed that 6PPDQ exposure altered proliferation, wound-closure capacity, and invasion in Caco-2 and HCT116 cells across the nanomolar concentration range, with non-monotonic and cell-line-dependent responses. In contrast, NCM460 cells showed no increase in EdU incorporation at 10 or 100 nM, whereas reduced proliferation at higher concentrations was accompanied by increased LDH release. 6PPDQ also altered the expression of several prognosis-associated candidate genes. These findings identify cellular phenotypes and candidate molecular responses associated with 6PPDQ exposure under the tested in vitro conditions, but do not establish their causal roles or in vivo relevance. Further mechanistic and in vivo studies are required.

Humans↗

PubMed2026

Genome mining of alkaliphilic cyanobacterial consortia: identification of biosynthetic gene clusters in Sodalinema and associated heterotrophs.

Alkaline soda lakes are high-pH environments that host specialized microbial communities with potential for biotechnology and natural product discovery. We characterized three Sodalinema-dominated cyanobacterial consortia enriched from Canadian soda lakes over 510 days. Using hybrid metagenomic sequencing and metatranscriptomics across pH, alkalinity, and temperature gradients, we reconstructed high-quality metagenome-assembled genomes and assessed functional activity. All consortia converged toward cyanobacteria dominance and exhibited temperature optima between 21°C and 30°C. Phylogenetic analysis placed Sodalinema genomes within a distinct clade affiliated with Candidatus Sodalinema alkaliphilum. Genomic analysis indicated complete biosynthetic pathways for vitamin B5, vitamin B7, and the molybdenum cofactor, but incomplete pathways for vitamins B1, B9, and B12, consistent with patterns observed in Sodalinema yuhuli. Metatranscriptomic profiles showed increased expression of genes involved in phycocyanin and carotenoid biosynthesis at pH 10.2 relative to pH 8.5. Biosynthetic gene cluster analysis revealed that most secondary metabolic potential resided in heterotrophic community members. Roseinatronobacter encoded pathways for N-acyl homoserine lactones, osmoprotectants, betalactones, and prodigiosin, while Alkalimonas, Wenzhouxiangella, and members of the Kiloniellales encoded clusters for lanthipeptides, cyclodipeptides, hydrogen cyanide, and pyrroloquinoline quinone. These findings indicate functional partitioning within the consortia and highlight the contribution of heterotrophs to secondary metabolism.IMPORTANCEAlkaline soda lakes contain microbial communities adapted to high pH that remain underexplored for biotechnology. This study focuses on Sodalinema, a filamentous cyanobacterium that dominates enriched consortia from Canadian soda lakes, and its associated heterotrophic partners. We show that while Sodalinema drives primary productivity, heterotrophic bacteria encode most of the pathways for antimicrobial and signaling compounds. These interactions may support community stability and defense against competing microorganisms. By linking genomic potential with gene expression, this work identifies alkaline cyanobacterial consortia as a source of bioactive compounds and provides a framework for exploring extremophilic microbial communities for natural product discovery.

Sodalinema↗

PubMed2026

Searching for New Genes That Cause Usher Syndrome.

PURPOSE: The purpose of this project was to identify novel Usher syndrome (USH) candidate genes from phenotyping data of 9139 knockout (KO) mouse lines. METHODS: We evaluated phenotype data for concurrent retinopathy and hearing abnormalities in single-gene KO mice generated by the International Mouse Phenotyping Consortium (IMPC). A search was performed to determine whether each gene had been previously associated with retinopathy and/or deafness in humans. Bioinformatic tools were used to predict protein interactions, molecular functions, signaling pathways, and the expression of human orthologues of candidate genes in the retina and inner ear. RESULTS: We identified 18 single-gene KO lines exhibiting hearing abnormality and retinopathy after ear and eye examinations, respectively, and/or by histopathology. The molecular functions and signaling pathways of the human orthologues of the 18 candidate genes partially overlapped with those of USH genes. Particularly, FER and DYRK1B proteins were predicted to interact with proteins encoded by known ciliopathy genes. ADIPOR1, ATP8B1, and MPDZ were associated with retinal degeneration in humans. CHSY1 and IDUA may be pathogenic causes of hearing impairment in people. Furthermore, CHSY1, CSTB, and SPRED1 were located adjacent to unsolved genetic loci related to USH. CONCLUSIONS: A screen of 9139 KO mouse lines revealed 18 candidate genes exhibiting both retinal and inner ear abnormalities consistent with the principal clinical features associated with USH. As the observed phenotypes are attributed to gene deletion in mice, these genes warrant further study to determine the causation of retinal degeneration and hearing loss in patients.

Animals↗

PubMed2026

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

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

Animals↗

PubMed2026

Male accessory gland proteins in Grapholita molesta: Identification and reproductive functional validation of four accessory gland-specific lipases.

Accessory gland proteins (Acps), synthesized in the male accessory glands (AGs), are transferred to females via spermatophores during mating and elicit diverse post-mating physiological and behavioral responses. However, Acps have not been comprehensively characterized in Grapholita molesta, a cosmopolitan orchard pest. Here, using data-independent acquisition mass spectrometry, we describe an integrated proteomic approach combining comparative AG analyses (virgin vs. newly mated) with spermatophore profiling to identify Acps in G. molesta. According to the established screening criteria, we identified 83 confirmed Acps, which were classified into nine categories. Tissue-specific expression patterns of 20 randomly selected Acp genes were evaluated, revealing that these genes were specifically or highly expressed in male AGs. Among the 83 confirmed Acps, four Acps harbored the PLN02872 superfamily domain and were classified into the canonical lipase family. Notably, their transcripts were all highly expressed in the AGs during the pre-maturation stage. These four Acps were selected for preliminary validation of their male reproductive functions. RNAi-mediated knockdown of three out of four lipase genes in G. molesta males significantly decreased the fertility of mated females, with phenotypes including a significant reduction in egg production and egg hatching rate. This study provides a comprehensive catalog of high-confidence Acps, lays a foundation for subsequent in-depth functional characterization of these reproductive proteins, and offers promising molecular targets for the development of novel genetic regulation-based integrated pest management strategies.

Animals↗

PubMed2026

Genome-wide identification, structural characterization, and evolutionary analysis of growth-related gene families in African catfish (Clarias gariepinus).

The somatotropic axis encompassing growth hormone (GH), insulin-like growth factor (IGF), myostatin (MSTN), and prolactin (PRL) signalling cascades is the master regulator of somatic growth, metabolism, and development in vertebrates. African catfish (Clarias gariepinus), a commercially pivotal aquaculture species, now possesses a chromosome-level reference genome (CGAR_prim_01v2); however, a systematic, genome-wide characterization spanning all five interconnected growth-related gene families has not previously been undertaken in this species. Here, we identified and characterized 15 growth-related genes spanning gh1, ghra, ghrb, Igf1, Igf2a, Igf2b, igf1ra, Igf1rb, Igf2r, Mstna, Mstnb, prl, prlra, prlrb, and smtlb distributed across 13 chromosomes. Complete one-to-one orthology with zebrafish confirmed strong dosage-balance conservation across >120 million years of teleost divergence. Physicochemical analysis resolved a clear biochemical dichotomy between compact, basic secreted ligands (19.88-45.81 kDa; pI up to 10.02) and large, acidic, heavily glycosylated membrane receptors (56.82-270.80 kDa; pI 4.85-5.97). Phylogenetic analysis confirmed 3R whole-genome duplication origins for all paralog pairs, while synteny analysis revealed a disruption of the ancestral gh1-prl chromosomal block in C. gariepinus, a finding that warrants further comparative and functional investigation. This genomic atlas provides the sequence and structural information including exon-intron boundaries, domain architecture, and chromosomal coordinates needed as a prerequisite for future marker-assisted selection and CRISPR-based myostatin-editing efforts in African catfish aquaculture, though translation into applied breeding outcomes will require subsequent functional and expression studies.

Animals↗

PubMed2026

Characterization and functional insights of histone deacetylases in bivalves: implications for temperature and immune response in Chlamys nobilis.

Histone deacetylases serve as pivotal epigenetic regulators that modulate chromatin remodeling and gene transcription, playing critical roles in immune defense and environmental stress responses in aquatic organisms. However, the evolutionary characteristics and functional roles of the HDAC family in bivalves remain poorly understood. In this study, genome-wide identification of the HDAC family across 30 bivalve species yielded 558 HDAC genes. Phylogenetic reconstruction categorized these genes into four conserved groups and revealed a unique, bivalve-specific SIRT8 clade. Using the noble scallop Chlamys nobilis as a representative model, expression profiling revealed distinct expression patterns among CnHDAC members. Class I and most Class III members were predominantly expressed in the gonads, while Class II members were enriched in immune-related tissues, implying their potential involvement in bivalve immunity. Upon temperature stress, CnHDAC1/2, CnHDAC11-1, CnHDAC11-2, CnSIRT2-1, CnSIRT4, CnSIRT6, and CnSIRT8-3 were significantly induced, highlighting their critical roles in temperature adaptation. Upon Vibrio exposure, CnHDAC1/2, CnHDAC8, CnSIRT4, and CnSIRT6 were upregulated, while CnHDAC4/5/7/9, CnHDAC6/10, CnSIRT2-2, CnSIRT5, CnSIRT7, and CnSIRT8-3 were downregulated, suggesting a coordinated epigenetic regulatory mechanism underlying host immune defense. In conclusion, this study systematically elucidates the evolutionary landscape of the HDAC family and underscores its potential involvement in environmental resilience and host immunity, providing a theoretical basis for the breeding of disease-resistant and stress-tolerant aquaculture bivalves.

Animals↗