Skip to main content

Mapping the widespread distribution and transmission dynamics of linezolid resistance in humans, animals, and the environment

Abstract

Background

The rise of linezolid resistance has been widely observed both in clinical and non-clinical settings. However, there were still data gaps regarding the comprehensive prevalence and interconnections of linezolid resistance genes across various niches.

Results

We screened for potential linezolid resistance gene reservoirs in the intestines of both humans and animals, in meat samples, as well as in water sources. A total of 796 bacteria strains out of 1538 non-duplicated samples were identified to be positive for at least one linezolid resistance gene, optrA, poxtA, cfr, and cfr(D). The prevalence of optrA reached 100% (95% CI 96.3–100%) in the intestines of pigs, followed by fish, ducks, and chicken at 77.5% (95% CI 67.2–85.3%), 62.0% (95% CI 52.2–70.9%), and 61.0% (95% CI 51.2–70.0%), respectively. The meat and water samples presented prevalences of 80.0% (95% CI 70.6–87.0%) and 38.0% (95% CI 25.9–51.9%), respectively. The unreported prevalence of the cfr(D) gene was also relatively higher at 13.0% (95% CI 7.8–21.0%) and 19.0% (95% CI 10.9–25.6%) for the feces samples of ducks and pigs, respectively. Enterococci were the predominant hosts for all genes, while several non-enterococcal species were also identified. Phylogenetic analysis revealed a significant genetic distance among linezolid resistance gene reservoirs, with polyclonal structures observed in strains within the same niche. Similar genetic arrays harboring assorted insertion sequences or transposons were shared by reservoirs displaying heterogeneous backgrounds, though large diversity in the genetic environment of linezolid resistance genes was also observed.

Conclusions

The linezolid resistance genes were widespread among various niches. The horizontal transfer played a crucial role in driving the circulation of linezolid resistance reservoirs at the human-animal-environment interfaces.

Video Abstract

Background

The increasing occurrence of antimicrobial resistance (AMR) is currently recognized as a global public health crisis. This issue expands beyond healthcare facilities, as interconnected human, animal, and environmental habitats collectively contribute to the emergence, evolution, and spread of AMR [1]. Consequently, a coordinated and multisectoral approach, such as One Health, has become imperative to comprehend the role of interconnected ecosystems in the complex dynamics of AMR and effectively address the resistance problem [2, 3]. From this perspective, tracking the circulation of AMR genes at the human-animal-environment interface holds profound significance and is considered one of the pivotal issues in One Health [4], particularly concerning the genes that bestow resistance to clinically important antibiotics, such as colistin, carbapenem, vancomycin, and oxazolidinones [5].

Linezolid, the pioneering member of the oxazolidinone family, stands as one of the few agents that remain effective against multidrug-resistant Gram-positive bacteria, including methicillin-resistant Staphylococcus aureus (MRSA) and vancomycin-resistant enterococci (VRE) [6]. However, in recent years, the emergence of linezolid resistance in Gram-positive communities worldwide has significantly compromised the clinical utility of this antimicrobial agent. Resistance to linezolid is commonly mediated by mutations in domain V of 23 s rRNA [7], modifications in ribosomal proteins L3 and L4 [8], or acquisition of transferable resistance determinants [9]. These determinants include optrA, poxtA, cfr, as well as its variants. The cfr gene [10] and cfr-like genes, specifically designated as cfr(B) [11], cfr(C) [12], cfr(D) [13], and cfr(E) [14], confer cross-resistance to oxazolidinones, phenicols, lincosamides, pleuromutilins, and streptogramin A (PhLOPSA phenotype) by inhibiting bacterial protein biosynthesis [15]. To date, members of the cfr family have been detected in at least two genera of Gram-positive bacteria, except for cfr(C) and cfr(E), which have been exclusively described in Clostridium and/or Campylobacter [9].

The novel plasmid-borne optrA and poxtA genes encode ATP-binding cassette (ABC)-F proteins that confer either resistance or decreased susceptibility to oxazolidinones and phenicols through a ribosomal protection mechanism [16, 17]. Since the identification of the optrA gene in enterococci of both human and animal origin in 2015 [16], it has subsequently been detected in various bacteria, including Staphylococcus, Streptococcus, Clostridium, Campylobacter, Aerococcus, Lactococcus, and Vagococcus [9, 18,19,20,21], exhibiting a remarkable diversity of nucleotide sequences [9, 22, 23]. Initially described in a clinical MRSA in 2018 [17], the poxtA gene has also been identified in Enterococcus and Ligilactobacillus [24]. Recently, a variant of the poxtA gene, poxtA2, has been discovered to be co-location with the cfr(D) gene in a conjugative plasmid of enterococci [25].

Previous studies have demonstrated that linezolid resistance genes are embedded into a plethora of variable genetic environments, involving multiple mobile genetic elements (MGEs). These elements play a crucial role in facilitating the horizontal transfer and further dissemination of these genes across strains, species, genus boundaries, or various potential niches [21, 23, 26,27,28,29]. Although linezolid has thus far been exclusively approved for human use, the cfr, optrA, and poxtA genes have been reported to drive the proliferation of linezolid resistance in Gram-positive bacteria within healthcare settings (22, 30,31,32,33,34), among animals and their products [35,36,37,38,39,40,41,42,43,44,45], as well as in various environmental contexts [9, 46,47,48,49,50,51]. This underscores the significance of monitoring these linezolid resistance determinants through the One Health approach.

Despite concerted efforts to gain a better understanding of the prevalence of linezolid resistance in numerous countries, existing studies generally concentrate on single or sporadic One Health sectors or the occurrence in specific genera, primarily enterococci. Considering the promiscuous nature of MGEs and conjugative plasmids, the rapid dissemination of the resistance locus is highly plausible, and the prevalence is still likely to be underestimated. In 2022, we reported an optimized screening approach with higher sensitivity, resulting in a higher prevalence of the optrA gene in fecal samples among healthy individuals at 19.3% in Hangzhou, China [52]. However, there are still data gaps concerning the epidemiological situation of these genes in Gram-positive communities from various niches and the relationships between these linezolid resistance reservoirs. In this study, we aimed to gain comprehensive knowledge of epidemiologic and genetic levels for linezolid resisnce among different niches and depict the potential routes of the transmissions through the One Health approach.

Methods

Sample collection and bacterial isolation

A total of 1538 samples were collected from Hangzhou City in China, spanning the period from February to December 2022. The samples encompassed various sources, including feces from healthy individuals (n = 1018), feces from farmed animals (100 each for chicken, duck, and pig), fish samples (n = 80), pork and poultry meat obtained from the supermarkets (n = 90), as well as water samples from five rivers and lakes in Hangzhou City (n = 50) (refer to Supplementary Table S1 for details).

As part of our research, we isolated 245 strains from 1018 human feces in our previous study [52]. The fecal samples for pigs, ducks, and chickens were collected using cloacal swabs. The fecal samples for fish were obtained in the form of intestinal content samples. Non-duplicated samples were collected. For meat samples, a standardized protocol was followed, wherein 50 g of meat from each sample was thoroughly vortexed with 40 mL of 0.1 M phosphate buffer. Subsequently, the mixture was homogenized to ensure accurate representation. In the case of water samples, a meticulous procedure was employed, wherein 500 mL of water was filtered through 0.45 µm filter membranes (Pall, USA). These filter membranes were then subjected to three rounds of washing with 10 mL of sterile saline, ensuring optimal sample preparation. The acquired liquid products were then combined and mixed. The isolation of Gram-positive cocci carrying linezolid resistance genes was conducted using our optimized approach, as previously described [52]: pretreated products and feces samples were inoculated into 5 ml of Luria–Bertani (LB) broth and incubated at 37 °C for 24 h. Subsequently, 100 μl of each enrichment was transferred to a subculture of 5 ml fresh LB broth containing 5% NaCl and 10 mg/L florfenicol. The Columbia agar base media, supplemented with 5% (v/v) sheep blood and 10 mg/L of florfenicol, was used to select and purify the putative target isolates.

Species identification and PCR analysis

Species identification was performed using MALDI-TOF MS (Bruker Daltonik GmbH, Germany) with the Biotyper Version 3.0 (DB2969 database). Additionally, confirmation of species identification was obtained by comparing the isolates with reference strains through the online tool Average Nucleotide Identity (ANI) calculator (http://enve-omics.ce.gatech.edu/ani/index) for those isolates that underwent WGS. Subsequently, isolates determined to be of the same species according to DB2969 database (including Enterococcus dongliensis, E. hulanensis, E. lactis, E. viikkiensis, and E. xiangfangensis) were employed to construct an in-house database. This in-house database was then utilized for the re-confirmation of non-sequenced enterococci, ensuring accuracy and reliability in the species identification process.

The presence of optrA, poxtA, cfr, and cfr-like genes was ascertained through PCR analysis and Sanger sequencing. To identify OptrA variants, a comparison was made between the deduced amino acid sequences of the isolates and the original OptrA from E. faecalis E349, which had been previously designated as the wild type [16].

Antimicrobial susceptibility testing

The minimal inhibitory concentrations (MICs) of seven antimicrobial agents, namely linezolid, chloramphenicol, penicillin G, ciprofloxacin, erythromycin, tetracycline, and vancomycin, were determined using the broth microdilution method [53]. The interpretation of results was conducted in accordance with the standard of Clinical and Laboratory Standards Institute (CLSI) [54, 55]. For the Vagococcus and Globicatella genera, the susceptibility breakpoints of the seven antimicrobial agents, as established for enterococci, were applied. To ensure the accuracy and reliability of the analysis, E. faecalis ATCC 29212, S. aureus ATCC 29213, and Streptococcus pneumoniae ATCC 49619 were employed as the quality control strains.

Whole-genome sequencing and genome analysis

The cohort selection of isolates subjected to WGS was carried out through a systematic approach, involving stratification based on their sources, species, and carriage of linezolid resistance genes. All isolates identified as non-enterococcal species underwent sequencing. Additionally, all enterococci carrying any of the poxtA, cfr, or cfr-like genes were also subjected to WGS. To manage the considerable number of enterococci exclusively carrying optrA, 10% of strains from each species within each niche were randomly selected using a Microsoft Excel randomization program for WGS. This approach resulted in a final cohort of 294 isolates for WGS, comprising 208 enterococci and 86 non-enterococcal ones (Supplementary Table S1).

The genomic DNA underwent WGS utilizing the Illumina NovaSeq 6000 platform. Subsequently, the sequencing data were de novo assembled using SPAdes Version 3.13.1 [56]. The identification of antimicrobial resistance genes (ARGs) and the sequence types (STs) for the assembly scaffolds were performed with default settings, utilizing ResFinder Version 4.1 [57] and MLST Version 2.0 [58], respectively. The newly obtained STs were uploaded and assigned through PubMLST (https://pubmlst.org/). Furthermore, the 16S rRNA genes of all isolates were obtained using Barrnap Version 0.9 (https://github.com/tseemann/barrnap). The alignment of the 16S sequences was conducted using MAFFT Version 7 [59]. For the construction of the maximum likelihood phylogenetic tree based on the 16S rRNA of all strains subjected to WGS, IQ-TREE Version 1.6.12 was employed [60]. To facilitate core-genome alignment, single-nucleotide polymorphism (SNP) calling, and the construction of SNP-based phylogenetic trees, the Parsnp script [61] was utilized. Subsequently, all the phylogenetic trees were visualized presented and annotated using the online tool iTOL Version 3 [62]. BLASTN analysis was performed to compare and annotate the contigs containing the target genes with known sequences of the NCBI database. The sequences of strains with short linezolid resistance gene-carrying contigs were mapped to the known structures in this study, and the potential gaps between contigs were filled by PCR analysis and Sanger sequencing to generate complete structures. Additionally, the linear alignment of the genetic environment of the linezolid resistance genes in different isolates was visualized using EasyFig Version 2.2.2 [63].

Results

Prevalence and distribution of linezolid resistance genes in various niches

The prevalence and distribution of linezolid resistance genes were illustrated in Table 1. A total of 796 florfenicol-resistant strains were obtained, comprising 245 isolates of human origin, 78 isolates from meat, 74 isolates from cloacal swabs of chicken, 102 isolates from duck feces, 191 isolates from pig feces, and 24 isolates from water. The optrA gene was most frequently detected among each niche, with an overall prevalence of 37.2% (95% CI 34.8–39.6%). The highest prevalence of optrA was exhibited in the intestine of pigs (100.0%, 95% CI 96.3–100%), followed by meat samples (80.0%, 95% CI 70.6–87.0%). Among fish samples, a prevalence of 77.5% (95% CI 67.2–85.3%) was observed, which was about twice that for water samples at 38.0% (95% CI 25.9–51.9%). The fish showed the highest prevalence of the poxtA gene at 15.0% (95% CI 8.8–24.4%), followed by pigs at 12.0% (95% CI 7.0–19.8%). In water, the prevalence was 8.0% (95% CI 3.2–18.8%). The cfr gene was mainly detected in isolates from pig feces at 19.0% (95% CI 12.5–27.8%). Likewise, a higher prevalence of the cfr(D) gene in the intestine of pigs was observed at 17.0% (95% CI 10.9–25.6%), and an intestinal carriage rate of ducks was also observed at 13.0% (95% CI 7.8–21.0%). The cfr(B), cfr(C), or cfr(E) genes were not identified among these isolates.

Table 1 Prevalence and distribution of linezolid resistance genes in various niches

Enterococci were the primary hosts for the optrA gene in most of the niches, except for fish. A total of 689 optrA-positive enterococci of 19 species were collected from various niches (Supplementary Table S1), with E. faecalis (n = 402), E. faecium (n = 91), and E. avium (n = 64) being the most frequently identified species. In contrast, within the intestine of the fish, Lactococcus lactis (35/81) was found to be the more prevalent carrier of the optrA gene, followed by E. avium (17/81) and L. petauri (11/81). Additionaly, a variety of optrA-carrying non-enterococcal isolates were obtained from diverse niches, including Lactococcus (L. raffinolactis, L. garvieae, and L. formosensis), Vagococcus (Vagococcus lutrae, V. fluvialis, and V. carniphilus), Streptococcus (Streptococcus gallolyticus and S. parauberis), Aerococcus urinaeequi, and Globicatella sulfidifaciens. Notably, this study marks the first detection of the optrA gene in L. lactis, L. raffinolactis, V. fluvialis, V. carniphilus, S. parauberis, A. urinaeequi, and G. sulfidifaciens.

Regarding the poxtA gene, it was predominantly detected in enterococci, while one L. lactis from pig feces and one L. raffinolactis of fish origin also tested positive for this gene. The cfr gene was only present among enterococci. However, the cfr(D) gene was also detected in non-enterococcal isolates including V. carniphilus of porcine origin, V. fluvialis from chicken, and S. parauberis from ducks. All cfr- or cfr(D)-positive isolates showed concomitant carriage of the optrA and/or poxtA gene(s), with the exception of a single strain where these genes were absent.

Antimicrobial susceptibility results

A similar antimicrobial susceptibility pattern with high percentages of resistance to chloramphenicol, erythromycin, and tetracycline was observed across different genera of isolates (Supplementary Table S1). Among the 796 isolates, the overall percentage of resistance to linezolid was 62.8% (Fig. 1d), with a similar resistance rate of 67.3% observed for enterococci (Fig. 1a). The resistance rates for lactococci and vagococci were relatively lower at 28.8% and 23.8%, respectively (Fig. 1b, c). However, it is worth noting that 59.3% of lactococci exhibited decreased susceptibility to linezolid, leaving only 11.9% of them remained susceptible to this antibiotic. About half of the strains showed resistance to ciprofloxacin, but the rates were relatively higher for lactococci and vagococci, reaching 86.4% and 76.2%, respectively. All vagococci remained susceptible to penicillin G, and this biotic also showed high efficiency against 78.6% of the enterococci obtained. Notably, with the exception of intrinsically resistant strains, no vancomycin-resistant isolates were observed in this study.

Fig. 1
figure 1

Antimicrobial susceptibility profiles for linezolid resistance gene reservoirs. LZD, linezolid; C, chloramphenicol; CIP, ciprofloxacin; P, penicillin G; E, erythromycin; TE, tetracycline; R, resistance rate; I, intermediate rate; S, susceptibility rate

Identification of ARG and OptrA variants

Whole-genome sequencing (WGS) analysis was conducted to comprehensively characterize the genetic features of the selected cohort, comprising 208 enterococci and 86 nonenterococcal isolates. The detailed analysis revealed the presence of five poxtA-positive E. hulanensis strains (D5-2, D6-4, D19-2, D36-2, and D86-2), one E. faecalis strain (M6), and one E. raffinosus strain (B646-2) harboring the poxtA2 variant, which was recently reported in the literature [25]. Furthermore, a total of 36 ARGs were characterized (Supplementary Table S1). These genes collectively confer resistance to phenicols, macrolides, fosfomycin, tetracyclines, trimethoprim, and aminoglycosides, thus aligning with the multiple drug resistance phenotypes previously observed for these isolates.

The alignment of the amino sequence of OptrAE349, designated as the wild-type (WT), revealed substitutions in 34 positions of the amino acid sequence (Supplementary Figure S1), with Y176D, K3E, and G393D substitutions being commonly identified. These alterations formed a total of 34 OptrA variants among 217 optrA-positive isolates subjected to WGS. The most frequently identified variant was the EDM type (n = 43), followed by WT (n = 32), DD (n = 25), and EYDNDM (n = 25). A total of nine novel OptrA variants were designated as ED'DM, EYDDNI, EYDNDNM, EYRCDVKVDAMINI, FK, IDKTKGP, KLDKK, KLDKP, and V (Fig. 2). Corresponding MIC values for each OptrA variant fell within different ranges (Supplementary Figure S1). Higher linezolid MIC values were commonly associated with the WT, EDP, KD, KLDP, and RDK variants, regardless of the genera or species. No distinct genus boundary for the OptrA variants was found, as isolates from different species or sources shared the same or similar OptrA variants. Moreover, in nine strains, binary carriage of two types of OptrA variants within the same strain was observed, including two E. faecalis strains (M88 and P42-1), six L. lactis strains (M55, Y75, Y7, Y39, Y73, and Y74-2), and one E. avium strain (M44-2) (Supplementary Table S1).

Fig. 2
figure 2

The flow of the source, species, OptrA variants, and basic structure of genetic context of optrA-carrying isolates

Phylogenetic analysis of linezolid resistance gene-carrying isolates

To investigate the phylogeny of linezolid resistance reservoirs, multiple alignments of the 16S rRNA sequences of these strains were performed to construct the maximum likelihood tree (Fig. 3). The analysis revealed that the 294 strains from seven genera clustered in different branches, and even strains identified as the same genus displayed significant phylogenetic distance. Moreover, no apparent aggregation of strains from a single source or with the same linezolid resistance gene patterns was observed.

Fig. 3
figure 3

Phylogeny of strains carrying linezolid resistance genes. The phylogenetic tree was generated based on the 16S rRNA sequences of linezolid resistance gene reservoirs by the maximum likelihood method. The genera and source of each strain were illustrated as the legend indicated. The tiny squares with linear arrangements from the inside to the outside indicated the carriage of linezolid resistance genes

To further investigate the genetic backgrounds and relatedness of linezolid resistance gene-carrying isolates among different niches, SNPs for frequently isolated species were analyzed and mapped. A total of 13 SNP-based maximum-likelihood phylogenetic trees were generated for Enterococcus [E. faecalis (n = 53), E. faecium (n = 36), E. avium (n = 34), E. dongliensis (n = 15), E. gallinarum (n = 13), E. casseliflavus (n = 9), E. raffinosus (n = 9), E. hirae (n = 7), and E. hulanensis (n = 6)], Lactococcus [L. lactis (n = 38) and L. petauri (n = 15)], and Vagococcus [V. lutrae (n = 11) and V. carniphilus (n = 7)] (Supplementary Figure S2). Generally, large phylogenetic differences were observed for both enterococci and non-enterococci, even strains originating from the same source were dispersed across the different nodes on the tree, indicating the variety of genetic backgrounds of these linezolid resistance reservoirs. However, despite there being obvious clonal aggregation yet to be found, some clades still contained closely related isolates from the same or various sources. A distinct clustering profile of the optrA-positive L. lactis population in the intestine of fish was mapped (Supplementary Figure S2c), but six strains with different profiles of OptrA variants still showed close relations (Clade A). A similar scenario was also discovered in V. carniphilus from pigs (Supplementary Figure S2l). Three isolates of Clade A that displayed one amino acid alteration in OptrA (EYDNDM VS. EYDNDNM) were clonally related with a co-carriage of the cfr(D) gene. Similarly, it was of note that the inter-niche clonal spread was always associated with the binary carriage of linezolid resistance genes in our study. Three E. faecium isolates derived from water samples of two conjoint rivers were clustered together (Clade A) with the coexisting of the optrADD and poxtA genes (Supplementary Figure S2b). For E. avium, five highly related isolates from duck feces displayed the co-occurrence of the optrAEDM and cfr(D) genes (Clade A, Supplementary Figure S2d). The clonal dissemination of binary linezolid resistance gene carriers was also observed in Clade A of E. raffinosus from fish samples that were positive for optrAEDDM and poxtA genes (Supplementary Figure S2j), as well as E. hulanensis of duck origin that carried poxtA2 and cfr(D) (Supplementary Figure S2m). Not only the single clonal spread was observed, but two individual clades of L. petauri of fish origin were clustered as Clade A and Clade B, the isolates also showed a close genetic distance inner each clade (Supplementary Figure S2f).

The potential epidemiology linkages were also noticed across niches boundary. Two E. avium strains of human origin and four strains of porcine origin were clustered closely together on Clade B of the phylogenetic tree (Supplementary Figure S2d). All six of these strains showed the binary carriage of both optrAEDM and cfr(D) genes. In the case of E. faecalis, ST16 was the most frequently identified sequence type, displaying a polyclonal structure in the phylogenetic tree (Supplementary Figure S2a). However, this ST type was still shared by strains from various sources, including human feces (n = 2), duck feces (n = 3), pig feces (n = 1), pork and poultry meat (n = 3), and water samples (n = 1). Following that, the ST116 (n = 4) was also shared by the strains derived from the feces of humans (n = 1) and pigs (n = 1) as well as water samples (n = 2), indicating ST16 and ST116 types of E. faecalis may act as the vehicle for the transmission of linezolid resistance among different niches. Despite the linkage for the inter-niche clonal spread was not as evident as that within a single niche, closely related strains from various sources still demonstrated the phylogenomic relationships and potential transmission routes of these linezolid resistance reservoirs.

Genetic environment of linezolid resistance genes

The arrays for the genetic environment of the optrA, poxtA, poxtA2, cfr, and cfr(D) genes were represented as illustrations in Supplementary Figure S3. A high variation of optrA-carrying genetic platforms with abundance and diversity in MGEs and ARGs was observed, which yielded plenty of different arrangements of genetic environment. IS1216E (n = 41) was mostly identified upstream and/or downstream of the optrA gene as previous studies extensively described [9]. Further, one or multiple copies of assorted other ISs including ISVlu1 (n = 17), ISEnfa1 (n = 5), ISEnfa5 (n = 4), IS-LL6 (n = 3), IS256 (n = 1), IS1251 (n = 1), and ISS1N (n = 1) were also found to flank the optrA-core regions. Transposons were found in the vicinity of optrA as well, including the complete Tn554 and Tn558 or their relics in different sizes. However, the sequence at the junctions of the transposons varied with respect to the target sites, indicating the different origin and integration routes by individual mobilization events. The co-location of other ARGs of these optrA-carrying contigs was characterized, with fexA (n = 62) and erm(A) (n = 53) being most commonly detected, further contributing to the co-selection of optrA.

To better clarify the corresponding among these genetic environments with the OptrA variants, species of hosts, and sources, the 105 arrangements were simplified to 22 basic types as Fig. 2 listed. The arrangement of Tn558-araC-optrA (n = 84) was most prevalent among the optrA-positive isolates. Following that, the arrangements of IS1216-fexA-optrA-IS1216 (n = 48), and IS1216- optrA-IS1216 (n = 25) were also common. The same or similar arrangements were shared by different optrA variants, suggesting the possible interconnection of these strains and a lack of conservation of optrA during dissemination.

The genetic environment for the poxtA gene was relatively consistent (Supplementary Figure S3j), while IS1216E was found to flank these contigs in the same or opposite orientations. The recently reported poxtA2 gene, whose 3′ end was not disrupted by IS1216E but directly associated with the cfr (D) was discovered in seven enterococci of different species and sources, and the arrangement of these contigs was homologous to that identified in porcine manure [25]. Consistent with that reported elsewhere, the cfr(D) gene was commonly associated with a complete or truncated guaA gene downstream and a truncated ISSeq2 upstream, flanked or not by IS1216E (Supplementary Figure S3j). However, it was of note that the insertion of optrAWT into the common arrangement of cfr(D) thus generated the array of IS1216E-ΔISSeq2-cfr(D)-optrA-IS1216E (Supplementary Figure S3h). It was observed in S. parauberis D84-2, an unreported reservoir for these two linezolid resistance genes. The ISEnfa5-cfr-ISEnfa5 genetic array was most frequently identified in the cfr-carrying enterococci (n = 16), followed by that of IS256-cfr-IS256 (n = 6), both of which were dissimilar with that characterized in staphylococci in our earlier study [30]. Overall, though showed diversity, these genetic contexts of linezolid resistance genes were shared by the heterogeneous reservoirs in the levels of species and sources, indicating the horizontal transfer and potential interconnections of these reservoirs.

Discussion

Representative studies on surveillance of the epidemiological situation of linezolid resistance via the One Health approach in various niches have been conducted on a global scale. However, comprehensive and multifield research on the resistance to this last-line resort antibiotic is still rare, and its prevalence has been underestimated. This might be attributed to the urgent need to optimize sensitive and standardized methodologies for the detection of ARGs. Our previous study reported a human intestinal carriage rate of optrA at 19.3% with an optimized high-sensitivity screening approach, which displayed a more than tripled rate compared to the traditional method currently widely used [52]. We used this method to screen for samples from multiple niches to comprehensively perceive the epidemiologic phenomenon of all linezolid resistance genes, beyond only focusing on clinical settings. As a result, we discovered an unexpectedly high prevalence and a large diversity of linezolid resistance reservoirs at multiple levels.

The wide distribution of linezolid resistance genes across various ecological niches has been emphasized, and the remarkably high prevalence of the optrA gene has raised significant concern. In addition, our study provided comprehensive knowledge on the prevalence of large-scale samples of the cfr(D) gene, which was also relatively high among these resistance genes. The animal-origin samples were identified as important reservoirs for linezolid resistance in our study, suggesting that these resistance determinants might have accumulated in the intestinal microbiota of animals. Despite linezolid not being approved for use in the livestock or poultry industries, our findings demonstrated that the distribution of linezolid resistance genes could be affected by the residual concentration of florfenicol in livestock manures [64]. This indicates that the extensive use of veterinary antibiotics may lead to the enrichment and dissemination of acquired linezolid resistance genes in the gut microbiota of livestock through co-selection. The high prevalence of linezolid resistance genes among farm animals is a cause for concern, as these resistance determinants may enter the food chain via various routes, such as spillover at the slaughterhouse or contamination of the circulated water as a form of excretion [65]. This could potentially result in human intestinal colonization and therapeutic failure, which poses a potential risk to public health [66].

Our study also unveiled the remarkable diversity in linezolid resistance, exemplified by the heterogeneous backgrounds of the resistance reservoirs. Firstly, we made a valuable addition to the database of host species carrying linezolid resistance genes. Though enterococci showed a predominance in harboring the four linezolid resistance genes, we revealed the previously neglected role of multiple non-enterococcal species, indicating the widespread distribution of these genes in Gram-positive communities. Notably, the previously unnoticed L. lactis emerged as the most abundant species carriering the optrA gene in the intestinal of fish. Additionally, optrA and poxtA were observed in other Lactococcus spp., suggesting the relevance of linezolid resistance dissemination in the ecosystem for these bacteria as commensal intestinal flora of humans and animals [67]. The Vagococcus spp. also behaved as an important linezolid reservoir, which was equally represented in a recent study in Europe [44]. The occurrence of linezolid resistance genes in less frequently studied species indicates that the reservoirs we have currently discovered likely represent only a fraction of the larger picture, and many uncharacterized individuals with widespread phenomena remain to be identified.

Even within Enterococcus spp., which is considered one of the best-studied Gram-positive genera, we observed a diverse array of species carrying linezolid resistance genes, including less identified enterococci such as E. hulanensis, E. viikkiensis, and E. xiangfangensis, among others. The carriage of linezolid resistance genes in these infrequent enterococci is noteworthy, particularly considering that enterococci have been proposed for monitoring antibiotic resistance in food animals [65]. These less-common enterococci are easily overlooked during surveillance, making their identification crucial for comprehensive monitoring efforts. Further, not only did we observe species diversity, but strains of the same species displayed substantial phylogenetic distances or possessed different linezolid resistance genes. The polyclonal structure of each species population further underscores the highly diverse phylogenetic lineages of the embedded linezolid resistance genes, consistent with findings from previous studies [34, 41, 68].

Another dimension of diversity in linezolid resistance was manifested in the highly variable genetic platforms of these resistance genes. We discovered a wealth of MGEs in the flanking regions of linezolid resistance genes, potentially playing a crucial role in the horizontal transfer of these resistance determinants. Among those elements, the significance of IS1216 was emphasized, as it was most frequently identified and has been shown to facilitate the mobilization of optrA- and poxtA- carrying contigs in the form of IS1216-based translocatable units [23, 26]. Furthermore, these linezolid resistance loci were also commonly associated with the co-selection of multiple additional ARGs, potentially promoting the persistence of these linezolid resistance genes within the Gram-positive population. The abundance and diversity of MGEs and ARGs are indicative of genetic structural plasticity, which is also evident within the linezolid resistance genes themselves, as reflected in the presence of multiple gene variants.

The greatest genetic flexibility was observed in the highly variable amino acid sequence of OptrA, a distinctive feature of the optrA gene contributing to the diverse linezolid resistance levels we observed. However, the carriage of OptrA profiles did not consistently correspond to specific resistance phenotypes, suggesting the involvement of potential additional mechanisms.

Despite the extensive diversity observed, there were still certain similarities and potential relatedness between these linezolid resistance reservoirs. Shared gene variants and genetic contexts among heterogeneous reservoirs of linezolid resistance, underscored the flow of linezolid resistance genes in the ecosystem. The absence of represented epidemic clones we observed suggested that these linezolid resistance determinants could be readily acquired or exchanged through horizontal transfer, facilitated by the abundance of various insertion sequences and transposons that enable mobilization and transposition events. Clonal spread also contributed to the dissemination of linezolid resistance, primarily within specific niches. Although its significance was less pronounced across diverse niches, some strains of different origins still displayed limited phylogenetic distance, implying potential epidemiological linkages among multiple niches.

One of the limitations of this study pertains to the sample collection, which was conducted in a single-center manner, and as such, the epidemic features observed only represent the situation in Hangzhou City. To gain a more comprehensive understanding of the circulation of these linezolid resistance genes, further work and investigations on a larger geographical scale are warranted. As optrA, poxtA, cfr, and cfr-like genes confer a wide range of MIC values to linezolid (varied from 2 to 32 mg/L) while higher values to florfenicol at least for 16 mg/L; thus, we chose florfenicol to isolate the target isolates. The unexpectedly high prevalence reflected by this study methodology may showed slight discrepancy but was much closer to the true prevalence.

Conclusions

In conclusion, our study has provided comprehensive knowledge of linezolid resistance at the human-animal-environment interface and highlighted a remarkably high prevalence of the optrA gene in multiple niches. Additionally, we have reported the previously unexplored prevalence of poxtA, cfr, and cfr(D) genes, collectively demonstrating that animal intestines and associated settings serve as important linezolid resistance gene reservoirs. These resistance determinants have embedded into a wide range of genetic backgrounds and represent a prevalent situation among Gram-positive communities, not only in enterococci as the predominant host, but also in a variety of non-enterococcal species, including Lactococcus, Vagococcus, Streptococcus, Aerococcus, and Globicatella. The horizontal transfer facilitated by MGEs in the flanking region of linezolid resistance genes has been identified as the primary mechanism responsible for the wide dissemination of linezolid resistance on a large scale. Additionally, the clonal spread has also contributed to the spread of resistance. Our study provides evidence of potential clonal relatedness and epidemiology linkage among heterogenous niches, it is important to note that further research is needed to fully understand the circulation of linezolid resistance genes in the ecosystem. Nevertheless, the data obtained in this study can serve as a valuable baseline for assessing the capacity of linezolid resistance genes in multiple niches, demonstrating the power of One Health surveillance in informing policies aimed at improving public health.

Availability of data and materials

The whole-genome sequencing data set was retrieved with BioProject accession numbers PRJNA948548 and PRJNA848114.

Abbreviations

AMR:

Antimicrobial resistance

MRSA:

Methicillin-resistant Staphylococcus aureus

VRE:

Vancomycin-resistant enterococci

ABC:

ATP-binding cassette

MGE:

Mobile genetic element

ANI:

Average Nucleotide Identity

MIC:

Minimal inhibitory concentrations

CLSI:

Clinical and Laboratory Standards Institute

WGS:

Whole-genome sequencing

ARG:

Antimicrobial resistance gene

ST:

Sequence type

SNP:

Single-nucleotide polymorphism

WT:

Wild-type

References

  1. Hernando-Amado S, Coque TM, Baquero F, Martínez JL. Defining and combating antibiotic resistance from One Health and Global Health perspectives. Nat Microbiol. 2019;4:1432–42. https://doi.org/10.1038/s41564-019-0503-9.

    Article  CAS  PubMed  Google Scholar 

  2. Iriti M, Vitalini S, Varoni EM. Humans, animals, food and environment: One health approach against global antimicrobial resistance. Antibiotics. 2020;9:346. https://doi.org/10.3390/antibiotics9060346.

    Article  PubMed  PubMed Central  Google Scholar 

  3. McEwen SA, Collignon PJ. Antimicrobial resistance: a One Health perspective. Aarestrup FM, Schwarz S, Shen J, Cavaco L, editors. Microbiol Spectr. 2018;6:6 2 10. https://doi.org/10.1128/microbiolspec.ARBA-0009-2017.

    Article  Google Scholar 

  4. Van Puyvelde S, Deborggraeve S, Jacobs J. Why the antibiotic resistance crisis requires a One Health approach. Lancet Infect Dis. 2018;18:132–4. https://doi.org/10.1016/S1473-3099(17)30704-1.

    Article  PubMed  Google Scholar 

  5. Koutsoumanis K, Allende A, Álvarez-Ordóñez A, Bolton D, Bover-Cid S, Chemaly M, et al. Role played by the environment in the emergence and spread of antimicrobial resistance (AMR) through the food chain. EFSA J. 2021;19:e06651. https://doi.org/10.2903/j.efsa.2021.6651.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Meka VG, Gold HS. Antimicrobial Resistance to Linezolid. Clin Infect Dis. 2004;39:1010–5. https://doi.org/10.1086/423841.

    Article  CAS  PubMed  Google Scholar 

  7. Tsiodras S, Gold HS, Sakoulas G, Eliopoulos GM, Wennersten C, Venkataraman L, et al. Linezolid resistance in a clinical isolate of Staphylococcus aureus. Lancet. 2001;358:207–8. https://doi.org/10.1016/S0140-6736(01)05410-1.

    Article  CAS  PubMed  Google Scholar 

  8. Mendes RE, Deshpande LM, Farrell DJ, Spanu T, Fadda G, Jones RN. Assessment of linezolid resistance mechanisms among Staphylococcus epidermidis causing bacteraemia in Rome. Italy J Antimicrob Chemother. 2010;65:2329–35. https://doi.org/10.1093/jac/dkq331.

    Article  CAS  PubMed  Google Scholar 

  9. Schwarz S, Zhang W, Du X-D, Krüger H, Feßler AT, Ma S, et al. Mobile oxazolidinone resistance genes in Gram-positive and Gram-negative bacteria. Clin Microbiol Rev. 2021;34:e00188-e220. https://doi.org/10.1128/CMR.00188-20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Schwarz S, Werckenthin C, Kehrenberg C. Identification of a plasmid-borne chloramphenicol-florfenicol resistance gene in Staphylococcus sciuri. Antimicrob Agents Chemother. 2000;44:2530–3. https://doi.org/10.1128/AAC.44.9.2530-2533.2000.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Deshpande LM, Ashcraft DS, Kahn HP, Pankey G, Jones RN, Farrell DJ, et al. Detection of a New cfr-Like Gene, cfr(B), in Enterococcus faecium isolates recovered from human specimens in the United States as part of the SENTRY antimicrobial surveillance program. Antimicrob Agents Chemother. 2015;59:6256–61. https://doi.org/10.1128/AAC.01473-15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Tang Y, Dai L, Sahin O, Wu Z, Liu M, Zhang Q. Emergence of a plasmid-borne multidrug resistance gene cfr(C) in foodborne pathogen Campylobacter. J Antimicrob Chemother. 2017;72:1581–8. https://doi.org/10.1093/jac/dkx023.

    Article  CAS  PubMed  Google Scholar 

  13. Pang S, Boan P, Lee T, Gangatharan S, Tan SJ, Daley D, et al. Linezolid-resistant ST872 Enteroccocus faecium harbouring optrA and cfr(D) oxazolidinone resistance genes. Int J Antimicrob Agents. 2020;55:105831. https://doi.org/10.1016/j.ijantimicag.2019.10.012.

    Article  CAS  PubMed  Google Scholar 

  14. Stojković V, Ulate MF, Hidalgo-Villeda F, Aguilar E, Monge-Cascante C, Pizarro-Guajardo M, et al. cfr (B), cfr (C), and a New cfr-Like Gene, cfr(E), in Clostridium difficile strains recovered across Latin America. Antimicrob Agents Chemother. 2019;64:e01074-e1119. https://doi.org/10.1128/AAC.01074-19.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Long KS, Poehlsgaard J, Kehrenberg C, Schwarz S, Vester B. The Cfr rRNA methyltransferase confers resistance to phenicols, lincosamides, oxazolidinones, pleuromutilins, and streptogramin A antibiotics. Antimicrob Agents Chemother. 2006;50:2500–5. https://doi.org/10.1128/AAC.00131-06.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Wang Y, Lv Y, Cai J, Schwarz S, Cui L, Hu Z, et al. A novel gene, optrA, that confers transferable resistance to oxazolidinones and phenicols and its presence in Enterococcus faecalis and Enterococcus faecium of human and animal origin. J Antimicrob Chemother. 2015;70:2182–90. https://doi.org/10.1093/jac/dkv116.

    Article  CAS  PubMed  Google Scholar 

  17. Antonelli A, D’Andrea MM, Brenciani A, Galeotti CL, Morroni G, Pollini S, et al. Characterization of poxtA, a novel phenicol–oxazolidinone–tetracycline resistance gene from an MRSA of clinical origin. J Antimicrob Chemother. 2018;73:1763–9. https://doi.org/10.1093/jac/dky088.

    Article  CAS  PubMed  Google Scholar 

  18. Cai J, Chen J, Schwarz S, Wang Y, Zhang R. Detection of the plasmid-borne oxazolidinone/phenicol resistance gene optrA in Lactococcus garvieae isolated from faecal samples. Clin Microbiol Infect. 2021;27:1358–9. https://doi.org/10.1016/j.cmi.2021.04.027.

    Article  CAS  PubMed  Google Scholar 

  19. Zhu Y, Yang W, Schwarz S, Xu Q, Yang Q, Wang L, et al. Characterization of the novel optrA-carrying pseudo-compound transposon Tn7363 and an Inc18 plasmid carrying cfr(D) in Vagococcus lutrae. J Antimicrob Chemother. 2022;77:921–5. https://doi.org/10.1093/jac/dkab478.

    Article  CAS  PubMed  Google Scholar 

  20. Coccitto SN, Fioriti S, Cinthi M, Morroni G, Di Giannatale E, Mingoia M, et al. Detection of phenicol-oxazolidinone resistance gene optrA in Aerococcus viridans from bovine faeces. Italy J Antimicrob Chemother. 2021;76:2479–81. https://doi.org/10.1093/jac/dkab162.

    Article  CAS  PubMed  Google Scholar 

  21. Zhou Y, Li J, Schwarz S, Zhang S, Tao J, Fan R, et al. Mobile oxazolidinone/phenicol resistance gene optrA in chicken Clostridium perfringens. J Antimicrob Chemother. 2020;75:3067–9. https://doi.org/10.1093/jac/dkaa236.

    Article  CAS  PubMed  Google Scholar 

  22. Cai J, Schwarz S, Chi D, Wang Z, Zhang R, Wang Y. Faecal carriage of optrA-positive enterococci in asymptomatic healthy humans in Hangzhou. China Clin Microbiol Infect. 2019;25(630):e1-630.e6. https://doi.org/10.1016/j.cmi.2018.07.025.

    Article  CAS  Google Scholar 

  23. Brenciani A, Morroni G, Schwarz S, Giovanetti E. Oxazolidinones: mechanisms of resistance and mobile genetic elements involved. J Antimicrob Chemother. 2022;77:2596–621. https://doi.org/10.1093/jac/dkac263.

    Article  CAS  PubMed  Google Scholar 

  24. Zhu Y, Yang W, Schwarz S, Xu Q, Yang Q, Wang L, et al. Characterization of an MDR Lactobacillus salivarius isolate harbouring the phenicol-oxazolidinone-tetracycline resistance gene poxtA. J Antimicrob Chemother. 2022; dkac169. doi:https://doi.org/10.1093/jac/dkac169.

  25. Cinthi M, Coccitto SN, Fioriti S, Morroni G, Simoni S, Vignaroli C, et al. Occurrence of a plasmid co-carrying cfr(D) and poxtA2 linezolid resistance genes in Enterococcus faecalis and Enterococcus casseliflavus from porcine manure. Italy J Antimicrob Chemother. 2022;77:598–603. https://doi.org/10.1093/jac/dkab456.

    Article  CAS  PubMed  Google Scholar 

  26. He T, Shen Y, Schwarz S, Cai J, Lv Y, Li J, et al. Genetic environment of the transferable oxazolidinone/phenicol resistance gene optrA in Enterococcus faecalis isolates of human and animal origin. J Antimicrob Chemother. 2016;71:1466–73. https://doi.org/10.1093/jac/dkw016.

    Article  CAS  PubMed  Google Scholar 

  27. Dejoies L, Sassi M, Schutz S, Moreaux J, Zouari A, Potrel S, et al. Genetic features of the poxtA linezolid resistance gene in human enterococci from France. J Antimicrob Chemother. 2021;76:1978–85. https://doi.org/10.1093/jac/dkab116.

    Article  CAS  PubMed  Google Scholar 

  28. Bender JK. Rapid emergence of highly variable and transferable oxazolidinone and phenicol resistance gene optrA in German Enterococcus spp. clinical isolates. Int J Antimicrob Agents. 2018; 9.

  29. Xuan H, Xia L, Schwarz S, Jia H, Yao X, Wang S, et al. Various mobile genetic elements carrying optrA in Enterococcus faecium and Enterococcus faecalis isolates from swine within the same farm. J Antimicrob Chemother. 2023;78:504–11. https://doi.org/10.1093/jac/dkac421.

    Article  CAS  PubMed  Google Scholar 

  30. Cai JC, Hu YY, Zhou HW, Chen G-X, Zhang R. Dissemination of the same cfr-carrying plasmid among methicillin-resistant Staphylococcus aureus and Coagulase-negative staphylococcal isolates in China. Antimicrob Agents Chemother. 2015;59:3669–71. https://doi.org/10.1128/AAC.04580-14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Deshpande LM, Castanheira M, Flamm RK, Mendes RE. Evolving oxazolidinone resistance mechanisms in a worldwide collection of enterococcal clinical isolates: results from the SENTRY Antimicrobial Surveillance Program. J Antimicrob Chemother. 2018;73:2314–22. https://doi.org/10.1093/jac/dky188.

    Article  CAS  PubMed  Google Scholar 

  32. Cui L, Wang Y, Lv Y, Wang S, Song Y, Li Y, et al. Nationwide Surveillance of Novel Oxazolidinone Resistance Gene optrA in Enterococcus Isolates in China from 2004 to 2014. Antimicrob Agents Chemother. 2016;60:7490–3. https://doi.org/10.1128/AAC.01256-16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Wardenburg KE, Potter RF, D’Souza AW, Hussain T, Wallace MA, Andleeb S, et al. Phenotypic and genotypic characterization of linezolid-resistant Enterococcus faecium from the USA and Pakistan. J Antimicrob Chemother. 2019;74:3445–52. https://doi.org/10.1093/jac/dkz367.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Moure Z, Lara N, Marín M, Sola-Campoy PJ, Bautista V, Gómez-Bertomeu F, et al. Interregional spread in spain of linezolid-resistant Enterococcus spp. isolates carrying the optrA and poxtA genes. Int J Antimicrob Agents. 2020;55:105977. https://doi.org/10.1016/j.ijantimicag.2020.105977.

    Article  CAS  PubMed  Google Scholar 

  35. Shen J, Wang Y, Schwarz S. Presence and dissemination of the multiresistance gene cfr in Gram-positive and Gram-negative bacteria. J Antimicrob Chemother. 2013;68:1697–706. https://doi.org/10.1093/jac/dkt092.

    Article  CAS  PubMed  Google Scholar 

  36. Cavaco LM, Bernal JF, Zankari E, Léon M, Hendriksen RS, Perez-Gutierrez E, et al. Detection of linezolid resistance due to the optrA gene in Enterococcus faecalis from poultry meat from the American continent (Colombia). J Antimicrob Chemother. 2016; dkw490. doi:https://doi.org/10.1093/jac/dkw490.

  37. Tyson GH, Sabo JL, Hoffmann M, Hsu C-H, Mukherjee S, Hernandez J, et al. Novel linezolid resistance plasmids in Enterococcus from food animals in the USA. J Antimicrob Chemother. 2018 [cited 27 Apr 2022]. doi:https://doi.org/10.1093/jac/dky369.

  38. Elghaieb H, Freitas AR, Abbassi MS, Novais C, Zouari M, Hassen A, et al. Dispersal of linezolid-resistant enterococci carrying poxtA or optrA in retail meat and food-producing animals from Tunisia. J Antimicrob Chemother. 2019;74:2865–9. https://doi.org/10.1093/jac/dkz263.

    Article  CAS  PubMed  Google Scholar 

  39. Lei C-W, Kang Z-Z, Wu S-K, Chen Y-P, Kong L-H, Wang H-N. Detection of the phenicol–oxazolidinone–tetracycline resistance gene poxtA in Enterococcus faecium and Enterococcus faecalis of food-producing animal origin in China. J Antimicrob Chemother. 2019;74:2459–61. https://doi.org/10.1093/jac/dkz198.

    Article  CAS  PubMed  Google Scholar 

  40. Almeida LM, Lebreton F, Gaca A, Bispo PM, Saavedra JT, Calumby RN, et al. Transferable resistance gene optrA in Enterococcus faecalis from swine in Brazil. Antimicrob Agents Chemother. 2020;64:11.

    Article  Google Scholar 

  41. Timmermans M, Bogaerts B, Vanneste K, De Keersmaecker SCJ, Roosens NHC, Kowalewicz C, et al. Large diversity of linezolid-resistant isolates discovered in food-producing animals through linezolid selective monitoring in Belgium in 2019. J Antimicrob Chemother. 2021;77:49–57. https://doi.org/10.1093/jac/dkab376.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Huang Z, Bai Y, Wang Q, Yang X, Zhang T, Chen X, et al. Persistence of transferable oxazolidinone resistance genes in enterococcal isolates from a swine farm in China. Front Microbiol. 2022;13:1010513. https://doi.org/10.3389/fmicb.2022.1010513.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Abdullahi IN, Lozano C, Juárez-Fernández G, Höfle U, Simón C, Rueda S, et al. Nasotracheal enterococcal carriage and resistomes: detection of optrA-, poxtA- and cfrD-carrying strains in migratory birds, livestock, pets, and in-contact humans in Spain. Eur J Clin Microbiol Infect Dis. 2023;42:569–81. https://doi.org/10.1007/s10096-023-04579-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Nüesch-Inderbinen M, Heyvaert L, Treier A, Zurfluh K, Cernela N, Biggel M, et al. High occurrence of Enterococcus faecalis, Enterococcus faecium, and Vagococcus lutrae harbouring oxazolidinone resistance genes in raw meat-based diets for companion animals – a public health issue, Switzerland, September 2018 to May 2020. Eurosurveillance. 2023;28https://doi.org/10.2807/1560-7917.ES.2023.28.6.2200496.

  45. EFFORT Group, Munk P, Knudsen BE, Lukjancenko O, Duarte ASR, Van Gompel L, et al. Abundance and diversity of the faecal resistome in slaughter pigs and broilers in nine European countries. Nat Microbiol. 2018;3:898–908. https://doi.org/10.1038/s41564-018-0192-9.

    Article  CAS  Google Scholar 

  46. Zhao Q, Wang Y, Wang S, Wang Z, Du X, Jiang H, et al. Prevalence and abundance of florfenicol and linezolid resistance genes in soils adjacent to swine feedlots. Sci Rep. 2016;6:32192. https://doi.org/10.1038/srep32192.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  47. Freitas AR, Elghaieb H, León-Sampedro R, Abbassi MS, Novais C, Coque TM, et al. Detection of optrA in the African continent (Tunisia) within a mosaic Enterococcus faecalis plasmid from urban wastewaters. J Antimicrob Chemother. 2017;72:3245–51. https://doi.org/10.1093/jac/dkx321.

    Article  CAS  PubMed  Google Scholar 

  48. Gu J, Xie X-J, Liu J-X, Shui J-R, Zhang H-Y, Feng G-Y, et al. Prevalence and transmission of antimicrobial-resistant Staphylococci and Enterococci from shared bicycles in Chengdu. China Sci Total Environ. 2020;738:139735. https://doi.org/10.1016/j.scitotenv.2020.139735.

    Article  ADS  CAS  PubMed  Google Scholar 

  49. Biggel M, Nüesch-Inderbinen M, Jans C, Stevens MJA, Stephan R. Genetic context of optrA and poxtA in florfenicol-resistant enterococci isolated from flowing surface water in Switzerland. Antimicrob Agents Chemother. 2021;65:e01083-e1121. https://doi.org/10.1128/AAC.01083-21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Fioriti S, Coccitto SN, Cedraro N, Simoni S, Morroni G, Brenciani A, et al. Linezolid resistance genes in enterococci isolated from sediment and zooplankton in two Italian coastal areas. Björkroth J, editor. Appl Environ Microbiol. 2021;87:e02958-20. https://doi.org/10.1128/AEM.02958-20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Yu Y. Lilium spp., as unnoticed environmental vector, spreading OptrA-carrying Enterococcus spp. Sci Total Environ. 2022;816:151540. https://doi.org/10.1016/j.scitotenv.2021.151540.

    Article  ADS  CAS  PubMed  Google Scholar 

  52. Shen W, Huang Y, Cai J. An optimized screening approach for the oxazolidinone resistance gene optrA yielded a higher fecal carriage rate among healthy individuals in Hangzhou, China. Tuon FF, editor. Microbiol Spectr. 2022;10:e02974-22. https://doi.org/10.1128/spectrum.02974-22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Clinical and Laboratory Standards Institute. Methods for dilution antimicrobial susceptibility tests for bacteria that grow aerobically, 11th edition CLSI standard M07. 2018; Clinical and Laboratory Standards Institute, Wayne, PA.

  54. Clinical and Laboratory Standards Institute. Performance standards for antimicrobial susceptibility testing, 32nd edition. CLSI supplement M100. 2022: Clinical and Laboratory Standards Institute, Wayne, PA.

  55. Clinical and Laboratory Standards Institute, 2015. Methods for antimicrobial dilution and disk susceptibility testing of infrequently isolated or fastidious bacteria, approved standard, third edition M45. Wayne, PA.

  56. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19:455–77. https://doi.org/10.1089/cmb.2012.0021.

    Article  MathSciNet  CAS  PubMed  PubMed Central  Google Scholar 

  57. Bortolaia V, Kaas RS, Ruppe E, Roberts MC, Schwarz S, Cattoir V, et al. ResFinder 4.0 for predictions of phenotypes from genotypes. J Antimicrob Chemother. 2020;75:3491–500. https://doi.org/10.1093/jac/dkaa345.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Larsen MV, Cosentino S, Rasmussen S, Friis C, Hasman H, Marvig RL, et al. Multilocus sequence typing of total-genome-sequenced bacteria. J Clin Microbiol. 2012;50:1355–61. https://doi.org/10.1128/JCM.06094-11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Katoh K, Standley DM. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in performance and usability. Mol Biol Evol. 2013;30:772–80. https://doi.org/10.1093/molbev/mst010.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Chernomor O, Von Haeseler A, Minh BQ. Terrace aware data structure for phylogenomic inference from supermatrices. Syst Biol. 2016;65:997–1008. https://doi.org/10.1093/sysbio/syw037.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Treangen TJ, Ondov BD, Koren S, Phillippy AM. The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes. 2014; 15.

  62. Letunic I, Bork P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016;44:W242–5. https://doi.org/10.1093/nar/gkw290.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Sullivan MJ, Petty NK, Beatson SA. Easyfig: a genome comparison visualizer. Bioinformatics. 2011;27:1009–10. https://doi.org/10.1093/bioinformatics/btr039.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Wang Y. Association of florfenicol residues with the abundance of oxazolidinone resistance genes in livestock manures. J Hazard Mater. 2020;399:123059. https://doi.org/10.1016/j.jhazmat.2020.123059.

    Article  CAS  PubMed  Google Scholar 

  65. Torres C, Alonso CA, Ruiz-Ripa L, León-Sampedro R, Del Campo R, Coque TM. Antimicrobial Resistance in Enterococcus spp. of animal origin. Aarestrup FM, Schwarz S, Shen J, Cavaco L, editors. Microbiol Spectr. 2018;6: 6.4.24. https://doi.org/10.1128/microbiolspec.ARBA-0032-2018.

  66. Sørensen TL, Blom M, Monnet DL, Frimodt-Møller N, Poulsen RL, Espersen F. Transient Intestinal Carriage after Ingestion of Antibiotic-Resistant Enterococcus faecium from Chicken and Pork. N Engl J Med. 2001;345:1161–6. https://doi.org/10.1056/NEJMoa010692.

    Article  PubMed  Google Scholar 

  67. Klare I, Konstabel C, Werner G, Huys G, Vankerckhoven V, Kahlmeter G, et al. Antimicrobial susceptibilities of Lactobacillus, Pediococcus and Lactococcus human isolates and cultures intended for probiotic or nutritional use. J Antimicrob Chemother. 2007;59:900–12. https://doi.org/10.1093/jac/dkm035.

    Article  CAS  PubMed  Google Scholar 

  68. Chen H, Wang X, Yin Y, Li S, Zhang Y, Wang Q, et al. Molecular characteristics of oxazolidinone resistance in enterococci from a multicenter study in China. BMC Microbiol. 2019;19:162. https://doi.org/10.1186/s12866-019-1537-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Funding

This work was supported by the National Key Research and Development Program of China (grant no. 2022YFD1800403) and the National Natural Science Foundation of China (grant no. 81971988).

Author information

Authors and Affiliations

Authors

Contributions

JC (Jiachang Cai) and RZ conceived and designed the project. JC (Jiachang Cai) and JC (Jiawei Chen) collected bacterial strains. JC (Jiachang Cai) and WS conducted the experiments. WS and ND conducted the bioinformatic computation. WS drafted the manuscript. CC, ND, RZ, and JC (Jiachang Cai) supervised the project, edited the manuscript, and verified the data. All authors commented on a draft of the manuscript and approved the final version.

Corresponding authors

Correspondence to Rong Zhang or Jiachang Cai.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Supplementary Table S1. Characteristics of strains carrying linezolid resistance genes from multiple niches. Detailed information of linezolid resistance reservoirs isolated in this study regarding their source, species, antimicrobial susceptibility profiles, carriage of antimicrobial resistance genes, and genetic context of linezolid resistance genes were illustrated.

Additional file 2:

Supplementary Figure S1. Representation of alignment of the amino acid sequence and relevant linezolid MIC values for OptrA variants. Substitutions (green square) in amino sequence, when compared to the wild-type OptrA sequence (yellow square), were illustrated as indicated. The numbers in the blue square represented the numbers for strains that displayed the corresponding linezolid MICs of certain species, while the grey ones represented no corresponding strains in our study.

Additional file 3:

Supplementary Figure S2. Phylogenetic maximum-likelihood tree of linezolid resistance reservoirs. The SNP-based phylogenetic tree for E. faecalis (a), E. faecium (b), L. lactis (c), E. avium (d), E. dongliensis (e), L. petauri (f), E. gallinarum (g), V. lutrae (h), E. casseliflavus (i), E. raffinosus (j), E. hirae (k), V. carniphilus (l), and E. hulanensis (m). The sample type, source, and OptrA variants were illustrated as the legend indicated. The * represented the two types of OptrA variants in the strain as indicated. The tiny squares with linear arrangements from the inside to the outside indicated the carriage of antimicrobial resistance genes as the label indicated. The filled squares represented the presence of this gene in the genome of this isolate and vice versa. The strains closely related were labeled in the yellow or green shade.

Additional file 4:

Supplementary Figure S3. The genetic environment of the linezolid resistance genes. The arrows indicate the positions and directions of transcription of the different genes. All insertion sequences were labeled as orange, antimicrobial resistance genes were labeled as red, gene encoding proteins with known functions were labeled as blue, the rep genes were labeled as green, and the hypothesis proteins were labeled as grey. The letters at the lower right corner of optrA indicated the OptrA variants. N represented the number of strains that harbored the corresponding resistance locus encoding the gene or variants. Δ symbol indicates a truncated gene.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Shen, W., Cai, C., Dong, N. et al. Mapping the widespread distribution and transmission dynamics of linezolid resistance in humans, animals, and the environment. Microbiome 12, 52 (2024). https://doi.org/10.1186/s40168-023-01744-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40168-023-01744-2

Keywords