- 1College of Animal Science and Technology, Inner Mongolia Minzu University, Tongliao, China
- 2Zhalantun Vocational College, Hulunbuir, China
- 3Tongliao Agriculture and Animal Husbandry Development Center, Tongliao, China
- 4Animal Disease Prevention and Control Center of Horqin District, Tongliao, China
- 5Inner Mongolia Engineering Technology Research Center for Prevention and Control of Beef Cattle Diseases, Tongliao, China
- 6Beef Cattle Industry School of Inner Mongolia Autonomous Region, Tongliao, China
Pasteurella multocida (Pm) is a significant cause of respiratory disease in beef cattle, with serious consequences for the cattle industry. The bacterial quorum-sensing (QS) system is used for communication and regulation of various processes, including bacterial growth, virulence, biofilm formation and antimicrobial resistance. This study investigated the effects of the C6 signaling molecule on Pm resistance to the antimicrobial enrofloxacin (ENR). Bacteria were divided into groups, with one group (Pm-E1) treated with sub-inhibitory concentrations of ENR, a second treated with sub-inhibitory concentrations of ENR + 200 μM C6 (Pm-E2), and a control group (Pm-YQ). Transcriptomic, proteomic, and metabolomic sequencing of the bacteria was then performed. The results showed that C6 increased the minimum inhibitory concentration of ENR against Pm from 0.25 μg/mL to 1 μg/mL. Differentially expressed genes (DEGs), proteins (DEPs), and metabolites (DEMs), with 798 DEGs, 249 DEPs, and 499 DEMs identified in the Pm-YQ vs. Pm-E1 group, while the Pm-E1 vs. Pm-E2 group contained 784 DEGs, 301 DEPs, and 397 DEMs. Further analysis of the DEGs, DEPs, and DEMs suggested potential mechanisms by which C6 might induce Pm resistance to ENR through regulating the SOS response, CpxAR two-component system, and the ABC transporter system. These findings not only provide insight into the QS-mediated drug resistance mechanisms in Pm, but also highlight the potential of targeting the QS system for the development of novel interventions to control pasteurellosis and counteract antimicrobial resistance.
1 Introduction
Consumer demand for beef is growing rapidly in China. However, outbreaks of disease in beef cattle have serious effects on the beef industry. One such disease is bovine respiratory disease (BRD), which is common in intensive cattle farming conditions (Pardon and Buczinski, 2020; Myrenås et al., 2024). The development of BRD is associated with multiple factors, including pathogens such as viruses, bacteria, and mycoplasma, as well as predisposing conditions, such as the animal’s health and immune status, stress responses, and the husbandry environment (such as housing and management; Lachowicz-Wolak et al., 2025; Buczinski et al., 2021; Donlon et al., 2023). Among these, bacterial infections represent one of the core pathogenic components, posing significant challenges to health and leading to substantial economic losses (Monteiro et al., 2025). Pasteurella multocida (Pm), one of the primary opportunistic pathogens responsible for BRD, is frequently found in the upper respiratory tracts of healthy cattle. A subsequent decline in health, increase in stress, or poor husbandry management can induce transformation of the bacterium into a primary or secondary pathogen, leading to disease (Islam et al., 2023; E-Kobon et al., 2017). As a zoonotic pathogen, Pm can not only infect the human respiratory tract but may also act as a foodborne pathogen through contamination of beef and milk products. Consumption of undercooked beef or unpasteurized milk may lead to Pm infection via the digestive tract, although this appears to be relatively uncommon. Based on differences in its capsular antigen (K antigen), Pm can be classified into five serotypes (A, B, D, E, and F). Further subdivision according to lipopolysaccharide antigens (O antigens) yields 16 distinct serotypes (serovars 1–16) that vary in terms of host ranges and disease manifestations (Peng et al., 2019; Yu et al., 2016; Wang et al., 2023; Wang et al., 2021). Current research indicates that all five capsular serotypes can infect cattle, with serotype A primarily causing bovine pneumonia, while serotypes B and E are mainly responsible for hemorrhagic septicemia in cattle and buffalo. There are fewer reports of infections involving serotypes D and F in cattle (Wilson and Ho, 2013; Zhou et al., 2023; Elsayed et al., 2021). It is particularly noteworthy that the A:3 strain is the major serotype implicated in BRD and is often involved in co-infections with other pathogens such as Mannheimia haemolytica and Mycoplasma species, exacerbating disease severity (Yang et al., 2022).
Currently, antibiotics are the first-line drugs used for the prophylaxis and treatment of Pm infections in veterinary practice (Cuevas et al., 2020). Although antibiotics play a crucial role in controlling these infections, excessive and inappropriate use of the drugs has increased selective pressure on the expression of resistance-associated genes, leading to increased bacterial resistance and the emergence and spread of multidrug-resistant strains. The escalating problem of antimicrobial resistance has become a major concern for both veterinarians and livestock producers (Wang et al., 2025). In recent years, fluoroquinolone antibiotics, including ciprofloxacin and enrofloxacin (ENR), have been widely used in China to treat Pm infections. However, the widespread and prolonged use of these drugs has led to increased resistance to fluoroquinolone antibiotics in clinical strains (Li et al., 2023).
Quorum sensing (QS) is an important bacterial communication system. This involves the production and release of QS signaling molecules, their binding to specific receptors, and the induction of downstream signaling (Fan et al., 2022). Bacteria produce and sense signaling molecules through the QS mechanism mediated by autoinducers (AIs), enabling mutual information exchange. This type of bacterial communication is controlled by two types of QS systems, known as AI-1 and AI-2 (Pacheco et al., 2021), of which the system utilizing acyl-homoserine lactones (AHLs or AI-1) as signaling molecules in Gram-negative bacteria has been extensively studied (Coquant et al., 2020). The production and recognition of specific signaling molecules enable the coordination of bacterial gene expression and physiological activities. When the concentrations of these signaling molecules and the bacterial population density reach critical thresholds, the bacterial community initiates a series of cascade reactions to modulate the expression of specific genes, enabling bacterial adaptation to environmental changes (Zeng et al., 2023). Studies have found that QS, one of the key regulators of bacterial survival, can regulate a variety of cellular processes, including drug efflux pumps and biofilm formation (Wei and Zhao, 2018; Bäuerle et al., 2018; Wu et al., 2020). Therefore, an understanding of the bacterial QS system offers new directions for the in-depth investigation of the regulatory mechanisms underlying drug resistance, as well as the development of new drugs (Zhao et al., 2020).
In the type I QS system mediated by AHL signaling molecules, LuxI enzymes are responsible for AHL biosynthesis, with derivation of the lactone ring domain from S-adenosylmethionine (SAM). Interaction of AHL with LuxR-type receptors regulates the expression of downstream target genes. However, some bacterial strains cannot synthesize LuxI enzymes themselves but can instead express LuxR or its homologs, and are thus able to respond to exogenous AHL signaling molecules. Previous studies have confirmed that Escherichia coli, despite lacking a complete LuxR/LuxI-type QS system, can sense exogenous AHL signals through SdiA, a transcription factor belonging to the LuxR family, and can then regulate the expression of its own drug resistance-related genes (Wang et al., 2025). To date, there have been no investigations of AHL-mediated QS systems in Pm. Our preliminary studies demonstrated that exogenous supplementation of AHL signaling molecules significantly influenced drug resistance in Pm. However, analyses based on current literature and public genomic databases indicate that no typical LuxR-type receptors or their homologous proteins have been identified in Pm. This finding suggests that AHL-mediated regulation of drug resistance in Pm may not operate through the canonical QS pathway, but rather involves alternative molecular mechanisms that remain to be elucidated. As the parasitic environment of Pm contains a complex variety of strains, making it highly susceptible to interference from exogenous QS signaling molecules, the present study focused on analysis of bovine-derived Pm. We explored the effects of exogenous N-octanoyl-L-homoserine lactone (C6), a typical AHL signaling molecule, on enrofloxacin resistance in Pm. In addition, combined with transcriptomic, proteomic and metabolomic sequencing, we systematically analyzed the molecular mechanisms underlying C6 regulation of Pm resistance. The findings of the study offer a novel perspective for understanding QS-mediated antibiotic resistance in Pm, and also provide a theoretical foundation and promising interventional strategies for addressing bacterial antibiotic resistance in clinical practice.
2 Materials and methods
2.1 Bacterial strain
The capsular serotype A Pm was isolated from beef cattle with pneumonia in the Wulagai region of Inner Mongolia and is kept in the Laboratory of Preventive Veterinary Medicine, College of Animal Science and Technology, Inner Mongolia MINZU University (Wang et al., 2024a).
2.2 Drug sensitivity testing
The sensitivity of the isolated strains to ENR was determined using the microbroth dilution method recommended by Clinical and Laboratory Standards Institute (CLSI) (CLSI, 2015), with specific procedures performed as described by Wang et al. (2025). ENR (purity≥98%, Yuanye Biotechnology Co. Ltd., Shanghai, China, Product No. S17081) stock solution (2048 μg/mL) was prepared in 1 M NaOH. The experiment was performed with three biological replicates. The MIC was defined as the lowest concentration of the drug at which no visible bacterial growth was observed. The MIC values of ENR for the isolated strains were interpreted and categorized as susceptible (S), intermediate (I), or resistant (R) according to the breakpoint criteria of both CLSI for Pm and the National Antimicrobial Resistance Monitoring System (NARMS; Wang et al., 2024b).
2.3 Effect of AHL signaling molecules on drug resistance in pm
The L-type molecule N-octanoyl-L-homoserine lactone C6 was dissolved in dimethyl sulfoxide (DMSO) and physiological saline, and the effects of the solvent and C6 on the growth of the bacterial strain were evaluated. The results showed that neither the solvent nor C6 had a significant impact on bacterial growth. As described by Zhang et al., 2018, a preliminary experiment was conducted by adding the AHL C6 at concentrations of 40 μM (2 × 10−4 mol/L), 60 μM (3 × 10−4 mol/L), 80 μM (4 × 10−4 mol/L), 100 μM (5 × 10−4 mol/L), 200 μM (1 × 10−3 mol/L) and 400 μM (2 × 10−3 mol/L) to identify an appropriate concentration for use. Moreover the MIC values were identical for both the 200 μM (1 × 10−3 mol/L) and 400 μM (2 × 10−3 mol/L) concentrations. Building upon these preliminary findings and relevant literature, 200 μM (1 × 10−3 mol/L) was selected for subsequent experiments. This concentration represents the lowest concentration of C6 that showed a maximal effect on the MIC of ENR against Pm, thereby minimizing potential non-specific effects associated with higher doses. It is also the classic concentration used in previous studies on the regulation of AHL-mediated resistance in Gram-negative bacteria. The experimental procedure is described in Section 2.2. Bacterial cultures were prepared using BHI medium, and 200 μM (1 × 10−3 mol/L) of N-octanoyl-L-homoserine lactone C6 was added to the bacterial suspension containing ENR. The cultures were then incubated at 37 °C for 24 h, after which the MIC values were determined. All experiments were performed with three biological replicates.
2.4 Sample preparation
The bacteria were divided into three groups, namely, a group treated with a sub-inhibitory concentration of ENR, a second group treated with the sub-inhibitory ENR concentration plus C6, and a control group. Specifically, Group 1 (Pm-E1) was treated with 0.125 μg/mL ENR, Group 2 (Pm-E2) was treated with 0.5 μg/mL ENR with the addition of 200 μM (1 × 10−3 mol/L) C6, and the control group (Pm-YQ) was grown without either ENR or C6. All groups were incubated in BHI medium with shaking at a constant temperature of 37 °C, and the bacteria were harvested at the exponential growth phase (OD600 = 0.5).
2.5 Transcriptomic sequencing and identification of DEGs
Total RNA was extracted from the Pm-YQ, Pm-E1, and Pm-E2 groups using a Bacterial RNA Extraction Kit (Guangzhou Feiyang Bioengineering Co., Ltd., China) as specified in the manufacturer’s instructions. The RNA samples were sent to Guangzhou Gene Denovo Biotechnology Co. for high-throughput sequencing on an Illumina NovaSeq 6,000 platform (San Diego, CA, United States) and subsequent analysis. The raw reads were filtered to remove low-quality reads, including those containing adapters, nucleotides with quality scores lower than 20, and more than 10% unknown nucleotides, before mapping to the ribosomal RNA (rRNA) database using Bowtie 2 version 2.2.8. The rRNA-mapped reads were removed, and the remaining reads were mapped to the reference genome (accession number CP142011). Genes were annotated using Rockhopper, and expression levels were normalized using the fragments per kilobase of transcript per million mapped reads (FPKM) method. The edgeR package (version 3.12.1) in R1 was used for identification of the DEGs using the criteria of fold change ≥ 2 and false discovery rate (FDR) < 0.05. All DEGs were annotated using the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases.
2.6 Real-time quantitative PCR (RT-qPCR)
Eleven DEGs (Supplementary Table S1) were selected for verification of their expression levels using RT-qPCR. Each 20 μL reaction volume contained 10 μL of 2 × SYBR Green Master Mix, 0.4 μL each of the forward and reverse primers (10 μM), 1 μL of template, and 8.2 μL of ddH₂O. The conditions included an initial denaturation at 95 °C for 30 s, followed by 40 cycles of 95 °C for 10 s, and finally 60 °C for 30 s. The 16S rRNA gene was used as an internal control, and relative expression levels were calculated using the 2-∆∆Ct method.
2.7 Metabolomic sequencing and identification of DEMs
Bacterial samples prepared as described in Section 2.4 were sent to Shanghai Personalbio Biotechnology Co., Ltd. for metabolomic analysis. Chromatographic analysis was performed using an ACQUITY UPLC HSS T3 column (100 Å, 1.8 μm, 2.1 × 100 mm). An Orbitrap Exploris 120 mass spectrometer (Thermo Fisher Scientific, Waltham, MA, United States) was used to collect DDA mass spectrometry data in both positive and negative ion modes under the control of Xcalibur software (version 4.7, Thermo Fisher Scientific). Compound Discoverer 3.3 (version 3.3.2.31, Thermo Fisher Scientific) was used for peak extraction, alignment, correction, and other operations on the offline data. Peaks that were not detected in at least 50% of the QC samples were removed, while the missing values of undetected peaks were filled using Fill Gaps software, and the summed total peak area was normalized. p-values were calculated using statistical tests, the variable importance in projection (VIP) was determined by the OPLS-DA dimension reduction method, and the fold change was assessed to determine differences between two groups. Metabolites with p-values < 0.05 and VIP > 1 were considered statistically significant. Clustering analysis of the abundance of DEMs was performed using the pheatmap package (version 1.0.12) in R. Correlation analysis of DEMs was conducted using corrplot (version 4.0.3). Functional analysis of DEMs was undertaken by KEGG enrichment analysis using clusterProfiler (version 4.6.0), revealing significantly enriched metabolic pathways.
2.8 Proteomic sequencing and identification of DEPs
Bacterial samples prepared as described in Section 2.4 were sent to Shanghai Personalbio Biotechnology Co., Ltd. for protein extraction and proteomic analysis. Samples were separated using the Vanquish Neo ultra-high-performance liquid chromatography (UHPLC) system. Data-independent acquisition (DIA) analysis was undertaken with a Vanquish Neo system (Thermo Fisher Scientific) for chromatographic separation. Samples separated by nano-flow high-performance liquid chromatography underwent DIA mass spectrometry analysis using the Orbitrap Astral high-resolution mass spectrometer (Thermo Fisher Scientific). The raw mass spectrometry (MS) data were analyzed using DIA-NN (v1.8.1) with a library-free method. All data were based on 99% protein identification confidence with FDR ≤ 1%. Protein intensities were normalized and preliminarily filtered using MaxLFQ, followed by imputation of missing values using the KNN method. DEPs with fold change > 2 and p < 0.05 were identified using t-tests. Functional enrichment analyses were performed using the R package clusterProfiler (v4.10.0) for GO and KEGG analyses.
2.9 Statistical analysis
Statistical analysis was performed using GraphPad Prism 5 (Graph Software, San Diego, CA, United States). All data are expressed as the mean ± standard error based on three independent experiments. A p < 0.05 was considered statistically significant.
3 Results
3.1 Drug sensitivity results
The MIC value of ENR against Pm-YQ was determined to be 0.25 μg/mL. Treatment with 200 μM (1 × 10−3 mol/L) C6 resulted in a four-fold increase in the MIC of ENR against this strain, increasing from 0.25 to 1 μg/mL.
3.2 DEG analysis
3.2.1 Identification of DEGs
DEGs were identified in edgeR using the criteria |log2FC| > 1 and FDR < 0.05. A total of 798 DEGs were found in the Pm-YQ vs. Pm-E1 comparison group (p < 0.05), of which 294 genes were up-regulated and 504 were down-regulated. As shown in Table 1, the most significantly up-regulated DEGs were dsbA, rpoH, trmA, htpX, dnaK, clpB, and hslV, while the most significantly down-regulated DEGs were bioD, glpQ, glpK, adhE, fumC, frdC, and frdD. The Pm-E1 vs. Pm-E2 yielded 784 DEGs (p < 0.05), including 590 up-regulated and 194 down-regulated genes. As shown in Table 2, the most significantly up-regulated genes were rapA, recC, mlaC, dnaE, brnQ, fabG, and pnuC, while the most significantly down-regulated genes were mglB, galT, aspA, mglA, dhaL, rpiB, and galK. Volcano plots (Supplementary Figure S1) and heatmaps (Supplementary Figure S2) illustrate the identified DEGs.
3.2.2 GO functional annotation of DEGs
To compare the functions of the DEGs between the different groups, the up- and down-regulated DGEs identified via edgeR were analyzed using GO, with a threshold value p < 0.05. The GO categories were cellular component (CC), molecular function (MF), and biological process (BP). The DEGs in the Pm-YQ vs. Pm-E1 group were classified into 15 MFs, 3 CCs, and 24 BPs, while those in the Pm-E1 vs. Pm-E2 comparison included 15 MFs, 3 CCs, and 26 BPs (Figures 1A,B). Classification of the GO functional annotation results in Level 2 GO terms showed that most DEGs were significantly enriched in metabolic processes, cellular processes, response to stimulus, and localization in the BP category, cellular anatomical entity and protein-containing complex in the CC category, and transporter activity, catalytic activity, and binding in MFs.
Figure 1. GO functional annotation of DEGs. (A) Functional annotation of Pm-YQ vs. Pm-E1 GO DEGs. (B) Functional annotation of Pm-E1 vs. Pm-E2 GO DEGs.
3.2.3 KEGG enrichment analysis of DEGs
To investigate the pathways influencing drug resistance in Pm, the DEGs were further evaluated using the KEGG database. The results are shown in Figure 2, with higher RichFactors indicating higher degrees of enrichment. The number of DEGs in the individual pathways is indicated by the size of the dots, with larger dots representing greater numbers; dot colors indicate the q-value corrected by multiple hypothesis testing, with values ranging from 0 to 1, and the smaller the value, the more significant the enrichment. The results showed that the DEGs were primarily enriched in pathways related to metabolism, biochemical processes, and signal transduction, as well as in several specific pathways, including ATP-binding cassette (ABC) transporters, the tricarboxylic acid (TCA) cycle, two-component systems, ribosomes, RNA polymerases, and bacterial chemotaxis. To verify the reliability of the transcriptomic result, 11 DEGs (recX, recA, LexA, Lon, rpoC, dnaE, mglB, otnI, glpT, glpA, and nrdD) were randomly selected for RT-qPCR verification. Measurement of the expression levels of these 11 genes showed that their expression was consistent with the transcriptome sequencing results (Supplementary Figures S3A–D), indicating the reliability of the transcriptome sequencing results.
3.3 DEM analysis
3.3.1 Identification of DEMs
DEMs were identified using a combination of univariate and multivariate statistical analyses, using the criteria of VIP > 1 and p-value < 0.05. The Pm-YQ vs. Pm-E1 comparison identified 499 DEMs (p < 0.05), including 203 up accumulated metabolites and 296 down accumulated metabolites. The Pm-E1 vs. Pm-E2 group contained 397 DEMs (p < 0.05), of which 182 were up accumulated metabolites and 215 down accumulated metabolites. Tables 3, 4 show a selection of up accumulated and down accumulated DEMs in Pm-YQ vs. Pm-E1 and Pm-E1 vs. Pm-E2 groups, respectively. Volcano plots (Supplementary Figure S4) and heatmaps (Supplementary Figure S5) illustrate the identified DEMs.
3.3.2 KEGG enrichment analysis of DEMs
The results of the KEGG enrichment analysis revealed that the DEMs were primarily enriched in pathways involved in the synthesis, metabolism, and transport of various substances. Several of these pathways are associated with bacterial resistance, including pathways involving ABC transporters, nucleotide metabolism, purine metabolism, and pyrimidine metabolism (Figure 3). In the Pm-YQ vs. Pm-E1 group, DEMs enriched in ABC transporters and general metabolic pathways were significantly down-regulated, whereas those associated with purine metabolism were up-regulated. Conversely, in the Pm-E1 vs. Pm-E2 group, DEMs enriched in nucleotide metabolism, purine metabolism, and pyrimidine metabolism were significantly down-regulated, while up-regulation was observed for DEMs enriched in ABC transporters.
3.4 Analysis of DEPs
3.4.1 Identification of DEPs
DEPs were identified using t-tests with the criteria of fold change ≥ 1.5 or fold change ≤0.6667 and p-value < 0.05. A total of 249 DEPs (p < 0.05) were identified in the Pm-YQ vs. Pm-E1 group, of which 103 were up-regulated and 146 were down-regulated. In the Pm-E1 vs. Pm-E2 group, 301 DEPs (p < 0.05) were detected, including 152 up-regulated and 149 down-regulated proteins. Representative up- and down-regulated DEPs from the Pm-YQ vs. Pm-E1 and Pm-E1 vs. Pm-E2 comparisons are presented in Tables 5 and 6, respectively. Volcano plots (Supplementary Figure S6) and heatmaps (Supplementary Figure S7) illustrate the identified DEPs.
3.4.2 GO functional annotation and enrichment analysis of DEPs
The DEPs identified in the Pm-YQ vs. Pm-E1 group were found to be classified into 9 MFs, 2 CCs, and 11 BPs, while those in the Pm-E1 vs. Pm-E2 comparison included 11 MFs, 2 CCs, and 15 BPs (Figure 4). Classification of the GO functional annotation results in Level 2 GO terms showed that most of the DEGs were significantly enriched in metabolic processes, cellular processes, response to stimulus, and localization in the BP category, cellular anatomical entity and protein-containing complex in the CC category, and transporter activity, catalytic activity, and binding in MFs. The GO enrichment results indicated that the DEPs were significantly enriched in terms related to various damage repair processes and transport functions, including the SOS response, DNA repair, ABC transporter complex, and transmembrane transporter activity (Figure 5).
Figure 4. GO functional annotation of DEPs. (A) GO functional annotation of Pm-YQ vs. Pm-E1 DEPs. (B) GO functional annotation of Pm-E1 vs. Pm-E2 DEPs.
Figure 5. GO enrichment analysis of DEPs. (A) GO enrichment analysis of Pm-YQ vs. Pm-E1 DEPs. (B) GO enrichment analysis of Pm-E1 vs. Pm-E2 DEPs.
3.4.3 KEGG enrichment analysis of DEPs
KEGG enrichment analysis showed that in the Pm-YQ vs. Pm-E1 comparison, DEPs were enriched not only in various metabolic pathways but also in specific pathways, such as ABC transporters, nucleotide excision repair, two-component systems, and bacterial chemotaxis. In contrast, DEPs from the Pm-E1 vs. Pm-E2 comparison were primarily enriched in biosynthetic, metabolic and two-component system pathways, including carbon metabolism, pentose phosphate pathway, and microbial metabolism in different environments (Figure 6).
Figure 6. KEGG enrichment analysis of DEPs. (A) Pm-YQ vs. Pm-E1 group of DEPs. (B) Pm-E1 vs. Pm-E2 group of DEPs.
3.5 Combined transcriptomic, metabolomic, and proteomic analyses
3.5.1 Combined KEGG enrichment analysis
For better evaluation of the biological phenotypes and underlying regulatory mechanisms of the strains, the transcriptomic, metabolomic, and proteomic data were integrated. KEGG enrichment analysis revealed that, overall, the DEGs, DEMs, and DEPs were enriched in pathways involved in carbon and energy metabolism (Figure 7). Notably, the most significant enrichment was observed in pathways associated with amino sugar and nucleotide sugar metabolism, starch and sucrose metabolism, the pentose phosphate pathway, glycolysis/gluconeogenesis, the TCA cycle, and pyruvate metabolism. Furthermore, pathways associated with microbial adaptation, such as microbial metabolism in different environments, ABC transporters, and the two-component system, were also significantly enriched. The consistent enrichment patterns observed across all three omics layers indicate a coherent molecular response ranging from gene expression to metabolic flux. The heatmap (Figure 8) provides a visual summary of the differentially expressed molecules identified through the integrated analysis.
Figure 7. Combined KEGG enrichment analysis of data from the three omics analyses: (A) KEGG enrichment analysis of Pm-YQ vs. Pm-E1; (B) KEGG enrichment analysis of Pm-E1 vs. Pm-E2.
Figure 8. Combined heatmap of the results of the three omics analyses. (A) Significant DEGs, DEMs, and DEPs in Pm-YQ vs. Pm-E1. (B) Significant DEGs, DEMs, and DEPs in Pm-E1 vs. Pm-E2.
3.5.2 Correlation analysis
The transcriptomic, proteomic, and metabolomic data were integrated, enabling the construction of multi-omics association networks focusing on three key pathways, namely, ABC transporters, the two-component system, and QS, with close associations with bacterial antibiotic resistance. The aim was to examine the mechanism underlying the regulatory activity of the QS signaling molecule C6 on resistance to ENR in Pm. As shown in Figure 9, although the multi-omics data from the two groups exhibited varying degrees of association across all three pathways, the most significant interaction network was observed for the ABC transporter pathway. This suggests that bacterial strains initially perceive antibiotic stress through the two-component system, leading to the regulation of the expression of ABC transporters to induce drug efflux, ultimately resulting in the development of resistance. Notably, following exogenous supplementation of the C6 QS signaling molecule, a significant increase was observed in the number of molecules associated with the ABC transporter pathway. This further suggested that the C6 signaling molecule may specifically modulate the expression of ABC transporters, thereby promoting drug efflux activity and ultimately enhancing bacterial resistance to ENR.
Figure 9. Transcriptomic, metabolomic, and proteomic association network. Circles represent genes, triangles indicate proteins, and squares represent metabolites. (A) Association network of Pm-YQ vs. Pm-E1. (B) Association network of Pm-E1 vs. Pm-E2.
4 Discussion
Bovine respiratory diseases caused by Pm infection have had severe consequences for the recent development of the beef cattle industry in China, resulting in substantial economic losses. Fluoroquinolones are widely used as first-line antimicrobial agents in veterinary practice for treating Pm infections due to their potent bactericidal activity. However, driven by antibiotic misuse and natural selection, there has been a steady increase in Pm resistance to fluoroquinolone antibiotics, posing a serious threat to the beef cattle industry (Dhayalan et al., 2024). ENR, a fluoroquinolone antibacterial drug, interferes with the processes of DNA replication, transcription, and repair in bacteria by inhibiting the activity of DNA gyrase (topoisomerase II) and topoisomerase IV, ultimately resulting in DNA damage and death (Aboyadak and Ali, 2024; Xi, 2024). In this study, transcriptome sequencing revealed that treatment with sub-inhibitory concentrations of ENR led to marked down-regulation of genes involved in DNA synthesis, such as nrdD, consistent with the antibacterial mechanism of ENR (Torrents et al., 2000). Concurrently, the expression of genes involved in the bacterial SOS response, including recA, recX, and LexA, was significantly up-regulated. The SOS response is a critical mechanism associated with bacterial repair of DNA damage. Damage to bacterial DNA induces binding of the RecA protein to the exposed single-stranded DNA to form RecA-ssDNA nucleoprotein filaments, which mediate the autocleavage of the lexA repressor. This process derepresses the SOS response and initiates repair of the DNA damage (Mo et al., 2018; Cox, 2003; Xu et al., 2024; Cárdenas et al., 2012). The proteome sequencing showed significantly up-regulation expression of RecA, RecX, and other proteins, further verifying the activation of this pathway. Additionally, the expression of genes related to bacterial metabolism, including mglB, glpT, and glpA, was also significantly down-regulated, indicating an overall suppression of bacterial metabolic activities. The metabolome sequencing results showed general down-regulation of metabolic pathways involving various substances, corroborating the transcriptome sequencing findings. These results suggest that stress caused by sub-inhibitory concentrations of enrofloxacin induces the bacterial strain to reduce its metabolic level to conserve energy, thereby prioritizing energy allocation for the repair of DNA damage resulting from enrofloxacin exposure. In summary, we hypothesize that the bacterial strain enhances its resistance to ENR primarily through activation of the SOS response to initiate DNA damage repair, while simultaneously reducing energy consumption to allocate sufficient resources for DNA repair.
There have been numerous studies on the regulation of bacterial virulence, acid tolerance, and biofilm formation by Type I QS systems. However, to date, there has been limited investigation of their roles in the regulation of antibiotic resistance in bacteria, and almost nothing is known about their influence on resistance in Pm. Notably, the genome of Pm may lack typical LuxR-type receptors; however, studies have indicated that cells lacking LuxR homologs can still respond to AHL signaling molecules through alternative receptors. Therefore, the exogenous AHL signaling molecules supplemented in this study may not act via reliance on a complete canonical QS system, but are more likely attributable to the direct or indirect biological activities of the AHL molecules themselves (Delago et al., 2021). This study demonstrates that treatment of bacteria with the AHL signaling molecule C6 significantly modulated resistance to ENR. The results of the transcriptomic, metabolomic, and proteomic sequencing suggest that C6 may regulate bacterial resistance to ENR through three pathways.
First, C6 can modulate the bacterial SOS response. The SOS response system is known to play a critical role in the development of antibiotic resistance in bacteria. As a conserved regulatory network, it initiates appropriate responses to the presence of DNA damage, not only inducing DNA repair but also increasing the mutation frequency of target sites associated with drug resistance, thereby promoting the development of resistance (Ulrich et al., 2013; Kim et al., 2013). This study found that after the addition of sub-inhibitory concentrations of ENR, the expression of core SOS response genes, such as recA and recX, was significantly up-regulated. Furthermore, the subsequent introduction of the C6 signaling molecule led to a further up-regulation of these genes. Meanwhile, proteomic sequencing revealed significant up-regulation of proteins such as RecA and RecX, verifying the observed changes in gene expression. Additionally, the metabolomic sequencing results indicated significant down-regulation of multiple metabolic processes in the bacterial strain, resulting in reduced energy consumption. Based on these findings, it is therefore suggested that AHL signaling molecule C6 may positively regulate the expression of genes involved in DNA recombination repair such as recA, thereby activating and reinforcing the SOS repair pathway. In addition, the reduction in non-essential energy consumption provides sufficient energy for the bacterial SOS response, ultimately maintaining bacterial DNA integrity and enhancing bacterial resistance to ENR.
Second, C6 can impact the permeability of the bacterial outer membrane by activating the CpxAR two-component system. This system represents a regulatory hub involved in environmental adaptation found ubiquitously in Gram-negative bacteria, and is composed of the membrane-associated histidine kinase CpxA, the cytoplasmic response regulator CpxR, and the periplasmic accessory factor CpxP (Masi et al., 2020; Hunke et al., 2012). This system regulates various bacterial processes, including biofilm formation, environmental adaptation, antibiotic resistance, and virulence factor expression (Gahlot et al., 2022; Hirakawa et al., 2005; He et al., 2022). This study, through multi-omics sequencing, revealed that the addition of AHL signaling molecule C6 led to significant enrichment of genes, metabolites, and proteins associated with the two-component system in the bacterial strain, most notably up-regulation of the cpxA gene, indicating that AHL signaling molecule C6 may activate the CpxAR system either directly or indirectly. Additionally, significant differential expression was observed in outer membrane proteins and their associated regulatory genes, as seen in the marked down-regulation of the ompA gene. Studies have shown that the outer membrane protein OmpA not only participates in bacterium-host interactions, biofilm formation, and serum resistance but also represents an important channel enabling penetration of the antibiotic into the bacterium (Scribano et al., 2024). Furthermore, Hu et al. found that the CpxAR system can regulate Salmonella Typhimurium resistance to ceftriaxone by influencing the expression of the outer membrane proteins OmpD and STM1530 (Hu et al., 2011). Together, these findings suggest that the AHL signaling molecule C6 may enhance bacterial resistance to ENR by activating the CpxAR two-component system, which in turn modulates the expression of genes encoding outer membrane proteins, reduces the permeability of the outer membrane, and decreases the rate and total amount of antibiotic entry.
Third, C6 can contribute to the regulation of drug efflux mediated by ABC transporter systems. The ABC superfamily, one of the largest families of transmembrane proteins in bacteria, plays a significant role in the active efflux of antibiotics. When antibiotics penetrate the outer membrane and enter the cytoplasm, they may be recognized and bound directly by specific ABC transporters, utilizing energy derived from ATP hydrolysis to actively transport the drugs from the cell (Thomas and Tampé, 2020). The present study found that the addition of the C6 signaling molecule to bacteria led to significant up-regulation of multiple genes, metabolites, and proteins within the ABC transporter pathway, with notable increases in the expression of the bacA gene encoding the YddA transporter. YddA was initially proposed as a putative drug efflux transporter in E. coli. Subsequent research by Feng et al. confirmed that YddA is a multidrug efflux transporter belonging to the ABC family that relies on ATP hydrolysis for energy, and can recognize and transport a variety of structurally and functionally unrelated substrates, including quinolone drugs, and exhibits typical multidrug efflux characteristics. Moreover, that study observed markedly increased transcription of the YddA-associated gene when quinolone drug concentrations were raised (Feng et al., 2020). Combined with the results of this study, these findings suggest that the AHL signaling molecule C6 can enhance the active efflux of ENR by Pm by regulating the expression of the ABC family YddA drug efflux transporter, thereby mediating Pm resistance to ENR.
Based on multi-omics sequencing analysis, this study provides preliminary evidence that C6 may influence bacterial resistance to ENR through various pathways, including activation of the SOS response, modulation of the CpxAR two-component system mediating outer membrane permeability, and regulation of ABC transporter-mediated drug efflux. These findings not only elucidate the QS-mediated resistance mechanism in Pm but also provide a theoretical basis and suggest potential targets for the development of novel QSIs. However, it should be noted that the integrated multi-omics approach employed in this study is fundamentally correlative in nature, and its conclusions are subject to certain limitations. Multi-omics techniques typically provide static molecular snapshots at specific time points, making it difficult to fully capture the inherent dynamic complexity involved in temporal changes in QS regulation. The interpretation of the data also relied on the completeness and accuracy of existing annotation databases, which may have affected the reliability of the pathway inferences. Furthermore, while multi-omics data can effectively delineate global molecular trends and generate important hypotheses, they cannot provide direct confirmation of specific regulatory mechanisms. For example, the multiple drug resistance mechanisms hypothesized in this study still require functional validation through targeted follow-up experiments such as gene knockout, complementation assays, or protein function interference. Future studies should investigate the specific molecular mechanisms by which C6 regulates drug resistance in Pm to provide a theoretical foundation for the development of new resistance inhibitors and antimicrobial agents, as well as opening new avenues for the prevention and treatment of bovine pasteurellosis.
5 Conclusion
This study found that the addition of 200 μM (1 × 10−3 mol/L) AHL signaling molecule C6 to the test strains significantly enhanced their resistance to ENR. Based on multi-omics analyses, it is hypothesized that the AHL signaling molecule may regulate the resistance of the test strains to ENR through three potential pathways: First, it modulates the bacterial SOS response by upregulating the expression of DNA recombination and repair genes such as recA, thereby activating and enhancing the SOS repair pathway to protect bacterial DNA integrity and consequently strengthen bacterial resistance to ENR. Second, it activates the CpxAR two-component system to modulate outer membrane permeability, reducing both the rate and total amount of antibiotic entry into bacterial cells, thus increasing bacterial resistance to ENR. Third, it is involved in regulating the ABC transporter-mediated drug efflux system by modulating the expression of the ABC family efflux transporter YddA, enhancing the active efflux capacity of Pm against ENR, thereby mediating Pm resistance to ENR.
Data availability statement
All sequences have been deposited in the NCBI Sequence Read Archive, the MetaboLights and the ProteomeXchange database under accession numbers PRJNA1343251, PXD069645 and MTBLS13502, respectively, and are publicly accessible via the following links: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1343251, https://proteomecentral.proteomexchange.org/cgi/GetDataset?ID=PXD069645 and https://www.ebi.ac.uk/metabolights/editor/MTBLS13502/descriptors.
Author contributions
MS: Conceptualization, Writing – original draft. ZZ: Data curation, Writing – original draft. DH: Conceptualization, Writing – original draft. XG: Conceptualization, Writing – original draft. YW: Conceptualization, Writing – original draft. JW: Conceptualization, Writing – original draft. YL: Conceptualization, Writing – original draft. HC: Conceptualization, Writing – original draft. BX: Conceptualization, Writing – original draft. KA: Conceptualization, Writing – original draft. DY: Conceptualization, Funding acquisition, Project administration, Writing – original draft. KL: Conceptualization, Funding acquisition, Project administration, Supervision, Writing – original draft, Writing – review & editing. ZW: Conceptualization, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This study was supported by the Inner Mongolia Autonomous Region Central Guiding Local Development Fund Project (2022ZY0149, Inner Mongolia, China), the Inner Mongolia Autonomous Region Natural Science Foundation Project (2025MS03111 and 2023QN03010, Inner Mongolia, China).
Acknowledgments
We would like to thank the Laboratory of Pharmacology and Toxicology, School of Animal Medicine, Jilin Agricultural University for help with the experiments.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that Generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2026.1766173/full#supplementary-material
Footnotes
References
Aboyadak, I., and Ali, N. G. (2024). Enrofloxacin, effective treatment of Pseudomonas aeruginosa and Enterococcus faecalis infection in Oreochromis niloticus. Microorganisms. 12:901. doi: 10.3390/microorganisms12050901,
Bäuerle, T., Fischer, A., Speck, T., and Bechinger, C. (2018). Self-organization of active particles by quorum sensing rules. Nat. Commun. 9:3232. doi: 10.1038/s41467-018-05675-7,
Buczinski, S., Achard, D., and Timsit, E. (2021). Effects of calfhood respiratory disease on health and performance of dairy cattle: a systematic review and meta-analysis. J. Dairy Sci. 104, 8214–8227. doi: 10.3168/jds.2020-19941,
Cárdenas, P. P., Carrasco, B., Defeu Soufo, C., César, C. E., Herr, K., Kaufenstein, M., et al. (2012). RecX facilitates homologous recombination by modulating RecA activities. PLoS Genet. 8:e1003126. doi: 10.1371/journal.pgen.1003126,
CLSI (2015). Performance standards for antimicrobial disk and dilution susceptibility tests for Bacteria isolated from animals; clinical laboratory standard institution[CLSI]. Clinical and Laboratory Standards Institute: Wayne, PA.
Coquant, G., Grill, J. P., and Seksik, P. (2020). Impact of N-acyl-homoserine lactones, quorum sensing molecules, on gut immunity. Front. Immunol. 11:1827. doi: 10.3389/fimmu.2020.01827,
Cox, M. M. (2003). The bacterial RecA protein as a motor protein. Ann. Rev. Microbiol. 57, 551–577. doi: 10.1146/annurev.micro.57.030502.090953,
Cuevas, I., Carbonero, A., Cano, D., García-Bocanegra, I., Amaro, M. Á., and Borge, C. (2020). Antimicrobial resistance of Pasteurella multocida type B isolates associated with acute septicemia in pigs and cattle in Spain. BMC Vet. Res. 16:222. doi: 10.1186/s12917-020-02442-z,
Delago, A., Gregor, R., Dubinsky, L., Dandela, R., Hendler, A., Krief, P., et al. (2021). A bacterial quorum sensing molecule elicits a general stress response in Saccharomyces cerevisiae. Front. Microbiol. 12:632658. doi: 10.3389/fmicb.2021.632658,
Dhayalan, A., Prajapati, A., Yogisharadhya, R., Chanda, M. M., and Shivachandra, S. B. (2024). Anti-quorum sensing and anti-biofilm activities of Pasteurella multocida strains. Microb. Pathog. 197:107085. doi: 10.1016/j.micpath.2024.107085,
Donlon, J. D., McAloon, C. G., Hyde, R., Aly, S., Pardon, B., and Mee, J. F. (2023). A systematic review of the relationship between housing environmental factors and bovine respiratory disease in preweaned calves - part 1: Ammonia, air microbial count, particulate matter and endotoxins. Vet. J. 300-302:106031. doi: 10.1016/j.tvjl.2023.106031,
E-Kobon, T., Leeanan, R., Pannoi, S., Anuntasomboon, P., Thongkamkoon, P., and Thamchaipenet, A. (2017). OmpA protein sequence-based typing and virulence-associated gene profiles of Pasteurella multocida isolates associated with bovine haemorrhagic septicaemia and porcine pneumonic pasteurellosis in Thailand. BMC Vet. Res. 13:243. doi: 10.1186/s12917-017-1157-6,
Elsayed, M. S. A. E., Eldsouky, S. M., Roshdy, T., Said, L., Thabet, N., Allam, T., et al. (2021). Virulence determinants and antimicrobial profiles of Pasteurella multocida isolated from cattle and humans in Egypt. Antibiotics 10:480. doi: 10.3390/antibiotics10050480,
Fan, Q., Zuo, J., Wang, H., Grenier, D., Yi, L., and Wang, Y. (2022). Contribution of quorum sensing to virulence and antibiotic resistance in zoonotic bacteria. Biotechnol. Adv. 59:107965. doi: 10.1016/j.biotechadv.2022.107965,
Feng, Z., Liu, D., Liu, Z., Liang, Y., Wang, Y., Liu, Q., et al. (2020). Cloning and functional characterization of putative Escherichia coli ABC multidrug efflux transporter YddA. J. Microbiol. Biotechnol. 30, 982–995. doi: 10.4014/jmb.2003.03003,
Gahlot, D. K., Wai, S. N., Erickson, D. L., and Francis, M. S. (2022). Cpx-signalling facilitates Hms-dependent biofilm formation by Yersinia pseudotuberculosis. NPJ Biofilms Microbiomes. 8:13. doi: 10.1038/s41522-022-00281-4,
He, X., Zhang, J., Wang, S., Yang, Z., Zhang, H., and Zhou, X. (2022). Cloning, expression, purification, and biochemical characterization of CpxR protein from pectobacterium carotovorum. Biotechnol. Appl. Biochem. 69, 898–905. doi: 10.1002/bab.2161,
Hirakawa, H., Inazumi, Y., Masaki, T., Hirata, T., and Yamaguchi, A. (2005). Indole induces the expression of multidrug exporter genes in Escherichia coli. Mol. Microbiol. 55, 1113–1126. doi: 10.1111/j.1365-2958.2004.04449.x,
Hu, W. S., Chen, H. W., Zhang, R. Y., Huang, C. Y., and Shen, C. F. (2011). The expression levels of outer membrane proteins STM1530 and OmpD, which are influenced by the CpxAR and BaeSR two-component systems, play important roles in the ceftriaxone resistance of Salmonella enterica serovar typhimurium. Antimicrob. Agents Chemother. 55, 3829–3837. doi: 10.1128/AAC.00216-11,
Hunke, S., Keller, R., and Müller, V. S. (2012). Signal integration by the Cpx-envelope stress system. FEMS Microbiol. Lett. 326, 12–22. doi: 10.1111/j.1574-6968.2011.02436.x,
Islam, E. A., Fegan, J. E., Tefera, T. A., Curran, D. M., Waeckerlin, R. C., Ng, D., et al. (2023). Reverse vaccinology-based identification of a novel surface lipoprotein that is an effective vaccine antigen against bovine infections caused by Pasteurella multocida. PLoS Pathog. 19:e1011249. doi: 10.1371/journal.ppat.1011249,
Kim, J., Noh, J., and Park, W. (2013). Insight into norfloxacin resistance of Acinetobacter oleivorans DR1: target gene mutation, persister, and RNA-Seq analyses. J. Microbiol. Biotechnol. 23, 1293–1303. doi: 10.4014/jmb.1307.07059,
Lachowicz-Wolak, A., Chmielina, A., Przychodniak, I., Karwańska, M., Siedlecka, M., Klimowicz-Bodys, M., et al. (2025). Antimicrobial-resistance and virulence-associated genes of Pasteurella multocida and Mannheimia haemolytica isolated from polish dairy calves with symptoms of bovine respiratory disease. Microorganisms 13:491. doi: 10.3390/microorganisms13030491,
Li, X. S., Qi, Y., Xue, J. Z., Xu, G. Y., Xu, Y. X., Li, X. Y., et al. (2023). Transcriptomic changes and satP gene function analysis in Pasteurella multocida with different levels of resistance to Enrofloxacin. Vet Sci. 10:257. doi: 10.3390/vetsci10040257,
Masi, M., Pinet, E., and Pagès, J. M. (2020). Complex response of the CpxAR two-component system to β-lactams on antibiotic resistance and envelope homeostasis in Enterobacteriaceae. Antimicrob. Agents Chemother. 64, e00291–e00220. doi: 10.1128/AAC.00291-20,
Mo, C. Y., Culyba, M. J., Selwood, T., Kubiak, J. M., Hostetler, Z. M., Jurewicz, A. J., et al. (2018). Inhibitors of LexA autoproteolysis and the bacterial SOS response discovered by an academic-industry partnership. ACS Infect. Dis. 4, 349–359. doi: 10.1021/acsinfecdis.7b00122,
Monteiro, H. F., Hoyos-Jaramillo, A., Garzon, A., Fritz, H. M., Miramontes, C. C., Matos, I. M., et al. (2025). Antibiogram use on dairy cattle for bovine respiratory disease: factors associated with bacterial pathogen identification and prediction of bacterial recovery using machine learning. J. Dairy Sci. 108, 6284–6296. doi: 10.3168/jds.2024-26205,
Myrenås, M., Pringle, M., Harbom, B., and Bengtsson, B. (2024). Pasteurella multocida from deep nasal swabs and tracheobronchial lavage in bovine calves from Sweden. Acta Vet. Scand. 66:58. doi: 10.1186/s13028-024-00781-7,
Pacheco, T., Gomes, A. É. I., Siqueira, N. M. G., Assoni, L., Darrieux, M., Venter, H., et al. (2021). SdiA, a quorum-sensing regulator, suppresses fimbriae expression, biofilm formation, and quorum-sensing signaling molecules production in Klebsiella pneumoniae. Front. Microbiol. 12:597735. doi: 10.3389/fmicb.2021.597735,
Pardon, B., and Buczinski, S. (2020). Bovine respiratory disease diagnosis: what Progress has been made in infectious diagnosis? Vet. Clin. North Am. Food Anim. Pract. 36, 425–444. doi: 10.1016/j.cvfa.2020.03.005,
Peng, Z., Wang, X., Zhou, R., Chen, H., Wilson, B. A., and Wu, B. (2019). Pasteurella multocida: genotypes and genomics. Microbiol. Mol. Biol. Rev. 83, e00014–e00019. doi: 10.1128/MMBR.00014-19,
Scribano, D., Cheri, E., Pompilio, A., Di Bonaventura, G., Belli, M., Cristina, M., et al. (2024). Acinetobacter baumannii OmpA-like porins: functional characterization of bacterial physiology, antibiotic-resistance, and virulence. Commun. Biol. 7:948. doi: 10.1038/s42003-024-06645-0,
Thomas, C., and Tampé, R. (2020). Structural and mechanistic principles of ABC transporters. Annu. Rev. Biochem. 89, 605–636. doi: 10.1146/annurev-biochem-011520-105201,
Torrents, E., Buist, G., Liu, A., Eliasson, R., Kok, J., Gibert, I., et al. (2000). The anaerobic (class III) ribonucleotide reductase from Lactococcus lactis. Catalytic properties and allosteric regulation of the pure enzyme system. J. Biol. Chem. 275, 2463–2471. doi: 10.1074/jbc.275.4.2463,
Ulrich, R. L., Deshazer, D., Kenny, T. A., Ulrich, M. P., Moravusova, A., Opperman, T., et al. (2013). Characterization of the Burkholderia thailandensis SOS response by using whole-transcriptome shotgun sequencing. Appl. Environ. Microbiol. 79, 5830–5843. doi: 10.1128/AEM.00538-13,
Wang, Z., Shi, J.-c., Wang, Y.-q., Sun, M., Meng, L.-h., Geng, C., et al. (2024a). Whole genome sequencing and genome evolution analysis of capsular serotype a and D Pasteurella multocida of bovine[J]. Biotechnol. Bull. 40, 282–290. doi: 10.13560/j.cnki.biotech.bull.1985.2024-0371
Wang, Z., Sun, M., Guo, S., Wang, Y., Meng, L., Shi, J., et al. (2024b). Detection of drug resistance in Escherichia coli from calves with diarrhea in the Tongliao region: an analysis of multidrug-resistant strains. Front Vet Sci. 11:1466690. doi: 10.3389/fvets.2024.1466690,
Wang, Z., Sun, M., Wang, Y., Shi, J., Gao, W., Han, D., et al. (2025). Regulation of ofloxacin resistance in Escherichia coli strains causing calf diarrhea by quorum-sensing acyl-homoserine lactone signaling molecules. Front Vet Sci. 12:1540132. doi: 10.3389/fvets.2025.1540132,
Wang, F., Wang, X., Ai, W., Zeng, D., Liang, W., Hua, L., et al. (2021). Three novel immunogenic proteins determined through 2-dimensional electrophoresis and mass spectrometry with immune serum confer protection against challenge with porcine Pasteurella multocida in mouse models. Res. Vet. Sci. 136, 303–309. doi: 10.1016/j.rvsc.2021.03.013,
Wang, H., Xin, L., Wu, Y., Liu, Y., Yao, W., Zhang, H., et al. (2023). Construction of a one-step multiplex real-time PCR assay for the detection of serogroups a, B, and E of Pasteurella multocida associated with bovine pasteurellosis. Front Vet Sci. 10:1193162. doi: 10.3389/fvets.2023.1193162,
Wei, C., and Zhao, X. (2018). Induction of viable but Nonculturable Escherichia coli O157:H7 by low temperature and its resuscitation. Front. Microbiol. 9:2728. doi: 10.3389/fmicb.2018.02728,
Wilson, B. A., and Ho, M. (2013). Pasteurella multocida: from zoonosis to cellular microbiology. Clin. Microbiol. Rev. 26, 631–655. doi: 10.1128/CMR.00024-13,
Wu, S., Liu, J., Liu, C., Yang, A., and Qiao, J. (2020). Quorum sensing for population-level control of bacteria and potential therapeutic applications. Cell. Mol. Life Sci. 77, 1319–1343. doi: 10.1007/s00018-019-03326-8,
Xi, F. (2024). The enrofloxacin pollution control from fish to environment. Mar. Pollut. Bull. 199:115923. doi: 10.1016/j.marpolbul.2023.115923
Xu, G., Li, P., Xue, Z., Qi, Y., Li, X., Zhu, D., et al. (2024). RecA inhibitor epicatechin prolongs the development of fluoroquinolone resistance in Pasteurella multocida. Int. J. Biol. Macromol. 255:128026. doi: 10.1016/j.ijbiomac.2023.128026,
Yang, W., Li, M., Zhang, C., Zhang, X., Guo, M., and Wu, Y. (2022). Pathogenicity, colonization, and innate immune response to Pasteurella multocida in rabbits. BMC Vet. Res. 18:416. doi: 10.1186/s12917-022-03517-9,
Yu, C., Sizhu, S., Luo, Q., Xu, X., Fu, L., and Zhang, A. (2016). Genome sequencing of a virulent avian Pasteurella multocida strain GX-pm reveals the candidate genes involved in the pathogenesis. Res. Vet. Sci. 105, 23–27. doi: 10.1016/j.rvsc.2016.01.013. 27033902,
Zeng, X., Zou, Y., Zheng, J., Qiu, S., Liu, L., and Wei, C. (2023). Quorum sensing-mediated microbial interactions: mechanisms, applications, challenges and perspectives. Microbiol. Res. 273:127414. doi: 10.1016/j.micres.2023.127414,
Zhang, Y., Ma, Q., Su, B., Chen, R., Lin, J., Lin, Z., et al. (2018). A study on the role that quorum sensing play in antibiotic-resistant plasmid conjugative transfer in Escherichia coli. Ecotoxicology 27, 209–216. doi: 10.1007/s10646-017-1886-0,
Zhao, X., Yu, Z., and Ding, T. (2020). Quorum-sensing regulation of antimicrobial resistance in bacteria. Microorganisms 8:425. doi: 10.3390/microorganisms8030425,
Keywords: drug resistance, metabolome sequencing, Pasteurella multocida, proteome sequencing, quorum-sensing molecule, transcriptomic sequencing
Citation: Sun M, Zhang Z, Han D, Ge X, Wang Y, Wang J, Li Y, Chen H, Xu B, Ao K, Yang D, Liu K and Wang Z (2026) Regulation of drug resistance to enrofloxacin in Pasteurella multocida strains from cattle by quorum-sensing acyl-homoserine lactone signaling molecules. Front. Microbiol. 17:1766173. doi: 10.3389/fmicb.2026.1766173
Edited by:
Kumaragurubaran Karthik, Tamil Nadu Veterinary and Animal Sciences University, IndiaReviewed by:
Carlos Gabriel Nieto Peñalver, National Scientific and Technical Research Council (CONICET), ArgentinaPande Gde Sasmita Julyantoro, Udayana University, Indonesia
Copyright © 2026 Sun, Zhang, Han, Ge, Wang, Wang, Li, Chen, Xu, Ao, Yang, Liu and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kai Liu, bGl1a2FpNzIxMDI2QDE2My5jb20=; Zi Wang, d2FuZ3ppeGlhb3h1QDE2My5jb20=
†These authors have contributed equally to this work and share first authorship
Zipeng Zhang1†