ARGNet: using deep neural networks for robust identification and classification of antibiotic resistance genes from sequences.

BACKGROUND: Emergence of antibiotic resistance in bacteria is an important threat to global health. Antibiotic resistance genes (ARGs) are some of the key components to define bacterial resistance and their spread in different environments. Identification of ARGs, particularly from high-throughput sequencing data of the specimens, is the state-of-the-art method for comprehensively monitoring their spread and evolution. Current computational methods to identify ARGs mainly rely on alignment-based sequence similarities with known ARGs. Such approaches are limited by choice of reference databases and may potentially miss novel ARGs. The similarity thresholds are usually simple and could not accommodate variations across different gene families and regions. It is also difficult to scale up when sequence data are increasing. RESULTS: In this study, we developed ARGNet, a deep neural network that incorporates an unsupervised learning autoencoder model to identify ARGs and a multiclass classification convolutional neural network to classify ARGs that do not depend on sequence alignment. This approach enables a more efficient discovery of both known and novel ARGs. ARGNet accepts both amino acid and nucleotide sequences of variable lengths, from partial (30-50 aa; 100-150 nt) sequences to full-length protein or genes, allowing its application in both target sequencing and metagenomic sequencing. Our performance evaluation showed that ARGNet outperformed other deep learning models including DeepARG and HMD-ARG in most of the application scenarios especially quasi-negative test and the analysis of prediction consistency with phylogenetic tree. ARGNet has a reduced inference runtime by up to 57% relative to DeepARG. CONCLUSIONS: ARGNet is flexible, efficient, and accurate at predicting a broad range of ARGs from the sequencing data. ARGNet is freely available at https://github.com/id-bioinfo/ARGNet , with an online service provided at https://ARGNet.hku.hk . Video Abstract.

Characterization of Neowestiellopsis persica A1387 (Hapalosiphonaceae) based on the cpcA, psbA, rpoC1, nifH and nifD gene sequences.

BACKGROUND: Complex descriptions of new strains of cyanobacteria appear very frequently. The main importance of these descriptions concerns potential new substances that they could synthesise, as well as their different properties as a result of their different ecological niches. The main gene used for these descriptions is 16 S with ITS or whole genome sequencing. Neowestiellopsis persica represents a unique example of the influence of ecology on morphological changes, with almost identical 16 S identity. Although our previously described Neowestiellopsis persica strain A1387 was characterized by 16 S analysis, we used different molecular markers to provide a way to separate strains of this genus that are closely related at the genetic level. MATERIALS AND METHODS: In order to conduct an in-depth study, several molecular markers, namely psbA, rpoC1, nifD, nifH and cpcA were sequenced and studied in Neowestiellopsis persica strain A1387. RESULTS: The results of the phylogenetic analysis, based on cpcA, showed that the studied strain A 1387 falls into a separate clade than N. persica, indicating that this signature sequence could be a useful molecular marker for phylogenetic separation of similar strains isolated in the future. CONCLUSIONS: Analysis of strain A1387 based on gene differences confirmed that it is a Neowestiellopsis strain. The morphological changes observed in the previous study could be due to different ecological and cultivation conditions compared to the type species. At the same time, the sequences obtained have increased our understanding of this species and will help in the future to better identify strains belonging to the genus Neowestiellopsis.

Identification of an Enterococcus faecium strain isolated from raw bovine milk co-harbouring the oxazolidinone resistance genes optrA and poxtA in China.

Oxazolidinones are potent antimicrobial agents used to treat human infections caused by multidrug-resistant Gram-positive bacteria. The growing resistance to oxazolidinones poses a significant threat to public health. In August 2021, a linezolid-resistant Enterococcus faecium BN83 was isolated from a raw milk sample of cow in Inner Mongolia, China. This isolate exhibited a multidrug resistance phenotype and was resistant to most of drugs tested including linezolid and tedizolid. PCR detection showed that two mobile oxazolidinones resistance genes, optrA and poxtA, were present in this isolate. Whole genome sequencing analysis revealed that the genes optrA and poxtA were located on two different plasmids, designated as pBN83-1 and pBN83-2, belonging to RepA_N and Inc18 families respectively. Genetic context analysis suggested that optrA gene on plasmid pBN83-1 was located in transposon Tn6261 initially found in E. faecalis. Comprehensive analysis revealed that Tn6261 act as an important horizontal transmission vector for the spread of optrA in E. faecium. Additionally, poxtA-bearing pBN83-2 displayed high similarity to numerous plasmids from Enterococcus of different origin and pBN83-2-like plasmid represented a key mobile genetic element involved in movement of poxtA in enterococcal species. The presence of optrA- and poxtA-carrying E. faecium in raw bovine milk represents a public health concern and active surveillance is urgently warranted to investigate the prevalence of oxazolidinone resistance genes in animal-derived food products.

Co-localization of antibiotic resistance genes is widespread in the infant gut microbiome and associates with an immature gut microbial composition.

BACKGROUND: In environmental bacteria, the selective advantage of antibiotic resistance genes (ARGs) can be increased through co-localization with genes such as other ARGs, biocide resistance genes, metal resistance genes, and virulence genes (VGs). The gut microbiome of infants has been shown to contain numerous ARGs, however, co-localization related to ARGs is unknown during early life despite frequent exposures to biocides and metals from an early age. RESULTS: We conducted a comprehensive analysis of genetic co-localization of resistance genes in a cohort of 662 Danish children and examined the association between such co-localization and environmental factors as well as gut microbial maturation. Our study showed that co-localization of ARGs with other resistance and virulence genes is common in the early gut microbiome and is associated with gut bacteria that are indicative of low maturity. Statistical models showed that co-localization occurred mainly in the phylum Proteobacteria independent of high ARG content and contig length. We evaluated the stochasticity of co-localization occurrence using enrichment scores. The most common forms of co-localization involved tetracycline and fluoroquinolone resistance genes, and, on plasmids, co-localization predominantly occurred in the form of class 1 integrons. Antibiotic use caused a short-term increase in mobile ARGs, while non-mobile ARGs showed no significant change. Finally, we found that a high abundance of VGs was associated with low gut microbial maturity and that VGs showed even higher potential for mobility than ARGs. CONCLUSIONS: We found that the phenomenon of co-localization between ARGs and other resistance and VGs was prevalent in the gut at the beginning of life. It reveals the diversity that sustains antibiotic resistance and therefore indirectly emphasizes the need to apply caution in the use of antimicrobial agents in clinical practice, animal husbandry, and daily life to mitigate the escalation of resistance. Video Abstract.

Investigating antimicrobial resistance genes in Kenya, Uganda and Tanzania cattle using metagenomics.

Antimicrobial resistance (AMR) is a growing problem in African cattle production systems, posing a threat to human and animal health and the associated economic value chain. However, there is a poor understanding of the resistomes in small-holder cattle breeds in East African countries. This study aims to examine the distribution of antimicrobial resistance genes (ARGs) in Kenya, Tanzania, and Uganda cattle using a metagenomics approach. We used the SqueezeMeta-Abricate (assembly-based) pipeline to detect ARGs and benchmarked this approach using the Centifuge-AMRplusplus (read-based) pipeline to evaluate its efficiency. Our findings reveal a significant number of ARGs of critical medical and economic importance in all three countries, including resistance to drugs of last resort such as carbapenems, suggesting the presence of highly virulent and antibiotic-resistant bacterial pathogens (ESKAPE) circulating in East Africa. Shared ARGs such as aph(6)-id (aminoglycoside phosphotransferase), tet (tetracycline resistance gene), sul2 (sulfonamide resistance gene) and cfxA_gen (betalactamase gene) were detected. Assembly-based methods revealed fewer ARGs compared to read-based methods, indicating the sensitivity and specificity of read-based methods in resistome characterization. Our findings call for further surveillance to estimate the intensity of the antibiotic resistance problem and wider resistome classification. Effective management of livestock and antibiotic consumption is crucial in minimizing antimicrobial resistance and maximizing productivity, making these findings relevant to stakeholders, agriculturists, and veterinarians in East Africa and Africa at large.

Identification of the novel fosfomycin resistance gene fosSC in Staphylococcus capitis.

OBJECTIVES: Fosfomycin has regained attention for treating infections caused by methicillin-resistant Staphylococcus aureus and multidrug-resistant coagulase-negative staphylococci. In this research, our objective was to investigate the mechanisms underlying fosfomycin resistance in Staphylococcus capitis. METHODS: The minimum inhibitory concentrations (MICs) of fosfomycin were assessed in 109 clinical S. capitis isolates by the agar dilution method. By cloning the fos-like genes into the shuttle vector, pTSSCm-Pcap, and observing the change in fosfomycin MICs, the gene function was verified. Core genome multilocus sequence typing and comparative genomics analysis were conducted to determine the population characteristics of S. capitis isolates and analyse the genetic environment of the fos-like genes. RESULTS: We identified a novel fosfomycin resistance gene, fosSC, on the chromosome in 58 out of 109 (53.2%) S. capitis isolates. The deduced products of the fosSC genes shared 67.15-67.88% amino acid sequence identity with FosB. The RN-pT-fosSC transformants carrying fosSC showed a 512-fold increase in the fosfomycin MICs. The fosSC gene was embedded in a conserved genetic context, but IS431mec was located to the left of the fosSC gene in cluster L due to the insertion of staphylococcal cassette chromosome mec. CONCLUSIONS: The chromosomal fosSC genes in some lineages of S. capitis explained their high-level fosfomycin resistance. Ongoing surveillance is crucial for monitoring the potential threat of horizontal transfer, which could be facilitated by the presence of mobile genetic elements surrounding the fosSC gene.

Community coalescence and plant host filtering determine the spread of tetracycline resistance genes from pig manure into the microbiome continuum of the soil-plant system.

The spread of livestock manure-borne antibiotic resistance genes (ARGs) into agroecosystems through manure application poses a potential threat to human health. However, there is still a knowledge gap concerning ARG dissemination in coalescing manure, soil and plant microbiomes. Here, we examined the fate of tetracycline resistance genes (TRGs) originating from pig manure microbiomes and spread in the soil-A thaliana system and explored the effects of microbial functions on TRGs spread at different interfaces. Our results indicate that the TRGs abundances in all microbiome continuum of the soil-A. thaliana system were significantly increased with the application of a living manure microbiome, although the addition of manure with both an active and inactive microbiome caused a shift in the microbial community composition. This was attributed to the increasing relative abundances of tetA, tetL, tetM, tetO, tetW and tolC in the system. The application of living manure with DOX residues resulted in the highest relative abundance of total TRGs (3.30×10-3 copies/16S rRNA gene copies) in the rhizosphere soil samples. Community coalescence of the manure and soil microbiomes increased the abundance of Firmicutes in the soil and root microbiome, which directly explains the increase in TRG abundance observed in these interfaces. In contrast, the leaf microbiome differed markedly from that of the remaining samples, indicating strong plant host filtering effects on Firmicutes and TRGs from pig manure. The random forest machine learning model revealed microbial functions and their significant positive correlation with TRG abundance in the microbiome continuum of the system. Our findings revealed that community coalescence is the main driver of TRG spread from manure to the soil and root microbiomes. Plant host filtering effects play a crucial role in allowing certain microbial groups to occupy ecological niches in the leaves, thereby limiting the establishment of manure-borne TRGs in aboveground plant tissues.

What role do biocontrol agents with Mg2+ play in the fate of antibiotic resistome and pathogenic bacteria in the phyllosphere?

The contamination of the plant phyllosphere with antibiotics and antibiotic resistance genes (ARGs), caused by application of antibiotics, is a significant environmental issue in agricultural management. Alternatively, biocontrol agents are environmentally friendly and have attracted a lot of interest. However, the influence of biocontrol agents on the phyllosphere resistome remains unknown. In this study, we applied biocontrol agents to control the wildfire disease in the Solanaceae crops and investigated their effects on the resistome and the pathogen in the phyllosphere by using metagenomics. A total of 250 ARGs were detected from 15 samples, which showed a variation in distribution across treatments of biocontrol agents (BA), BA with Mg2+ (T1), BA with Mn2+ (T2), and kasugamycin (T3) and nontreated (CK). The results showed that the abundance of ARGs under the treatment of BA-Mg2+ was lower than that in the CK group. The abundance of cphA3 (carbapenem resistance), PME-1 (carbapenem resistance), tcr3 (tetracycline antibiotic resistance), and AAC (3)-VIIIa (aminoglycoside antibiotic resistance) in BA-Mg2+ was significantly higher than that in BA-Mn2+ (P < 0.05). The abundance of cphA3, PME_1, and tcr3 was significantly negatively related to the abundance of the phyllosphere pathogen Pseudomonas syringae (P < 0.05). We also found that the upstream and downstream regions of cphA3 were relatively conserved, in which rpl, rpm, and rps gene families were identified in most sequences (92%). The Ka/Ks of cphA3 was 0 in all observed sequences, indicating that under the action of purifying selection, nonsynonymous substitutions are often gradually eliminated in the population. Overall, this study clarifies the effect of biocontrol agents with Mg2+ on the distribution of the phyllosphere resistome and provides evolutionary insights into the biocontrol process. IMPORTANCE: Our study applied metagenomics analysis to examine the impact of biocontrol agents (BAs) on the phyllosphere resistome and the pathogen. Irregular use of antibiotics has led to the escalating dissemination of antibiotic resistance genes (ARGs) in the environment. The majority of BA research has focused on the effect of monospecies on the plant disease control process, the role of the compound BA with nutrition elements in the phyllosphere disease, and the resistome is still unknown. We believe BAs are eco-friendly alternatives for antibiotics to combat the transfer of ARGs. Our results revealed that BA-Mg2+ had a lower relative abundance of ARGs compared to the CK group, and the phyllosphere pathogen Pseudomonas syringae was negatively related to three specific ARGs, cphA3, PME-1, and tcr3. These three genes also present different Ka/Ks. We believe that the identification of the distribution and evolution modes of ARGs further elucidates the ecological role and facilitates the development of BAs, which will attract general interest in this field.

Microbial biofilms on macroalgae harbour diverse integron gene cassettes.

Integrons are genetic platforms that capture, rearrange and express mobile modules called gene cassettes. The best characterized gene cassettes encode antibiotic resistance, but the function of most integron gene cassettes remains unknown. Functional predictions suggest that many gene cassettes could encode proteins that facilitate interactions with other cells and with the extracellular environment. Because cell interactions are essential for biofilm stability, we sequenced gene cassettes from biofilms growing on the surface of the marine macroalgae Ulva australis and Sargassum linearifolium. Algal samples were obtained from coastal rock platforms around Sydney, Australia, using seawater as a control. We demonstrated that integrons in microbial biofilms did not sample genes randomly from the surrounding seawater, but harboured specific functions that potentially provided an adaptive advantage to both the bacterial cells in biofilm communities and their macroalgal host. Further, integron gene cassettes had a well-defined spatial distribution, suggesting that each bacterial biofilm acquired these genetic elements via sampling from a large but localized pool of gene cassettes. These findings suggest two forms of filtering: a selective acquisition of different integron-containing bacterial species into the distinct biofilms on Ulva and Sargassum surfaces, and a selective retention of unique populations of gene cassettes at each sampling location.

Aminoglycoside resistance genes in early members of the Acinetobacter baumannii ST78A (SMAL, Italian clone) reside in an IS26-bounded island in the chromosome.

BACKGROUND: The Acinetobacter baumannii isolate called SMAL, previously used to determine the structures of capsular polysaccharide and lipooligosaccharide, was recovered in Pavia, Italy in 2002 among the collection of aminoglycoside-resistant isolates designated as SMAL type. This type was later called the Italian clone, then ST78. ST78 isolates are now widely distributed. OBJECTIVES: To establish the resistance gene complement and the location and structure of acquired resistance regions in early members of the Italian/ST78 clone. METHODS: The draft genome of SMAL2002 was assembled from Illumina MiSeq reads. Contigs containing resistance genes were joined and located in the chromosome using PCR with custom primers. The resistance profile was determined using disc diffusion. RESULTS: SMAL2002 is an ST78A isolate and includes three aminoglycoside resistance genes, aadB (gentamicin, kanamycin, tobramycin) aphA1 (kanamycin, neomycin) and aac(6')-Ian (amikacin, kanamycin, tobramycin). The aadB gene cassette is incorporated at a secondary site in a relative of the aphA1-containing, IS26-bounded pseudo-compound transposon, PTn6020. The aac(6')-Ian gene is in an adjacent IS26-bounded structure that includes sul2 (sulphonamide) and floR (florfenicol) resistance genes. The two pseudo-compound transposons overlap and are in the chromosomal hutU gene flanked by an 8 bp target site duplication. Although aac(6')-Ian was not noticed previously, the same genes and structures were found in several available draft genomes of early ST78A isolates. CONCLUSIONS: This study highlights the importance of correlating resistance profiles with resistance gene content. The location of acquired resistance genes in the SMAL2002 chromosome represents the original location in the ST78A lineage of ST78.

Development of an ultrafast PCR to detect clinically relevant acquired vancomycin-resistance genes from cultured enterococci.

BACKGROUND: VRE are increasingly described worldwide. Screening of hospitalized patients at risk for VRE carriage is mandatory to control their dissemination. Here, we have developed the Bfast [VRE Panel] PCR kit, a rapid and reliable quantitative PCR assay for detection of vanA, vanB, vanD and vanM genes, from solid and liquid cultures adaptable to classical and ultrafast real-time PCR platforms. METHODS: Validation was carried out on 133 well characterized bacterial strains, including 108 enterococci of which 64 were VRE. Analytical performances were determined on the CFX96 Touch (Bio-Rad) and Chronos Dx (BforCure), an ultrafast qPCR machine. Widely used culture plates and broths for enterococci selection/growth were tested. RESULTS: All targeted van alleles (A, B, D and M) were correctly detected without cross-reactivity with other van genes (C, E, G, L and N) and no interference with the different routinely used culture media. A specificity and sensitivity of 100% and 99.7%, respectively, were determined, with limits of detection ranging from 21 to 238 cfu/reaction depending on the targets. The Bfast [VRE Panel] PCR kit worked equally well on the CFX and Chronos Dx platforms, with differences in multiplexing capacities (five and four optical channels, respectively) and in turnaround time (45 and 16 minutes, respectively). CONCLUSIONS: The Bfast [VRE Panel] PCR kit is robust, easy to use, rapid and easily implementable in clinical microbiology laboratories for ultra-rapid confirmation of the four main acquired van genes. Its features, especially on Chronos Dx, seem to be unmatched compared to other tools for screening of VRE.

Identification of ISVlu1-derived translocatable units containing optrA and/or fexA genes generated by homologous or illegitimate recombination in Lactococcus garvieae of porcine origin.

The optrA gene encodes an ABC-F protein which confers cross-resistance to oxazolidinones and phenicols. Insertion sequence ISVlu1, a novel ISL3-family member, was recently reported to be involved in the transmission of optrA in Vagococcus lutrae. However, the role of ISVlu1 in mobilizing resistance genes has not yet fully explored. In this study, two complete and three truncated copies of ISVlu1 were found on plasmid pBN62-optrA from Lactococcus garvieae. Analysis of the genetic context showed that both optrA and the phenicols resistance gene fexA were flanked by the complete or truncated ISVlu1 copies. Moreover, three different-sized ISVlu1-based translocatable units (TUs) carrying optrA and/or fexA, were detected from pBN62-optrA. Sequence analysis revealed that the TU-optrA was generated by homologous recombination while TU-fexA and TU-optrA+fexA were the products of illegitimate recombinations. Importantly, conjugation assays confirmed that pBN62-optrA was able to successfully transfer into the recipient Enterococcus faecalis JH2-2. To our knowledge, this is the first report about an optrA-carrying plasmid in L. garvieae which could horizontally transfer into other species. More importantly, the ISVlu1-flanked genetic structures containing optrA and/or fexA were also observed in bacteria of different species, which underlines that ISVlu1 is highly active and plays a vital role in the transfer of some important resistance genes, such as optrA and fexA.

Antibiotics and antibiotic resistance genes removal in biological aerated filter.

Two laboratory-level biological aerated filters (BAF) were constructed to explore their treatment capacity for simulated antibiotic wastewater at high (1 - 16 mg/L) and low (0 - 0.5 mg/L) concentrations. Results showed that BAF was capable of removing both sulfonamides and tetracyclines with an efficiency of over 90 % at 16 mg/L. The main mechanism for removing antibiotics was found to be biodegradation followed by adsorption. Paenarthrobacter was identified as the key genus in sulfonamides degradation, while Hydrogenophaga played a crucial role in tetracyclines degradation. Antibiotics resistant genes such as intI1, sul1, sul2, tetA, tetW and tetX were frequently detected in the effluent, with interception rates ranging from 105 - 106 copies/mL. The dominated microorganisms obtained in the study could potentially be utilized to enhance the capacity of biological processes for treating antibiotics contaminated wastewater. These findings contribute to a better understanding of BAF treating wastewater containing antibiotics and resistant genes.

Duplicated antibiotic resistance genes reveal ongoing selection and horizontal gene transfer in bacteria.

Horizontal gene transfer (HGT) and gene duplication are often considered as separate mechanisms driving the evolution of new functions. However, the mobile genetic elements (MGEs) implicated in HGT can copy themselves, so positive selection on MGEs could drive gene duplications. Here, we use a combination of modeling and experimental evolution to examine this hypothesis and use long-read genome sequences of tens of thousands of bacterial isolates to examine its generality in nature. Modeling and experiments show that antibiotic selection can drive the evolution of duplicated antibiotic resistance genes (ARGs) through MGE transposition. A key implication is that duplicated ARGs should be enriched in environments associated with antibiotic use. To test this, we examined the distribution of duplicated ARGs in 18,938 complete bacterial genomes with ecological metadata. Duplicated ARGs are highly enriched in bacteria isolated from humans and livestock. Duplicated ARGs are further enriched in an independent set of 321 antibiotic-resistant clinical isolates. Our findings indicate that duplicated genes often encode functions undergoing positive selection and horizontal gene transfer in microbial communities.

Modeling the limits of detection for antimicrobial resistance genes in agri-food samples: a comparative analysis of bioinformatics tools.

BACKGROUND: Although the spread of antimicrobial resistance (AMR) through food and its production poses a significant concern, there is limited research on the prevalence of AMR bacteria in various agri-food products. Sequencing technologies are increasingly being used to track the spread of AMR genes (ARGs) in bacteria, and metagenomics has the potential to bypass some of the limitations of single isolate characterization by allowing simultaneous analysis of the agri-food product microbiome and associated resistome. However, metagenomics may still be hindered by methodological biases, presence of eukaryotic DNA, and difficulties in detecting low abundance targets within an attainable sequence coverage. The goal of this study was to assess whether limits of detection of ARGs in agri-food metagenomes were influenced by sample type and bioinformatic approaches. RESULTS: We simulated metagenomes containing different proportions of AMR pathogens and analysed them for taxonomic composition and ARGs using several common bioinformatic tools. Kraken2/Bracken estimates of species abundance were closest to expected values. However, analysis by both Kraken2/Bracken indicated presence of organisms not included in the synthetic metagenomes. Metaphlan3/Metaphlan4 analysis of community composition was more specific but with lower sensitivity than the Kraken2/Bracken analysis. Accurate detection of ARGs dropped drastically below 5X isolate genome coverage. However, it was sometimes possible to detect ARGs and closely related alleles at lower coverage levels if using a lower ARG-target coverage cutoff (< 80%). While KMA and CARD-RGI only predicted presence of expected ARG-targets or closely related gene-alleles, SRST2 (which allows read to map to multiple targets) falsely reported presence of distantly related ARGs at all isolate genome coverage levels. The presence of background microbiota in metagenomes influenced the accuracy of ARG detection by KMA, resulting in mcr-1 detection at 0.1X isolate coverage in the lettuce but not in the beef metagenome. CONCLUSIONS: This study demonstrates accurate detection of ARGs in synthetic metagenomes using various bioinformatic methods, provided that reads from the ARG-encoding organism exceed approximately 5X isolate coverage (i.e. 0.4% of a 40 million read metagenome). While lowering thresholds for target gene detection improved sensitivity, this led to the identification of alternative ARG-alleles, potentially confounding the identification of critical ARGs in the resistome. Further advancements in sequencing technologies providing increased coverage depth or extended read lengths may improve ARG detection in agri-food metagenomic samples, enabling use of this approach for tracking clinically important ARGs in agri-food samples.

Surviving the Potato Stems: Differences in Genes Required for Fitness by Dickeya dadantii and Dickeya dianthicola.

Bacteria belonging to the genus Dickeya cause blackleg and soft rot symptoms on many plant hosts, including potato. Although there is considerable knowledge about the genetic determinants that allow Dickeya to colonize host plants, as well as the genes that contribute to virulence, much is still unknown. To identify the genes important for fitness in potato stems, we constructed and evaluated randomly barcoded transposon mutant (RB-TnSeq) libraries of Dickeya dadantii and Dickeya dianthicola. We identified 169 and 157 genes important for growth in D. dadantii and D. dianthicola in stems, respectively. This included genes related to metabolic pathways, chemotaxis and motility, transcriptional regulation, transport across membranes, membrane biogenesis, detoxification mechanisms, and virulence-related genes, including a potential virulence cluster srfABC, c-di-GMP modulating genes, and pectin degradation genes. When we compared the results of the stem assay with other datasets, we identified genes important for growth in stems versus tubers and in vitro conditions. Additionally, our data showed differences in fitness determinants for D. dadantii and D. dianthicola. These data provide important insights into the mechanisms used by Dickeya when interacting with and colonizing plants and thus might provide targets for management.

[High-throughput qPCR and Amplicon Sequencing as Complementary Methods for Profiling Antibiotic Resistance Genes in Urban Wetland Parks].

Urban wetland parks are an important practice for urban wetland protection and utilization due to the vast ecosystem service value. As emerging contaminants, antibiotic resistance genes (ARGs) are great attractions for environmental research and public concerns. Based on high-throughput qPCR and high-throughput amplicon sequencing techniques, we investigated the occurrence, abundance, and distribution profiles of antibiotic resistance genes in the aquatic environment of Xiamen urban wetland parks (five sites). The influencing factors and driving mechanisms of antibiotic resistance genes were deciphered on the basis of microbial community structure and water quality. Diverse and abundant ARGs were observed and coexisted in urban wet parks. A total of 217 ARGs were detected in the water body of urban wetland parks, with an abundance up to 6.48×109 copies·L-1. Urban wetland parks were important hotspots and repositories of the antibiotic resistome. A total of nine bacterial genera, including Marivivens, NS5_marine_group, and Planktomarina, were identified as the potential carriers of diverse resistance genes (41 ARGs). The microbial communities could alone explain 51% of alterations in the antibiotic resistome in the aquatic environment of the urban wetland parks. Therefore, the microbial community was the key driving force for the occurrence and evolution of ARGs in urban wetland parks. Based on the results, with the presence of ARGs and antibiotic resistance bacteria, it is suggested that the water environments of urban wetland parks have potential risks of water ecological security and human health, and it is necessary to further enhance the research and control of microbial contaminants in the aquatic environment of urban wetland parks.

Genome-centric analyses of 165 metagenomes show that mobile genetic elements are crucial for the transmission of antimicrobial resistance genes to pathogens in activated sludge and wastewater.

Wastewater is considered a reservoir of antimicrobial resistance genes (ARGs), where the abundant antimicrobial-resistant bacteria and mobile genetic elements facilitate horizontal gene transfer. However, the prevalence and extent of these phenomena in different taxonomic groups that inhabit wastewater are still not fully understood. Here, we determined the presence of ARGs in metagenome-assembled genomes (MAGs) and evaluated the risks of MAG-carrying ARGs in potential human pathogens. The potential of these ARGs to be transmitted horizontally or vertically was also determined. A total of 5,916 MAGs (completeness >50%, contamination <10%) were recovered, covering 68 phyla and 279 genera. MAGs were dereplicated into 1,204 genome operational taxonomic units (gOTUs) as a proxy for species ( average nucleotide identity >0.95). The dominant ARG classes detected were bacitracin, multi-drug, macrolide-lincosamide-streptogramin (MLS), glycopeptide, and aminoglycoside, and 10.26% of them were located on plasmids. The main hosts of ARGs belonged to Escherichia, Klebsiella, Acinetobacter, Gresbergeria, Mycobacterium, and Thauera. Our data showed that 253 MAGs carried virulence factor genes (VFGs) divided into 44 gOTUs, of which 45 MAGs were carriers of ARGs, indicating that potential human pathogens carried ARGs. Alarmingly, the MAG assigned as Escherichia coli contained 159 VFGs, of which 95 were located on chromosomes and 10 on plasmids. In addition to shedding light on the prevalence of ARGs in individual genomes recovered from activated sludge and wastewater, our study demonstrates a workflow that can identify antimicrobial-resistant pathogens in complex microbial communities. IMPORTANCE: Antimicrobial resistance (AMR) threatens the health of humans, animals, and natural ecosystems. In our study, an analysis of 165 metagenomes from wastewater revealed antibiotic-targeted alteration, efflux, and inactivation as the most prevalent AMR mechanisms. We identified several genera correlated with multiple ARGs, including Klebsiella, Escherichia, Acinetobacter, Nitrospira, Ottowia, Pseudomonas, and Thauera, which could have significant implications for AMR transmission. The abundance of bacA, mexL, and aph(3")-I in the genomes calls for their urgent management in wastewater. Our approach could be applied to different ecosystems to assess the risk of potential pathogens containing ARGs. Our findings highlight the importance of managing AMR in wastewater and can help design measures to reduce the transmission and evolution of AMR in these systems.