SEARCH · Search PubMed

Results for “Sequence Analysis, RNA”

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,667 records · Page 12Linked to original sources

PubMed2026

Exercise-associated epigenetic remodeling and TCR repertoire dynamics in Lynch syndrome carriers.

Lynch syndrome (LS) carriers are at elevated cancer risk. Emerging evidence suggests that exercise may serve as a non-pharmacologic preventive strategy, yet the epigenetic and immunological mechanisms underlying its protective effects in this population remain unclear. Here, we perform integrative multi-omics profiling of DNA methylation, gene expression, and the T cell receptor (TCR) repertoire in LS carriers undergoing a 52-week aerobic cycling intervention. We identify compartment-specific DNA methylation changes, including innate immune activation in cfDNA and oncogenic pathway repression in tissue. Integrative transcriptomic analysis highlights ISL1 as a key exercise-repressed, epigenetically regulated gene, and identifies FLCN as a colorectal cancer (CRC)-associated methylation target. TCR analysis reveals an exercise-associated increase in systemic repertoire diversity and tissue-specific clonal convergence, thus suggesting antigen-driven recruitment. Collectively, these findings uncover epigenetic and immune remodeling as potential mechanisms of exercise-mediated protection in LS.

Lynch syndrome↗

PubMed2026

Integrating genomic distance analyses in the description of a new family, genus, and species of sponge-associated antipatharians (black corals).

Antipatharians (black corals) are among the least studied coral groups, with much of their diversity still undescribed. Here, we present an integrative morphological, phylogenomic and genomic distance study of deep-sea antipatharians sampled in high seas areas of the North Pacific Ocean and from New Zealand's Exclusive Economic Zone. These corals grow on hexactinellid sponges - a unique characteristic in the order Antipatharia. Using a dataset of ultra-conserved elements and exons, combined with morphological analyses, we reconstruct phylogenomic relationships and formally describe a new family (Eidikopathidae fam. nov.), a new genus (Eidikopathesgen. nov.), and two new species (E. korallispongiasp. nov., E. zealandkoralliasp. nov.). Morphologically, the new family is distinguished by a corallum consisting of a network of loose branches that fuse with the sponge skeletal framework. Phylogenomic analyses recovered consistent topologies with strong nodal support, corroborating the distinct evolutionary placement of this sponge-associated lineage. Pairwise genomic distances estimated using the Tamura-Nei model were concordant with patristic genomic distances, identifying Pteridopathidae as the genetically closest family to Eidikopathidae fam. nov., followed by Myriopathidae and Stylopathidae, which were recovered as sister families in the phylogeny. This pattern shows that genomic distance complements, rather than simply mirrors, tree topology by quantifying accumulated sequence divergence among lineages. Together, these results provide the first genomic distance framework for Antipatharia, offering a baseline for future systematic, evolutionary, and biodiversity studies on this fundamental shallow, mesophotic and deep-sea coral group.

Animals↗

PubMed2026

Tigecycline-resistant Staphylococcus in waiting pens of a pig slaughterhouse: genomic insights into a food safety alert.

BACKGROUND: The waiting pens of slaughterhouses represent a critical control point in the 'farm-to-fork' continuum, yet their role in the emergence and dissemination of antimicrobial resistance remains understudied. This study investigated tigecycline-resistant Staphylococcus (TRS) in these high-risk zones to assess their prevalence, resistance mechanisms, and transmission dynamics. METHODS: 400 samples were collected from the waiting pens of a pig slaughterhouse in Guangzhou, China. Antimicrobial susceptibility testing, whole-genome sequencing, phylogenetic analysis, and molecular cloning were employed to characterize resistance mechanisms and transmission patterns. RESULTS: 78 TRS strains were isolated and classified into three species, including S. borealis, S. ureilyticus, and S. pasteuri. These isolates exhibited multidrug-resistant phenotypes and carried new mutations in rpsJ and tet(M), which were functionally confirmed to reduce tigecycline susceptibility. Phylogenetic evidence demonstrated clonal transmission between pig farms and the slaughterhouse. The tet(M) gene was located within Staphylococcal cassette chromosome mec elements mediated by IS257, while tet(L) was carried by plasmids formed through IS256/IS257-mediated recombination. CONCLUSIONS: Waiting pens serve as crucial reservoirs for the amplification and dissemination of antimicrobial resistance. Our findings underscore the urgent need for enhanced biosecurity measures, improved waste management, and routine molecular surveillance in these high-risk zones to mitigate the spread of resistance along the food production chain.

Animals↗

PubMed2026

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

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

Animals↗

PubMed2026

Lower-limb electromyographic responses across four recumbent cycling configurations in healthy young men.

Recumbent cycling provides an adjustable configuration in which backrest angle and crank-axis height can alter lower-limb loading. Backrest angle changes the orientation of the active musculature relative to gravity, whereas crank-axis height modifies the foot-pedal interface and the ankle lever arm. Whether such positioning is associated with differences in lower-limb electromyographic (EMG) activation across both factors together has not been examined. Using a randomized crossover design, fifty-nine healthy young men cycled at 100 W and 50 rpm under four conditions combining two backrest angles (30°, 75°) and two crank-axis heights (Rh, high, 77 cm; Rl, low, 45 cm). Surface EMG from the rectus femoris (RF), biceps femoris (BF), and gastrocnemius medialis (GAS) was analyzed for root mean square (RMS) amplitude, absolute onset sequence, and relative recruitment timing (%Sequence, derived using Telescan software and referenced to RF). RMS amplitude did not differ across conditions for any muscle (all p ≥ 0.356), and absolute onset sequence showed no condition effect (all p ≥ 0.197). For GAS %Sequence, the Greenhouse-Geisser-corrected omnibus test was significant (F(1.32, 76.56) = 28.45, p < 0.001), driven by a within-30° contrast in which %Sequence was lower under Rl than Rh (paired t(58) = 10.39, Holm-adjusted p < 0.0001; mean difference 22.0 percentage points, 95 % CI [17.8, 26.3]; dz = 1.35). Because %Sequence is referenced to RF, this exploratory shift should be interpreted cautiously. Lower-limb activation may be constrained by the closed-chain nature of pedaling, whereas distal recruitment timing may be sensitive to crank-axis height and the foot-pedal interface.

Humans↗

PubMed2026

Genome-wide SNP data support species boundaries in sympatric Polylepis Ruiz & Pav. (Rosaceae) species from Bolivia and Ecuador.

Species delimitation in the South American genus Polylepis is notoriously challenging due to high morphological similarity and phenotypic plasticity, likely driven by hybridization and gene flow. Previous phylogenetic studies suggested that genetic structure aligns more strongly with geography than with taxonomy, questioning existing species concepts and hampering conservation efforts. We used double-digest RAD sequencing (ddRADseq) to generate genome-wide SNP data for 11 Polylepis species sampled across multiple localities in Bolivia and Ecuador. Population genetic analyses, phylogenetic inference, and network approaches were combined to assess whether genetic structure aligns more closely with taxonomy or geography. Morphologically defined species formed largely cohesive genetic lineages across regions, with species identity explaining substantially more genetic variation than locality. While localized admixture and reticulation were detected among closely related taxa, widespread species showed strong genetic cohesion and clear separation from congeners. Our results indicate that the sampled Polylepis species from Bolivia and Ecuador maintain distinct genetic identities despite localized signals consistent with gene flow. This genome-wide support for current taxonomy highlights Polylepis as a valuable model for studying speciation under gene flow and indicates that multiple geographic sampling will be essential in reconstructing a robust phylogeny of the genus, with important implications for conservation planning in Andean montane forests.

Bolivia↗

PubMed2026

Molecular Landscape and Advanced Diagnostic Technologies for BRAF Mutations in Cancer: From Quantitative PCR and ddPCR to CRISPR-Based Platforms.

BRAF mutations are key oncogenic alterations across multiple malignancies, including melanoma, thyroid carcinoma, colorectal cancer, non-small cell lung cancer, glioma, and hairy cell leukemia. The most prevalent variant, BRAF-V600E, induces constitutive activation of the MAPK signaling pathway, promoting tumor progression and influencing therapeutic responsiveness. Accurate detection of BRAF alterations is therefore essential for molecular classification, prognostic assessment, treatment selection, and resistance surveillance. This review summarizes the molecular heterogeneity of BRAF mutations and critically evaluates current diagnostic methodologies. Conventional approaches such as allele-specific PCR and Sanger sequencing are compared with advanced quantitative platforms, including high-resolution melting analysis, droplet digital PCR, and next-generation sequencing, with emphasis on analytical sensitivity, mutation coverage, and clinical applicability. Emerging technologies such as CRISPR-based assays, rolling circle amplification systems, and nanoparticle-based biosensors and point-of-care diagnostic platforms are also discussed for their potential to enhance ultra-sensitive detection, particularly in liquid biopsy settings. These emerging tools are highlighted for their potential to enable ultra-sensitive, rapid, and decentralized mutation detection, particularly in liquid biopsy settings. Key challenges, including intratumoral heterogeneity, low allele-frequency variants, FFPE-associated artifacts, and clonal evolution under therapeutic pressure, are examined within a translational framework. In addition, we examine critical barriers to clinical implementation, including standardization, cost, and global accessibility of molecular diagnostics, and outline potential solutions through scalable technologies and decentralized testing strategies. We propose that optimal BRAF testing requires a mutation subclass-informed and clinically integrated strategy combining comprehensive baseline profiling with longitudinal molecular monitoring. Future diagnostic paradigms will likely integrate multi-omics data and artificial intelligence (AI)-assisted interpretation to refine precision oncology implementation. Looking forward, we propose that optimal BRAF testing will require integration of multi-omics profiling with AI-assisted interpretation, enabling automated variant classification, real-time clinical decision support, and improved prediction of therapeutic response and resistance.

Humans↗

PubMed2026

Efficacy of the NMIC-150 system in identifying extended-spectrum beta-lactamases in clinical isolates.

Extended-spectrum beta-lactamases (ESBLs) are significant contributors to the growing global crisis of antimicrobial resistance. This study evaluated the performance of the NMIC-150 System for susceptibility testing of third-generation cephalosporins (3GCs) and assessed whether ceftazidime-avibactam and aztreonam-avibactam could identify ESBL-producing carbapenem-resistant Enterobacterales (CREs). A total of 278 non-duplicate clinical isolates (Klebsiella pneumoniae, E. coli, and Proteus mirabilis) were analyzed. Antimicrobial susceptibility was determined using reference broth microdilution (BMD) and the NMIC-150 System. ESBL production was defined as an ≥eight-fold reduction in the minimum inhibitory concentration (MIC) of 3GCs in the presence of clavulanic acid, according to CLSI criteria. Whole-genome sequencing was performed to characterize ESBL and carbapenemase genes among 3GC-resistant isolates. A Random Forest model was used to predict ESBL-producing isolates based on MIC values. The NMIC-150 System demonstrated over 90% categorical and essential agreement with BMD for ceftazidime and ceftriaxone, along with robust predictive performance via Random Forest analysis. These findings suggest that the NMIC-150 System is a reliable platform for 3GC susceptibility testing and that an ≥eight-fold MIC reduction with ceftazidime-avibactam or aztreonam-avibactam may serve as a phenotypic indicator of ESBL production in CRE isolates. In conclusion, the NMIC-150 System shows potential for routine antimicrobial resistance surveillance and may facilitate the rapid identification of ESBL-producing CREs in clinical settings.

Microbial Sensitivity Tests↗

PubMed2026

Molecular adaptation of caspase genes to salinity stress in the tropical sea cucumber Stichopus monotuberculatus: A comparative analysis across echinoderms.

Apoptosis is an essential physiological process that plays a critical role in development and tissue homeostasis. Caspases, as central regulators of apoptosis, are crucial in controlling inflammation and cell death. In this study, we investigated the caspase gene family in Stichopus monotuberculatus to explore their potential roles in salinity stress adaptation. Five caspase genes were identified from the genome of S. monotuberculatus, including Smcaspase3, Smcaspase6, Smcaspase8a, Smcaspase8b, and Smcaspase8c. Phylogenetic analysis revealed that these Smcaspase genes clustered into distinct caspase subfamilies and showed high conservation with homologs from other echinoderms and representative vertebrates. Conserved motif and gene structure analyses showed relatively similar structural patterns within each clade, whereas divergence was observed among different subfamilies. Promoter analysis identified numerous cis-acting elements related to gene regulation, immune response, and growth and development. Expression profiling under salinity stress showed that Smcaspase8a was significantly upregulated, particularly under prolonged stress, whereas the other genes exhibited limited transcriptional responses. Our findings highlight caspase function in salinity stress and provide the foundation of molecular salinity adaptation mechanisms in S. monotuberculatus.

Animals↗

PubMed2026

Evolutionary characterization and expression profiling of ACC and FASN genes in Chinese mitten crab Eriocheir sinensis.

Acetyl-CoA carboxylase (ACC) and fatty acid synthase (FASN) are rate-limiting enzymes in the fatty acid biosynthetic pathway, yet their evolutionary relationships, sequence features, and expression profiles remain poorly understood in crustaceans, particularly in the economically important Chinese mitten crab (Eriocheir sinensis). Here, we identified and systematically analyzed ACC and FASN genes in E. sinensis using comparative genomic analyses across 43 species. ACC was highly conserved as a single-copy gene in invertebrates, in contrast to the multiple paralogs observed in vertebrates. Similarly, FASN was generally maintained as a single-copy gene across most taxa but exhibited lineage-specific expansion in certain insect groups. Phylogenetic and structural analysis revealed strong conservation of both genes within crustaceans, supported by multiple conserved motifs and canonical functional domains. Expression profiling showed predominant expression in the hepatopancreas and midgut, suggesting their potential involvement in crustacean lipid metabolism. During the molting cycle, ACC and FASN exhibited higher expression levels during stages C and D, suggesting an increased capacity for fatty acid biosynthesis before molting. In addition, dietary lipid levels experiment revealed that ACC and FASN expression responded dynamically to dietary lipid availability, with increased expression at moderate lipid levels but reduced expression under excessive lipid supplementation, indicating a possible adaptive transcriptional response to lipid status. Collectively, this study provides insights into the evolutionary conservation and expression dynamics of ACC and FASN and improves our understanding of lipid metabolic adaptation in crustaceans.

Animals↗

PubMed2026

Whole-Genome Deep Learning Predicts Chemotherapy Response in Colorectal Cancer.

Chemotherapy response in colorectal cancer (CRC) exhibits significant heterogeneity, with current clinical predictors failing to capture complex genomic determinants of resistance. We developed a hybrid deep learning framework integrating convolutional neural networks (CNNs) and bidirectional long short-term memory (BiLSTM) networks to analyze whole-genome somatic mutations, evolutionary conservation, chromatin accessibility, and 3D genome architecture in 2,546 TCGA patients. An attention mechanism identified predictive genomic regions. The model achieved an AUC of 0.92 (95% CI: 0.89-0.94) in cross-validation and 0.88 (95% CI: 0.85-0.91) in independent validation, outperforming clinical models (ΔAUC = +0.18, p < 0.001). Key predictors included non-coding variants in TP53, KRAS, and PIK3CA regulatory regions. Triple-positive patients (mutations in all 3 regions) had significantly worse progression-free survival (HR = 4.7, p < 0.001). Our framework enables accurate chemotherapy response prediction and reveals novel non-coding resistance mechanisms, advancing precision oncology in CRC.

Humans↗

PubMed2026

Identification of elements determining KIR gene demethylation at the CD56-bright stage of NK cell development.

The variegated expression of the KIR family of class I MHC receptors generates specialized natural killer (NK) cells capable of allele-specific HLA recognition. Understanding the mechanism of KIR gene activation will lead to improved methods for the generation of fully functional NK cells. A central RUNX-binding site in the KIR proximal promoter is required for gene activation. RUNX proteins recruit ten-eleven translocation (TET) proteins that generate 5-hydroxymethylcytosine (5hmC) and drive DNA demethylation. Assessment of 5-methylcytosine (5mC) and 5hmC residues at four stages of NK cell development reveals deposition of 5hmC primarily in a CREB site next to the RUNX site at the CD56Bright stage but not the subsequent CD56Dim stage representing fully mature NK cells. KIR promoter demethylation is delayed relative to other lineage-associated genes, indicating a high threshold for KIR gene demethylation in developing NK cells, and a window of opportunity for RUNX/TET-dependent KIR gene activation in CD56Bright NK cells.

6-base sequencing↗

PubMed2026

ChIP-seq profiling identifies diapause-regulated H3K27me3 targets in the fat body of Culex pipiens.

Culex pipiens, a principal vector of significant arboviruses, survives winter through diapause, a hormonally controlled inactive phase that enhances endurance under severe cold circumstances. Recent data suggests that epigenetic processes, namely histone post-translational modifications (hPTMs), play a crucial role in regulating seasonal dormancy. Prior studies from our laboratory indicated a decrease in the methylation of Histone 3 (H3K27me3) in diapausing fat body tissue, associated with elevated expression of the histone demethylase UTX. Nonetheless, the precise genomic areas impacted by these chromatin alterations remained unidentified. We used chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq) to delineate the genome-wide distribution of H3K27me3 across fat body chromatin in diapausing (D) and non-diapausing (ND) female Cx. pipiens. Notably, the higher signal at transcription start sites (TSSs) reflects localized redistribution rather than a global decrease, as diapausing fat bodies retain less H3K27me3 overall but concentrate it at promoters. To investigate the functional significance of these chromatin alterations, we confirmed a number of target loci via ChIP-qPCR and assessed gene expression with qRT-PCR. We identified many critical genes that were markedly increased in diapausing mosquitoes, exhibiting an inverse relation to H3K27me3 enrichment. Our data demonstrates different H3K27me3 chromatin landscapes between diapausing and non-diapausing Cx. pipiens, corroborating a hypothesis of selective, locus-specific repression in the non-diapause state and its targeted removal during diapause to permit activation of dormancy-associated genes. These results suggest that chromatin remodeling is a core driver of the diapause switch.

Animals↗

PubMed2026

Pervasive hybridization and introgression in Diervilleae (Caprifoliaceae).

Diervilleae (Caprifoliaceae) is a horticulturally important lineage with striking floral diversity and a long history of interspecific crossing, suggesting reticulate evolution. We integrated nuclear SNPs and whole plastome data to reconstruct a phylogenomic backbone for the tribe and to identify hybrids, cultivated accessions, and introgression among lineages. Nuclear and plastid phylogenies consistently recover Weigela and Diervilla as reciprocally monophyletic and resolve four major lineages within Weigela, providing a reproducible framework for revising sectional limits and species boundaries. Cultivated accessions form a well supported clade sister to W. florida and show predominantly W. florida ancestry while retaining contributions from multiple wild lineages, consistent with recurrent crossing, backcrossing, and selection. Analyses of wild populations reveal recurrent hybrids and enable plausible parental combinations to be inferred. Tests across the genome further indicate strong evidence for historical introgression across Diervilleae, with the strongest signals involving W. middendorffiana, W. maximowiczii, and Diervilla. Fossil evidence, divergence time estimation, and paleodistribution modelling together suggest range expansion during the Miocene and Pliocene followed by climate driven contraction, providing a spatiotemporal context for episodic contact, introgression, and the East Asia-North America disjunction.

Hybridization, Genetic↗

PubMed2026

AI-enabled viral genomics: from virus discovery to host prediction and emerging variant forecasting.

The rapid expansion of metagenomic sequencing has generated vast repositories of viral sequence data that far outpace our capacity to interpret them using conventional approaches. Highly divergent sequences, sparse functional annotation, and taxonomically uneven sampling present fundamental challenges for reference-dependent methods, which lose sensitivity precisely for novel and understudied viruses with high public health relevance. Artificial intelligence (AI) provides a new avenue to address these challenges by enabling predictive inference from viral genomes and proteins while reducing dependence on sequence similarity. In this Review, we discuss representative advances in AI for virus discovery, taxonomic classification and functional annotation, prediction of host range and zoonotic potential, and efforts toward forecasting emerging variants. These advances are transforming viral genomics from a largely descriptive discipline into one with increasing predictive capability. We also critically assess the major challenges that constrain current approaches, including the availability of high-quality and representative datasets, rigorous model evaluation, biological interpretability and responsible governance for increasingly capable AI models.

Artificial Intelligence↗

PubMed2026

Quality assessment, prognostic factors, and biomarkers for brain tumor analysis: a comprehensive systematic review.

The brain tumors possess different causative factors and properties, making their diagnosis and treatment difficult. Growth of these cancers usually leads to compression of the adjacent nerves and obstruction of the flow of cerebrospinal fluid, thus leading to increase in intracranial pressure. This affects the working of brain in many ways; thus, the difficulty involved in its treatment. With the improvements in technology in neuroimaging, including Diffusion Tensor Imaging (DTI), Positron Emission Tomography (PET), and multiparametric Magnetic Resonance Imaging (mpMRI), the diagnosis process has become easy. The effectiveness of any form of therapy in such patients depends primarily on their prognosis. While it is a common practice that physicians determine the prognosis of the disease by considering the age of the patient, histological grade of the tumor, and resection status, now this method has become more comprehensive by adding molecular signature and genetic analyses to the list of criteria. Next-generation sequencing (NGS) allows a reliable molecular classification. It increases the level of risk stratification, facilitating the application of therapies tailored to individual patients. Thus, molecular oncology has greatly changed our views on brain tumors' pathology and prognosis while neoadjuvant treatments aim at increasing the survival rate. On the other hand, radiogenomics is a field of study that combines non-invasive imaging phenotypes and genomic information in order to find unique molecular signatures of tumors without collecting samples from tumors. Molecular biomarkers are absolutely essential in the diagnosis of cancer, treatment monitoring, and recurrence of cancer. Advances in liquid biopsy technology, particularly the methods for circulating tumor DNA (ctDNA) and Extracellular Vesicle (EV) based analysis, have enabled the possibility of non-invasive monitoring of the progression of the tumors over time. This review highlights key studies and important scientific works about imaging technologies, biomarkers, and prognostic factors of malignant brain tumors.

Humans↗

PubMed2026

Boosting kynurenic acid in kombucha via substrate selection: metagenomic and biochemical insights.

Kombucha is gaining global popularity for its health benefits. This study explored the use of chestnut honey, a rich source of kynurenic acid (KYNA), to produce kombucha enriched with this metabolite. Five variants were prepared using different green/black tea blends and carbon sources: white sugar or acacia honey (controls) versus chestnut honey. Samples were analyzed for tryptophan metabolites, physicochemical properties, and microbial diversity. Komagataeibacter and Enterobacter were predominant bacterial genera in SCOBY. Candida and Aspergillus were predominated in the single sample analyzed for fungi. During fermentation, tryptophan decreased, while kynurenine increased. KYNA levels remained largely stable during fermentation and were mainly influenced by the fermentation substrate. No melatonin pathway derivatives were detected. On day 7, chestnut honey yielded kombucha with 381.680-739.915 μmol/L KYNA and elevated myricetin. Overall, chestnut honey-based kombucha represents a system in which substrate composition appears to be the main factor influencing KYNA levels in the final beverage.

Kynurenic Acid↗

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↗