Search PubMedSearch

SEARCH · Search PubMed

Results for “bacterial transcriptomics”

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.

292 records · Page 4Linked to original sources

A complete hlyCABD-like RTX operon marks a virulence-associated subset of trh-positive Vibrio parahaemolyticus from Hangzhou Bay, China.

Vibrio parahaemolyticus remains a major cause of seafood-associated gastroenteritis, yet routine surveillance still relies largely on the canonical hemolysin markers thermostable direct hemolysin (tdh) and tdh-related hemolysin (trh). To determine whether this framework overlooks accessory virulence determinants in trh-positive lineages, we analyzed 193 V. parahaemolyticus isolates collected between 2022 and 2025 from clinical, environmental, and seafood-associated sources in the Hangzhou Bay region of China. Serotyping identified 45 serotypes, with O10:K4 predominating among clinical isolates. Both clinical and non-clinical populations showed open pan-genomes, although the non-clinical group carried a larger accessory gene pool. We identified a complete hlyCABD-like RTX operon in 10 trh-positive isolates with T3SS2-associated virulence backgrounds. These RTX-positive isolates were distributed across seven sequence types and three of five phylogenetic groups. This distribution was lineage-restricted but non-clonal. In the representative hybrid-assembled genome, the operon occurred within a mosaic genomic region containing additional virulence- and mobility-associated genes, indicating a composite pathogenicity island-like element. In the tested subset, RTX-positive isolates showed significantly greater hemolytic activity than RTX-negative trh-positive isolates. This significant difference was consistently observed in both plate-based and liquid assays, and within the RTX-positive subset, hlyA expression correlated with hemolytic activity, whereas the trh gene and the tlh (thermolabile hemolysin) gene did not. A complete hlyCABD-like RTX operon therefore identifies a hemolysis-associated subset of trh-positive V. parahaemolyticus and supports its further evaluation as an additional target for food safety surveillance.

Vibrio parahaemolyticus

Genomic and Molecular Interaction Analysis of NodD1 in a Novel Bradyrhizobium yuanmingense sp. B64 Isolate for Nodulation and Symbiosis of Legume Plants.

Rhizobial bacteria are known for their ability to fix nitrogen for leguminous plants and their essential function for sustainable agriculture. This study characterizes the taxonomic status and functional potential of the Bradyrhizobium B64 isolate using integrated genomic and molecular approaches. The whole genome of the B64 isolate was sequenced via Illumina paired-end technology. Species delimitation was performed using average nucleotide identity (ANI) and digital DNA-DNA Hybridization (dDDH). The NodD1 protein structure was modeled using AlphaFold3 and validated by Ramachandran plot analysis. Molecular docking was then conducted to evaluate interactions between NodD1 and four signaling flavonoids: Apigenin, Daidzein, Genistein, and Naringenin. Genomic analysis revealed a maximum ANI of 94.4% and dDDH values between 51.4 and 62.4%. Since these values fall below the standard prokaryotic thresholds (ANI&#x2009;<&#x2009;95%; dDDH&#x2009;<&#x2009;70%), the B64 isolate is identified as a novel species. Physiological assays confirmed nitrogen fixation (1.97 ppm), IAA production (3.67 ppm), and phosphate solubilization (26.10 ppm). Structural validation showed 100% of NodD1 residues in allowed regions, ensuring high model reliability. Docking simulations demonstrated strong binding affinities across all flavonoids, with binding free energies ranging from -&#x2009;8.8 to -&#x2009;9.0&#xa0;kcal/mol. Daidzein exhibited the highest thermodynamic stability (-&#x2009;9.0&#xa0;kcal/mol), whereas apigenin showed the most extensive residue interaction network. The B64 isolate is a novel Bradyrhizobium species with a high symbiotic capacity. The stable NodD1-flavonoid interactions provide a molecular basis for efficient nodulation, positioning B64 as a promising candidate for developing lipo-chitooligosaccharide (LCO)-based biofertilizers.

Bradyrhizobium

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

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

Animals

From commensal to pathobiont: The emergence of virulence-enhanced Escherichia coli in China's food-animal systems - insights with future implications.

A fundamental shift in Escherichia coli epidemiology is being driven by convergence of virulence determinants and antimicrobial resistance within linked human-animal-environment systems. In China, the rapid growth of food-animal production, extensive antimicrobial use, and complex food networks are accelerating the emergence and dissemination of virulence-enhanced E. coli pathobionts. This review synthesizes recent epidemiological, genomics, and outbreak data to characterize China's evolving landscape of food-animal-associated E. coli. We highlight a significant shift from classical pathotypes to hybrid lineages that simultaneously carry virulence factors and last-resort antibiotic resistance determinants, including mcr-1, tet(X4), and blaNDM. These traits disseminate rapidly via plasmid-mediated horizontal gene transfer, facilitating rapid adaptation and enabling cross-sectoral One Health transmission. National surveillance, foodborne outbreak investigations, and whole-genome sequencing data show that food-animal reservoirs are active evolutionary niches that drive pathogen diversity and fitness, rather than serving merely as contamination sources. Whole-genome sequencing also pinpoints high-risk clones (e.g., ST394) and plasmid-mediated co-selection of virulence and AMR. The emergence of hybrid pathotypes (e.g., STEC/ETEC) and AMR-virulence co-selection challenges traditional classification and limits the effectiveness of conventional surveillance approaches. The 2017 colistin ban reduced mcr-1, yet ongoing resistance and emerging tet(X4) demand integrated surveillance. Collectively, these findings call for reconceptualizing E. coli as a dynamic genomic entity embedded within a unified ecological network. Addressing this threat requires an integrated One Health strategy including genomic surveillance, agricultural antimicrobial stewardship, and coordinated food-environment-clinical monitoring to prevent high-risk clone emergence and global spread.

Animals

Activity of Aztreonam-avibactam and Ceftazidime-Avibactam against Enterobacterales and Pseudomonas aeruginosa causing infections in patients hospitalized in hematology, oncology, and transplant units from United States medical centres (2019-2024).

Immunosuppression increases the risks and severity of infections and is associated with a higher incidence of infection with multidrug-resistant (MDR) pathogens. We evaluated the antimicrobial susceptibility of Enterobacterales and Pseudomonas aeruginosa from patients hospitalized in hospital units where the frequency of immunosuppressed patients is very high. Bacterial isolates were consecutively collected (1/patient) from 75 US medical centres in 2019-2024 and susceptibility tested by broth microdilution. Enterobacterales (n = 2,407) and P. aeruginosa (n = 485) from patients hospitalized in hematology, oncology, and transplant units were evaluated. Carbapenem-resistant Enterobacterales (CRE) were screened for &#x3b2;-lactamases by whole genome sequencing. Enterobacterales were mainly from bloodstream infection (BSI; 53.6%) and urinary tract infection (19.9%) and P. aeruginosa were mainly from BSI (37.9%) and pneumonia (35.0%). Aztreonam-avibactam, ceftazidime-avibactam, and meropenem-vaborbactam were highly active against Enterobacterales (99.9-99.4% susceptible), including MDR isolates (99.6-98.1% susceptible), but only aztreonam-avibactam exhibited good activity against CRE (95.8% susceptible). Ceftolozane-tazobactam showed good activity against Escherichia coli (95.7% S) and Klebsiella pneumoniae (92.8% S), but limited activity against Enterobacter cloacae species complex (75.9% susceptible). All (100.0%) carbapenemase (CBase)-producing CRE isolates were aztreonam-avibactam-susceptible while 77.4% were ceftazidime-avibactam-susceptible and 67.7% were meropenem-vaborbactam-susceptible. The most common CBases were KPC (41.7%), NDM (12.5%), and OXA-48 types (10.4%). Metallo-&#x3b2;-lactamases represented 23.5% of CBases and were identified in 16.7% of CREs. The most active agents against P. aeruginosa were ceftazidime-avibactam (95.7% susceptible), ceftolozane-tazobactam (94.8% susceptible), and tobramycin (91.5% susceptible). Piperacillin-tazobactam and meropenem were active against 81.4% and 82.5% of P. aeruginosa, respectively, and aztreonam-avibactam inhibited 78.6% of P. aeruginosa at &#x2264;8 mg/L.

Humans

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

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

Animals

Comparative genomics and full-length transcriptome profiling of wing morphs in Tetrix grossus (Orthoptera: Tetrigidae).

Wing polymorphism represents a paradigmatic dispersal-reproduction trade-off, yet its molecular basis remains uncharacterised in the phylogenetically distant pygmy grasshoppers (Tetrigidae). Here we integrate comparative genomics across ten orthopteran species with full-length transcriptomics of long-winged (FL) and short-winged (FS) Tetrix grossus. OrthoFinder recovered 118 orthogroups specific to T. grossus. Against a backdrop of pronounced gene-family contraction (36 expansions versus 222 contractions; net -186, mirrored at the ancestral Tetrix node, +37/-140), we identified an ancestral, Tetrix-specific expansion of hormone-regulation (12 genes; fold enrichment 7.93) and lipid/carbohydrate-metabolic families organised into syntenic clusters, alongside 513 positively selected genes enriched for integrin-mediated cell adhesion (6 genes), a process relevant to epithelial and appendage morphogenesis. Full-length transcriptomics of one long-winged (FL) and one short-winged (FS) adult female detected 7530 (FL) and 7515 (FS) expressed genes, with 794 FL- and 776 FS-restricted transcriptome-derived SNP-associated genes. The FL morph was enriched for an EGFR/Ras-Rho developmental-patterning axis and neuromuscular flight genes, whereas the FS morph was enriched for insulin/peptide-hormone response and growth-regulatory loci. Overall, we present genomic resources and testable hypotheses concerning the evolution and regulation of wing morphs in Tetrigidae rather than a validated genetic architecture of wing-morph determination.

Animals

Epigenetics and In Silico Transcriptome Analysis of Pediatric Acute Myeloid Leukemia.

Pediatric acute myeloid leukemia (AML) is a heterogeneous hematologic malignancy that accounts for about 15%-20% of childhood leukemias. Despite therapeutic advances, relapses remain common, and survival for high-risk patients is below 60%. Unlike adult AML, pediatric AML displays distinct genetic mutations, including FLT3-ITD, NPM1, KMT2A rearrangements, and core-binding factors (CBF) fusions, as well as extensive epigenetic dysregulation. Aberrant DNA methylation, histone modifications, and altered non-coding RNA expressions disrupt hematopoietic differentiation and activate oncogenic transcriptional networks. Recent advances in silico transcriptomic analysis have transformed the study of pediatric AML by integrating gene expression and epigenetic data to identify molecular drivers and regulatory networks. Computational RNA-seq pipelines and pathway analyses have highlighted key epigenetic regulators, including DNMT3A, TET2, and HDACs, as potential therapeutic targets. Multi-omics approaches combining transcriptomic, methylomic, and chromatin accessibility data are increasingly used to define biomarkers for diagnosis, prognosis, and therapeutic response. This review provides a comprehensive overview of the molecular and epigenetic landscape of pediatric AML, emphasizing the power of in silico transcriptome analysis to uncover disease mechanisms, refine patient stratification, and guide the development of precision-based epigenetic therapies aimed at improving long-term outcomes in children with AML.

Humans

Contribution of the Arg-33-His replacement in the histidine kinase CpxA to carbapenem and cefiderocol resistance in Serratia marcescens.

OBJECTIVE: The Serratia marcescens ROT_R clinical isolate, which was resistant to almost all &#x3b2;-lactams, including cefiderocol (4 mg/L), was recovered from a neonate 2 months after the isolation of the S. marcescens ROT_S strain that was susceptible to extended-spectrum cephalosporins (ESCs). In this study, we attempted to decipher the mechanism of resistance displayed by the ROT_R isolate. METHODS: The genomes of ROT_S and ROT_R were sequenced using the Illumina and the Oxford Nanopore Technologies. Long and short reads were assembled together, giving rise to a circularized hybrid genome. RESULTS: Genomic comparison between ROT_S and ROT_R disclosed only one mutation (G98A) in the cpxA gene of ROT_R, which led to the Arg-33-His substitution in the histidine kinase of the two-component system CpxA/CpxR. The cpxA alleles of ROT_S and ROT_R were amplified and cloned, thus giving rise to the pCpxA_WT and pCpxA_R33H recombinant plasmids, respectively, which were subsequently introduced into the S. marcescens HatR recipient strain, which lacks functional CpxA. The S. marcescens HatR (pCpxA_R33H) recombinant clone, which produced the altered CpxA_R33H variant, differed from the S. marcescens HatR (pCpxA_WT) recombinant clone, which produced the wild-type CpxA, by enhanced MICs of carbapenems and ESCs, including cefiderocol (1 mg/L). CONCLUSIONS: This study demonstrates that CpxA alteration, such as Arg-33-His substitution, can contribute to cefiderocol resistance. Although it increases slightly the MIC of cefiderocol without resulting per se in clinical resistance, it can contribute, in combination with other additional mechanisms, to achieve a high level of resistance to this siderophore cephalosporin.

Cefiderocol

A Randomized Clinical Trial to Compare Moxifloxacin Versus Azithromycin for the Treatment of Mycoplasma genitalium: The FARTHEST Study.

BACKGROUND: Mycoplasma genitalium (MG) is increasingly characterized by high rates of macrolide and fluoroquinolone resistance. International guidelines recommend resistance-guided therapy; however, access to genotypic testing is limited, and randomized trial evidence is lacking. We assessed the efficacy of moxifloxacin and azithromycin without resistance assays. METHODS: This monocentric, open-label, superiority, randomized controlled trial enrolled adults with MG infection detected by multiplex PCR, randomized 1:1 to receive moxifloxacin 400 mg daily for 10 days or azithromycin 500 mg daily for 6 days. A test of cure was performed &#x2265;28 days after treatment completion. The primary endpoint was microbiological cure in the intention-to-treat (ITT) and per-protocol (PP) populations. Subgroup analyses assessed symptomatic versus asymptomatic infections, doxycycline exposure, re-treatment, and sexual behavior. RESULTS: Among 358 randomized participants, 87.0% of those treated with moxifloxacin and 61.2% of those treated with azithromycin achieved microbiological cure in the ITT analysis (absolute risk difference 25.8%, 95% CI 16.5, 35.2). The superiority of moxifloxacin was confirmed in the ITT and PP populations. Moxifloxacin remained superior across most subgroups, whereas azithromycin showed comparable efficacy only among heterosexual individuals. Doxycycline coadministration did not improve outcomes. Both regimens were well tolerated, with only one case of discontinuation. CONCLUSIONS: Moxifloxacin demonstrated superior efficacy compared to azithromycin for treating MG infection in the absence of resistance testing. These randomized data support the use of moxifloxacin as a first-line option when resistance assays are unavailable and may inform treatment strategies.

Humans

Uncovering molecular regulatory networks of low-temperature stress response in Trachinotus ovatus via integrated transcriptome and metabolome analyses.

Golden pompano (Trachinotus ovatus) is one of the most economically important marine fish species in China. It is susceptible to low-temperature stress, which significantly challenges its production and supply. Nevertheless, study on the regulatory mechanisms underlying low-temperature stress responses in golden pompano remains limited. Here, we firstly performed a time-series transcriptome analysis to reconstruct dynamic response patterns under low-temperature stress in golden pompano. Transcriptome profiling identified common differentially expressed genes (DEGs), including fos, hlf, and hmgb1, as well as condition-specific DEGs across distinct low-temperature stress groups. Based on cluster analysis, all DEGs were classified into five distinct expression patterns, reflecting diversified regulation of expression in golden pompano during low-temperature stress. Furthermore, condition-specific regulatory modules were explored via weighted gene co-expression network analysis (WGCNA), highlighting that the two module hub genes, serbf2 and lipc, might respond to low-temperature stress by regulating the lipid catabolic process. Subsequently, untargeted metabolomic analysis revealed that glycerophospholipid metabolism was a significantly enriched common pathway, highlighting its crucial role in mediating the response to low-temperature stress. Finally, by integrating transcriptomic and metabolomic analyses, a gene-metabolite interaction network associated with glycerophospholipid metabolism under low-temperature stress was established. These findings underscore the significance of multiple candidate genes and glycerophospholipid metabolism in golden pompano's response to low-temperature stress, thereby laying a solid molecular foundation for the development of low-temperature-tolerant fish strains.

Animals

Diving Deeper Into Mechanisms of Acrylamide-Induced Toxicity: RNA Sequencing Reveals Transcriptomic Alteration and Retrotransposon Expression in Drosophila melanogaster.

Given the inevitability of human and animal exposure to acrylamide, there is increasing concern regarding its potential health risks. While a number of molecular mechanisms have been proposed, the complexity of acrylamide toxicological pathways and interactions remains incompletely characterized. In this study, we employed a transcriptomic approach to investigate the transcriptional responses of Drosophila melanogaster following exposure to acrylamide (100&#x2009;mg/kg). Our analysis identified 634 differentially expressed genes (DEGs), with 362 upregulated and 272 downregulated. Functional analysis revealed these DEGs are enriched in pathways related to reproduction, detoxification, cellular and metabolic processes, signaling, synaptic formation and organization. Notably, acrylamide exposure upregulated the expression of tau and beta-amyloid protein precursor-like genes, both implicated in Alzheimer's disease pathology. An aversive memory test further demonstrated that acrylamide impaired the short-term memory of treated flies. Additionally, acrylamide-induced toxicity altered the expression of nine long terminal repeat retrotransposons, belonging to the gypsy and pao superfamilies. By exploring the potential role of transposable element activity in acrylamide-mediated toxicity, this study provides novel insights into the molecular mechanisms underlying its effects. Collectively, these findings offer a more comprehensive understanding of the mechanisms and pathways associated with the toxic action and detoxification of acrylamide in D. melanogaster.

Animals

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

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

Adenosine Deaminase

Enrichment of Lysobacter in a long-term organically managed agricultural field with low soilborne disease incidence.

Disease-suppressive soils, in which soilborne pathogens are naturally suppressed, offer a promising model for sustainable crop protection, particularly in organic farming systems where chemical disease control options are limited. Although disease suppression in these soils is considered to rely on biological control, the underlying mechanisms remain poorly understood. In this study, we investigated soil from a long-term organically managed field in Shiga Prefecture, Japan, where soilborne disease incidence has remained consistently low, to identify bacterial community features potentially associated with this field. The 16S rRNA gene amplicon sequencing indicated that this soil harbored a bacterial community distinct from those of nearby agricultural soils. Following the application of organic compounds, the genus Lysobacter, a taxon with known antagonistic activity against plant pathogens, was markedly enriched in response to proteinaceous organic inputs. This enrichment was consistent across sampling times and specific to certain proteinaceous organic inputs, whereas minimal effects were observed on chitin, N-acetyl-d-glucosamine, or cysteine. Broader soil surveys indicated that Lysobacter enrichment was not strictly associated with whether soils had been managed under organic or conventional farming practices. Stepwise multiple regression analysis identified 10 co-occurring bacterial genera that were strongly associated with Lysobacter abundance. These findings highlight condition-dependent Lysobacter enrichment as a characteristic microbial response to proteinaceous organic amendments in this low-disease-incidence field and provide microbial insights that may inform microbiome-based strategies for sustainable soil management.

Lysobacter

Spatial transcriptomics of Ciona adult brains reveals functional zonalization and insights into neural gland function.

The ascidian Ciona is a pivotal chordate model for illuminating the evolutionary origins of the vertebrate brain. Here, spatial transcriptomics of the adult Ciona neural complex, combined with image-based computational super-resolution mapping, resolved distinct tissue domains including the cerebral ganglion, neural gland, ciliated funnel, neural gland duct/dorsal strand, and body wall muscle. Within the cerebral ganglion, high-resolution mapping revealed clear molecular zonalization separating the cortex and medulla, alongside regional specialization within the cortex itself. The neural gland exhibited localized enrichment of genes associated with extracellular matrix and cell-cell interactions. These spatial features suggest that the neural gland functions as a homeostatic and signaling interface, reminiscent of primitive vertebrate meninges or choroid plexus. Overall, this spatially defined gene expression map provides a foundational framework for understanding functional regionalization in the tunicate brain and its evolutionary relationship to vertebrate nervous systems.

Ciona

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

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

Animals