PLSDB: a resource of complete bacterial plasmids

Two New Rhizobiales Species Isolated from Root Nodules of Common Sainfoin (Onobrychis viciifolia) Show Different Plant Colonization Strategies

Root nodules of legume plants are primarily inhabited by rhizobial nitrogen-fixing bacteria. Here, we propose two new Rhizobiales species isolated from root nodules of common sainfoin (Onobrychis viciifolia), as shown by core-gene phylogeny, overall genome relatedness indices, and pan-genome analysis. Mesorhizobium onobrychidis sp. nov. actively induces nodules and achieves atmospheric nitrogen and carbon dioxide fixation. This species appears to be depleted in motility genes and is enriched in genes for direct effects on plant growth performance. Its genome reveals functional and plant growth-promoting signatures, like a large unique chromosomal genomic island with high density of symbiotic genetic traits. Onobrychidicola muellerharveyae gen. nov. sp. nov. is described as a type species of the new genus Onobrychidicola in Rhizobiaceae. This species comprises unique genetic features and plant growth-promoting traits (PGPTs), which strongly indicate its function in biotic stress reduction and motility. We applied a newly developed bioinformatics approach for in silico prediction of PGPTs (PGPT-Pred), which supports the different lifestyles of the two new species and the plant growth-promoting performance of M. onobrychidis in the greenhouse trial. IMPORTANCE The intensive use of chemical fertilizers has a variety of negative effects on the environment. Increased utilization of biological nitrogen fixation (BNF) is one way to mitigate those negative impacts. In order to optimize BNF, suitable candidates for different legume species are required. Despite intensive search for new rhizobial bacteria associated with legumes, no new rhizobia have recently been identified from sainfoin (Onobrychis viciifolia). Here, we report on the discovery of two new rhizobial species associated with sainfoin, which are of high importance for the host and may help to increase sustainability in agricultural practices. We employed the combination of in silico prediction and in planta experiments, which is an effective way to detect promising plant growth-promoting bacteria.

Comparative Genomics of Listeria Species Recovered from Meat and Food Processing Facilities

Listeria species (spp.) are contaminants that can survive in food, on equipment, and on food processing premises if appropriate hygiene measures are not used. Homologous stress tolerance genes, virulence gene clusters such as the prfA cluster, and clusters of internalin genes that contribute to the pathogenic potential of the strains can be carried by both pathogenic and nonpathogenic Listeria spp. To enhance understanding of the genome evolution of virulence and virulence-associated properties, a comparative genome approach was used to analyze 41 genome sequences belonging to L. innocua and L. welshimeri isolated from food and food processing facilities. Genetic determinants responsible for disinfectant and stress tolerance were identified, including the efflux cassette bcrABC and Tn6188_qac_1 disinfectant resistance determinant, and stress survival islets. These disinfectant-resistant genes were more frequently found in L. innocua (12%) than in L. welshimeri (2%). Several isolates representing the presumed nonpathogenic L. innocua still carried virulence-associated genes, including LGI2, LGI3, LIPI-3, and LIPI-4 which were absent in all L. welshimeri isolates. The mobile genetic elements identified were plasmids (pLGUG1 and J1776) and prophages (PHAGE_Lister_vB_LmoS_188, PHAGE_Lister_LP_030_3, PHAGE_Lister_A118, PHAGE_Lister_B054, and PHAGE_Lister_vB_LmoS_293). The results suggest that the presumed nonpathogenic isolates especially L. innocua can carry genes relevant to the strain's virulence and stress tolerance in the food and food processing facilities. IMPORTANCE This study provides genomic insights into the recently expanded genus in order to gain valuable information about the evolution of the virulence and stress tolerance properties of the genus Listeria and the distribution of these genetic elements pertinent to the pathogenic potential across Listeria spp. and clonal lineages in South Africa (SA).

Pathogenic strains of Shewanella putrefaciens contain plasmids that are absent in the probiotic strain Pdp11

Shewanella putrefaciens Pdp11 is a strain described as a probiotic for use in aquaculture. However, S. putrefaciens includes strains reported to be pathogenic or saprophytic to fish. Although the probiotic trait has been related to the presence of a group of genes in its genome, the existence of plasmids that could determine the probiotic or pathogenic character of this bacterium is unknown. In the present work, we searched for plasmids in several strains of S. putrefaciens that differ in their pathogenic and probiotic character. Under the different conditions tested, plasmids were only found in two of the five pathogenic strains, but not in the probiotic strain nor in the two saprophytic strains tested. Using a workflow integrating Sanger and Illumina reads, the complete consensus sequences of the plasmids were obtained. Plasmids differed in one ORF and encoded a putative replication initiator protein of the repB family, as well as proteins related to plasmid stability and a toxin-antitoxin system. Phylogenetic analysis showed some similarity to functional repB proteins of other Shewanella species. The implication of these plasmids in the probiotic or pathogenic nature of S. putrefaciens is discussed.

Comparative genome analysis of mycobacteria focusing on tRNA and non-coding RNA

BACKGROUND

The Mycobacterium genus encompasses at least 192 named species, many of which cause severe diseases such as tuberculosis. Non-tuberculosis mycobacteria (NTM) can also infect humans and animals. Some are of emerging concern because they show high resistance to commonly used antibiotics while others are used and evaluated in bioremediation or included in anticancer vaccines.

RESULTS

We provide the genome sequences for 114 mycobacterial type strains and together with 130 available mycobacterial genomes we generated a phylogenetic tree based on 387 core genes and supported by average nucleotide identity (ANI) data. The 244 genome sequences cover most of the species constituting the Mycobacterium genus. The genome sizes ranged from 3.2 to 8.1 Mb with an average of 5.7 Mb, and we identified 14 new plasmids. Moreover, mycobacterial genomes consisted of phage-like sequences ranging between 0 and 4.64% dependent on mycobacteria while the number of IS elements varied between 1 and 290. Our data also revealed that, depending on the mycobacteria, the number of tRNA and non-coding (nc) RNA genes differ and that their positions on the chromosome varied. We identified a conserved core set of 12 ncRNAs, 43 tRNAs and 18 aminoacyl-tRNA synthetases among mycobacteria.

CONCLUSIONS

Phages, IS elements, tRNA and ncRNAs appear to have contributed to the evolution of the Mycobacterium genus where several tRNA and ncRNA genes have been horizontally transferred. On the basis of our phylogenetic analysis, we identified several isolates of unnamed species as new mycobacterial species or strains of known mycobacteria. The predicted number of coding sequences correlates with genome size while the number of tRNA, rRNA and ncRNA genes does not. Together these findings expand our insight into the evolution of the Mycobacterium genus and as such they establish a platform to understand mycobacterial pathogenicity, their evolution, antibiotic resistance/tolerance as well as the function and evolution of ncRNA among mycobacteria.

A low-cost genomics workflow enables isolate screening and strain-level analyses within microbiomes

Earth's environments harbor complex consortia of microbes that affect processes ranging from host health to biogeochemical cycles. Understanding their evolution and function is limited by an inability to isolate genomes in a high-throughput manner. Here, we present a workflow for bacterial whole-genome sequencing using open-source labware and the OpenTrons robotics platform, reducing costs to approximately $10 per genome. We assess genomic diversity within 45 gut bacterial species from wild-living chimpanzees and bonobos. We quantify intraspecific genomic diversity and reveal divergence of homologous plasmids between hosts. This enables population genetic analyses of bacterial strains not currently possible with metagenomic data alone.

Effects of Nutrient Level and Growth Rate on the Conjugation Process That Transfers Mobile Antibiotic Resistance Genes in Continuous Cultures

Bacteria in the effluent of wastewater treatment plants (WWTPs) can transfer antibiotic resistance genes (ARGs) to the bacteria in receiving water through conjugation; however, there is a lack of quantitative assessment of this phenomenon in continuous cultures. Our objective was to determine the effects of background nutrient levels in river water column and growth rates of bacteria on the conjugation frequency of ARGs from effluent bacteria to river bacteria, as well as on the resulting resistance level (i.e., MICs) of the river bacteria. Chemostats were employed to simulate the discharge points of WWTPs into rivers, where effluent bacteria (donor cells) meet river bacteria (recipient cells). Both donor and recipient cells were Escherichia coli cells, and the donor cells were constructed by filter mating with bacteria in the effluent of a local WWTP. Results showed that higher bacterial growth rate (0.45 h^-1 versus 0.15 h^-1) led to higher conjugation frequencies (10^-4 versus 10^-6 transconjugant per recipient). The nutrient level also significantly affected the conjugation frequency, albeit to a lesser extent than the growth rate. The MIC against tetracycline increased from 2 mg/L in the recipient to 64 to 128 mg/L in transconjugants. In comparison, the MIC only increased to as high as 8 mg/L in mutants. Whole-genome sequencing showed that the tet-containing plasmid in both the donor and the transconjugant cells also occur in other fecal bacterial genera. The quantitative information obtained from this study can inform hazard identification related to the proliferation of wastewater-associated ARGs in surface water. IMPORTANCE WWTPs have been regarded as an important hot spot of ARGs. The discharge point of WWTP effluent, where ARGs may be horizontally transferred from bacteria of treated wastewater to bacteria of receiving water, is an important interface between the human-dominated ecosystem and the natural environment. The use of batch cultures in previous studies cannot adequately simulate the nutrient conditions and growth rates in receiving water. In this study, chemostats were employed to simulate the continuous growth of bacteria in receiving water. Furthermore, the experimental setup allowed for separate investigations on the effects of nutrient levels (i.e., simulating background nutrients in river water) and bacterial growth rates on conjugation frequencies and resulting resistance levels. The study generates statistically sound ecological data that can be used to estimate the risk of wastewater-originated ARGs as part of the One Health framework.

Revealing bacteria-phage interactions in human microbiome through the CRISPR-Cas immune systems

The human gut microbiome is composed of a diverse consortium of microorganisms. Relatively little is known about the diversity of the bacteriophage population and their interactions with microbial organisms in the human microbiome. Due to the persistent rivalry between microbial organisms (hosts) and phages (invaders), genetic traces of phages are found in the hosts' CRISPR-Cas adaptive immune system. Mobile genetic elements (MGEs) found in bacteria include genetic material from phage and plasmids, often resultant from invasion events. We developed a computational pipeline (BacMGEnet), which can be used for inference and exploratory analysis of putative interactions between microbial organisms and MGEs (phages and plasmids) and their interaction network. Given a collection of genomes as the input, BacMGEnet utilizes computational tools we have previously developed to characterize CRISPR-Cas systems in the genomes, which are then used to identify putative invaders from publicly available collections of phage/prophage sequences. In addition, BacMGEnet uses a greedy algorithm to summarize identified putative interactions to produce a bacteria-MGE network in a standard network format. Inferred networks can be utilized to assist further examination of the putative interactions and for discovery of interaction patterns. Here we apply the BacMGEnet pipeline to a few collections of genomic/metagenomic datasets to demonstrate its utilities. BacMGEnet revealed a complex interaction network of the Phocaeicola vulgatus pangenome with its phage invaders, and the modularity analysis of the resulted network suggested differential activities of the different P. vulgatus' CRISPR-Cas systems (Type I-C and Type II-C) against some phages. Analysis of the phage-bacteria interaction network of human gut microbiome revealed a mixture of phages with a broad host range (resulting in large modules with many bacteria and phages), and phages with narrow host range. We also showed that BacMGEnet can be used to infer phages that invade bacteria and their interactions in wound microbiome. We anticipate that BacMGEnet will become an important tool for studying the interactions between bacteria and their invaders for microbiome research.

Nodules of wild legumes as unique natural hotspots of antibiotic resistance genes

Root nodules (RN) of legumes have distinct microenvironment from their symbiotic roots and surrounding soils. The rhizobia can withstand the host-produced phytoalexins and antimicrobial compounds. We thus hypothesize that the wild legume RN may develop unique natural resistome and be antibiotic resistance gene (ARG) hotspots. In this study, in comparison with rhizosphere soil (RS) and bulk soil (BS), we characterized the feature of antibiotic resistance in the RN of two wild legumes, Medicago polymorpha and Astragalus sinicus, by metagenomics. It was shown that the total relative abundance of ARGs followed the order of RN > RS > BS for both legumes. ARGs encoding antibiotic efflux pump predominated in all samples with increased proportion from BS to RN samples for both legumes. Totally 275 ARG subtypes were detected, and diversity of ARGs in RN was significantly lower than in BS samples for both legumes. 32 and 25 unique ARGs subtypes were detected in RN of both legumes. Bacterial community played a key role in shaping nodule-associated resistome because both ARG profiles and bacterial community differed greatly among BS, RS and RN. Rhizobia potentially hosted 10 and 15 ARGs subtypes for both legumes. The number and proportion of plasmid- and ARG-carrying contigs (ACCs) were higher in RN than in BS. Host tracking analysis of plasmid-ACCs suggests that proportion of rhizobial bacteria identified as their hosts decreased from BS to RN samples. No plasmid-ACCs with multiple ARGs were observed in BS samples, whereas they were detected in RN samples of both legumes. Our study showed that even wild legume nodules are unique natural ARG hotspots and enough attention should be paid to the dissemination risk of ARGs posed by globally produced legume crops.

Long-term ecological and evolutionary dynamics in the gut microbiomes of carbapenemase-producing Enterobacteriaceae colonized subjects

Long-term colonization of the gut microbiome by carbapenemase-producing Enterobacteriaceae (CPE) is a growing area of public health concern as it can lead to community transmission and rapid increase in cases of life-threatening CPE infections. Here, leveraging the observation that many subjects are decolonized without interventions within a year, we used longitudinal shotgun metagenomics (up to 12 timepoints) for detailed characterization of ecological and evolutionary dynamics in the gut microbiome of a cohort of CPE-colonized subjects and family members (n = 46; 361 samples). Subjects who underwent decolonization exhibited a distinct ecological shift marked by recovery of microbial diversity, key commensals and anti-inflammatory pathways. In addition, colonization was marked by elevated but unstable Enterobacteriaceae abundances, which exhibited distinct strain-level dynamics for different species (Escherichia coli and Klebsiella pneumoniae). Finally, comparative analysis with whole-genome sequencing data from CPE isolates (n = 159) helped identify substrain variation in key functional genes and the presence of highly similar E. coli and K. pneumoniae strains with variable resistance profiles and plasmid sharing. These results provide an enhanced view into how colonization by multi-drug-resistant bacteria associates with altered gut ecology and can enable transfer of resistance genes, even in the absence of overt infection and antibiotic usage.

Assessing the genomic composition, putative ecological relevance and biotechnological potential of plasmids from sponge bacterial symbionts

Plasmid-mediated transfer of genes can have direct consequences in several biological processes within sponge microbial communities. However, very few studies have attempted genomic and functional characterization of plasmids from marine host-associated microbial communities in general and those of sponges in particular. In the present study, we used an endogenous plasmid isolation method to obtain plasmids from bacterial symbionts of the marine sponges Stylissa carteri and Paratetilla sp. and investigated the genomic composition, putative ecological relevance and biotechnological potential of these plasmids. In total, we isolated and characterized three complete plasmids, three plasmid prophages and one incomplete plasmid. Our results highlight the importance of plasmids to transfer relevant genetic traits putatively involved in microbial symbiont adaptation and host-microbe and microbe-microbe interactions. For example, putative genes involved in bacterial response to chemical stress, competition, metabolic versatility and mediation of bacterial colonization and pathogenicity were detected. Genes coding for enzymes and toxins of biotechnological potential were also detected. Most plasmid prophage coding sequences were, however, hypothetical proteins with unknown functions. Overall, this study highlights the ecological relevance of plasmids in the marine sponge microbiome and provides evidence that plasmids of sponge bacterial symbionts may represent an untapped resource of genes of biotechnological interest.

Methicillin Resistant Staphylococci Isolated from Goats and Their Farm Environments in Saudi Arabia Genotypically Linked to Known Human Clinical Isolates: a Pilot Study

We conducted a pilot whole genome sequencing (WGS) study to characterize the genotypes of nine methicillin resistant staphylococci (MRS) isolates recovered from goats and their farm environments in Eastern Province, Saudi Arabia, between November 2019 to August 2020. Seven out of nine isolates were methicillin resistant Staphylococcus aureus (MRSA), and two were methicillin resistant Staphylococcus epidermidis (MRSE). All MRSA isolates possessed genotypes previously identified to infect humans, including isolates harboring ST6-SCCmec IV-t304 (n = 4), ST5-SCCmec VI- t688 (n = 2) and ST5-SCCmec V-t311 (n = 1). 2 MRSA isolates possessed plasmids that were genetically similar to those identified in S. aureus isolates recovered from humans and poultry. In contrast, plasmids found in three MRSA isolates and one MRSE isolate were genetically similar to those recovered from humans. All MRSA isolates harbored the host innate modulate genes sak and scn previously associated with human infections. The genotypes of MRSE isolates were determined as ST35, a well-known zoonotic sequence type and ST153, which has been associated with humans. However, the MRSE isolates were untypeable due to extra ccr complexes identified in their SCCmec elements. Moreover, we identified in ST153 isolate SCCmec element also harbored the Arginine Catabolic Mobile Element (ACME) IV. All MRS isolates were phenotypically resistant to trimethoprim-sulfamethoxazole, an antibiotic for the decolonization of MRS. Three isolates carried antibiotic resistance genes in their SCCmec elements that were not previously described, including those encoding fusidic acid resistance (fusC) and trimethoprim resistance (dfrC) incorporated in the MRSA SCCmec VI. IMPORTANCE Our findings demonstrate a possible cross-transmission of methicillin resistant staphylococci between goats and their local environments and between goats and humans. Due to ever increasing resistance to multiple antibiotics, the burden of MRS has a significant impact on livestock farming, public health, and the economy worldwide. This study highlights that implementing a holistic approach to whole genome sequencing surveillance in livestock and farm environments would aid our understanding of the transmission of methicillin resistant staphylococci and, most importantly, allow us to implement appropriate infection control and hygiene practices.

Extensive diversity and rapid turnover of phage defense repertoires in cheese-associated bacterial communities

BACKGROUND

Phages are key drivers of genomic diversity in bacterial populations as they impose strong selective pressure on the evolution of bacterial defense mechanisms across closely related strains. The pan-immunity model suggests that such diversity is maintained because the effective immune system of a bacterial species is the one distributed across all strains present in the community. However, only few studies have analyzed the distribution of bacterial defense systems at the community-level, mostly focusing on CRISPR and comparing samples from complex environments. Here, we studied 2778 bacterial genomes and 188 metagenomes from cheese-associated communities, which are dominated by a few bacterial taxa and occur in relatively stable environments.

RESULTS

We corroborate previous laboratory findings that in cheese-associated communities nearly identical strains contain diverse and highly variable arsenals of innate and adaptive (i.e., CRISPR-Cas) immunity systems suggesting rapid turnover. CRISPR spacer abundance correlated with the abundance of matching target sequences across the metagenomes providing evidence that the identified defense repertoires are functional and under selection. While these characteristics align with the pan-immunity model, the detected CRISPR spacers only covered a subset of the phages previously identified in cheese, providing evidence that CRISPR does not enable complete immunity against all phages, and that the innate immune mechanisms may have complementary roles.

CONCLUSIONS

Our findings show that the evolution of bacterial defense mechanisms is a highly dynamic process and highlight that experimentally tractable, low complexity communities such as those found in cheese, can help to understand ecological and molecular processes underlying phage-defense system relationships. These findings can have implications for the design of robust synthetic communities used in biotechnology and the food industry. Video Abstract.

Diversity and dynamics of the CRISPR-Cas systems associated with Bacteroides fragilis in human population

Background

CRISPR-Cas (clustered regularly interspaced short palindromic repeats—CRISPR-associated proteins) systems are adaptive immune systems commonly found in prokaryotes that provide sequence-specific defense against invading mobile genetic elements (MGEs). The memory of these immunological encounters are stored in CRISPR arrays, where spacer sequences record the identity and history of past invaders. Analyzing such CRISPR arrays provide insights into the dynamics of CRISPR-Cas systems and the adaptation of their host bacteria to rapidly changing environments such as the human gut.

Results

In this study, we utilized 601 publicly available Bacteroides fragilis genome isolates from 12 healthy individuals, 6 of which include longitudinal observations, and 222 available B. fragilis reference genomes to update the understanding of B. fragilis CRISPR-Cas dynamics and their differential activities. Analysis of longitudinal genomic data showed that some CRISPR array structures remained relatively stable over time whereas others involved radical spacer acquisition during some periods, and diverse CRISPR arrays (associated with multiple isolates) co-existed in the same individuals with some persisted over time. Furthermore, features of CRISPR adaptation, evolution, and microdynamics were highlighted through an analysis of host-MGE network, such as modules of multiple MGEs and hosts, reflecting complex interactions between B. fragilis and its invaders mediated through the CRISPR-Cas systems.

Conclusions

We made available of all annotated CRISPR-Cas systems and their target MGEs, and their interaction network as a web resource at https://omics.informatics.indiana.edu/CRISPRone/Bfragilis. We anticipate it will become an important resource for studying of B. fragilis, its CRISPR-Cas systems, and its interaction with mobile genetic elements providing insights into evolutionary dynamics that may shape the species virulence and lead to its pathogenicity.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08770-8).

Antibiotic resistance genes in the gut microbiota of mothers and linked neonates with or without sepsis from low- and middle-income countries

Early development of the microbiome has been shown to affect general health and physical development of the infant and, although some studies have been undertaken in high-income countries, there are few studies from low- and middle-income countries. As part of the BARNARDS study, we examined the rectal microbiota of 2,931 neonates (term used up to 60 d) with clinical signs of sepsis and of 15,217 mothers screening for bla_CTX-M-15, bla_NDM, bla_KPC and bla_OXA-48-like genes, which were detected in 56.1%, 18.5%, 0% and 4.1% of neonates’ rectal swabs and 47.1%, 4.6%, 0% and 1.6% of mothers’ rectal swabs, respectively. Carbapenemase-positive bacteria were identified by MALDI-TOF MS and showed a high diversity of bacterial species (57 distinct species/genera) which exhibited resistance to most of the antibiotics tested. Escherichia coli, Klebsiella pneumoniae and Enterobacter cloacae/E. cloacae complex, the most commonly found isolates, were subjected to whole-genome sequencing analysis and revealed close relationships between isolates from different samples, suggesting transmission of bacteria between neonates, and between neonates and mothers. Associations between the carriage of antimicrobial resistance genes (ARGs) and healthcare/environmental factors were identified, and the presence of ARGs was a predictor of neonatal sepsis and adverse birth outcomes.

Analysis of gut microbiota of mothers and its neonates—as part of the BARNARDS study—reveals associations between β-lactamase gene carriage and neonatal sepsis risk in low-income settings.

Shooting hoops: globetrotting plasmids spreading more than just antimicrobial resistance genes across One Health

Our study provides novel insights into the global nature of antimicrobial resistance (AMR) plasmids across the food chain. We provide compelling evidence of the globetrotting nature of AMR plasmids and the need for surveillance to sequence plasmids with a template of analyses for others to expand these data. The AMR plasmids analysed were detected in 63 countries and in samples from humans, animals and the environment. They contained a combination of known and novel AMR genes, metal resistance genes, virulence factors, phage and replicon types.

Draft-genome analysis provides insights into the virulence properties and genome plasticity of Vibrio fluvialis organisms isolated from shrimp farms and Turag river in Bangladesh

Vibrio fluvialis is an opportunistic waterborne and seafood-borne enteric pathogen capable of causing severe diarrhea leading to death. This pathogen is endemic to Bangladesh, a country which is a major producer of cultured shrimp and wild-caught prawns. In this study, we carried out whole-genome sequencing of three V. fluvialis organisms isolated from shrimp farm and river sediment showing strong pathogenic characteristics in vivo and in vitro and compared their genomes against other V. fluvialis and related pathogenic species to glean insights into their potential as pathogens. Numerous virulence-associated genes including hemolysins, cytolysins, three separate Type IV pili, Types II and VI secretion systems, biofilm, and the V. cholerae pathogenesis regulating gene, toxR, were identified. Moreover, we found strain S-10 to have the propensity to acquire antibiotic resistance genes through horizontal gene transfer. These findings indicate that shrimp farms and rivers could be potential sources of V. fluvialis organisms which are an infection threat of public health concern.

Revealing Genomic Insights of the Unexplored Porcine Pathogen Actinobacillus pleuropneumoniae Using Whole Genome Sequencing

Actinobacillus pleuropneumoniae (APP) is the causative agent of pleuropneumonia in pigs, one of the most relevant bacterial respiratory diseases in the swine industry. To date, 19 serotypes have been described based on capsular polysaccharide typing with significant virulence dissimilarities. In this study, 16 APP isolates from Spanish origin were selected to perform antimicrobial susceptibility tests and comparative genomic analysis using whole genome sequencing (WGS). To obtain a more comprehensive worldwide molecular epidemiologic analyses, all APP whole genome assemblies available at the National Center for Biotechnology Information (NCBI) at the time of the study were also included. An in-house in silico PCR approach enabled the correct serotyping of unserotyped or incorrectly serotyped isolates and allowed for the discrimination between serotypes 9 and 11. A pangenome analysis identified the presence or absence of gene clusters to be serotype specific, as well as virulence profile analyses targeting the apx operons. Antimicrobial resistance genes were correlated to the presence of specific plasmids. Altogether, this study provides new insights into the genetic variability within APP serotypes, correlates phenotypic tests with bioinformatic analyses and manifests the benefits of populated databases for a better assessment of diversity and variability of relatively unknown pathogens. Overall, genomic comparative analysis enhances the understanding of transmission and epidemiological patterns of this species and suggests vertical transmission of the pathogen, including the resistance genes, within the Spanish integrated systems.

IMPORTANCE Pleuropneumonia is one of the most relevant respiratory infections in the swine industry. Despite Actinobacillus pleuropneumoniae (APP) being one of the most important pathogens in the pig production, this is the first comparative study including all available whole genome sequencing data from NCBI. Moreover, this study also includes 16 APP isolates of Spanish origin with known epidemiological relationships through vertical integrated systems. Genomic comparisons provided a deeper understanding of molecular and epidemiological knowledge between different APP serotypes. Furthermore, determination of resistance and toxin profiles allowed correlation with the presence of mobile genetic elements and specific serotype, respectively.

Genomic Analysis of an I1 Plasmid Hosting a sul3-Class 1 Integron and blaSHV-12 within an Unusual Escherichia coli ST297 from Urban Wildlife

Wild birds, particularly silver gulls (Chroicocephalus novaehollandiae) that nest near anthropogenic sites, often harbour bacteria resistant to multiple antibiotics, including those considered of clinical importance. Here, we describe the whole genome sequence of Escherichia coli isolate CE1867 from a silver gull chick sampled in 2012 that hosted an I1 pST25 plasmid with bla_SHV-12, a β-lactamase gene that encodes the ability to hydrolyze oxyimino β-lactams, and other antibiotic resistance genes. Isolate CE1867 is an ST297 isolate, a phylogroup B1 lineage, and clustered with a large ST297 O130:H11 clade, which carry Shiga toxin genes. The I1 plasmid belongs to plasmid sequence type 25 and is notable for its carriage of an atypical sul3-class 1 integron with mefB_∆260, a structure most frequently reported in Australia from swine. This integron is a typical example of a Tn21-derived element that captured sul3 in place of the standard sul1 structure. Interestingly, the mercury resistance (mer) module of Tn21 is missing and has been replaced with Tn2-bla_TEM-1 and a bla_SHV-12 encoding module flanked by direct copies of IS26. Comparisons to similar plasmids, however, demonstrate a closely related family of ARG-carrying plasmids that all host variants of the sul3-associated integron with conserved Tn21 insertion points and a variable presence of both mer and mefB truncations, but predominantly mefB_∆260.

UG/Abi: a highly diverse family of prokaryotic reverse transcriptases associated with defense functions

Reverse transcriptases (RTs) are enzymes capable of synthesizing DNA using RNA as a template. Within the last few years, a burst of research has led to the discovery of novel prokaryotic RTs with diverse antiviral properties, such as DRTs (Defense-associated RTs), which belong to the so-called group of unknown RTs (UG) and are closely related to the Abortive Infection system (Abi) RTs. In this work, we performed a systematic analysis of UG and Abi RTs, increasing the number of UG/Abi members up to 42 highly diverse groups, most of which are predicted to be functionally associated with other gene(s) or domain(s). Based on this information, we classified these systems into three major classes. In addition, we reveal that most of these groups are associated with defense functions and/or mobile genetic elements, and demonstrate the antiphage role of four novel groups. Besides, we highlight the presence of one of these systems in novel families of human gut viruses infecting members of the Bacteroidetes and Firmicutes phyla. This work lays the foundation for a comprehensive and unified understanding of these highly diverse RTs with enormous biotechnological potential.

High Prevalence of bla CTXM-1/IncI1-Iγ/ST3 Plasmids in Extended-Spectrum β-Lactamase-Producing Escherichia coli Isolates Collected From Domestic Animals in Guadeloupe (French West Indies)

Extended-spectrum β-lactamase-producing Enterobacteriaceae (ESBL-E) have been classified in the group of resistant bacteria of highest priority. We determined the prevalence of ESBL-E collected in feces from household and shelter pets in Guadeloupe (French West Indies). A single rectal swab was taken from 125 dogs and 60 cats between June and September 2019. The prevalence of fecal carriage of ESBL-E was 7.6% (14/185, 95% CI: 4.2-12.4), within the range observed worldwide. The only risk factor associated with a higher prevalence of ESBL-E rectal carriage was a stay in a shelter, suggesting that refuges could be hotspots for their acquisition. All but one (Klebsiella pneumoniae from a cat) were Escherichia coli. We noted the presence of a bla _CTX-M-1/IncI1-Iγ/sequence type (ST3) plasmid in 11 ESBL-producing E. coli isolates belonging to ST328 (n = 6), ST155 (n = 4) and ST953 (n = 1). A bla _CTX-M-15 gene was identified in the three remaining ESBL-E isolates. The bla _CTX-M-1 and most of the antimicrobial resistance genes were present in a well-conserved large conjugative IncI1-Iγ/ST3 plasmid characterized by two accessory regions containing antibiotic resistance genes. The plasmid has been detected worldwide in E. coli isolates from humans and several animal species, such as food-producing animals, wild birds and pets, and from the environment. This study shows the potential role of pets as a reservoir of antimicrobial-resistant bacteria or genes for humans and underlines the importance of basic hygiene measures by owners of companion animals.

Virulence plasmid pINV as a genetic signature for Shigella flexneri phylogeny

Shigella flexneri is a major health burden in low- and middle-income countries, where it is a leading cause of mortality associated with diarrhoea in children, and shows an increasing incidence among travellers and men having sex with men. Like all Shigella spp., S. flexneri has evolved from commensal Escherichia coli following the acquisition of a large plasmid pINV, which contains genes essential for virulence. Current sequence typing schemes of Shigella are based on combinations of chromosomal genetic loci, since pINV-encoded virulence genes are often lost during growth in the laboratory, making these elements inappropriate for sequence typing. By performing comparative analysis of pINVs from S. flexneri strains isolated from different geographical regions and belonging to different serotypes, we found that in contrast to plasmid-encoded virulence genes, plasmid maintenance genes are highly stable pINV-encoded elements. For the first time, to our knowledge, we have developed a S. flexneri plasmid multilocus sequence typing (pMLST) method based on different combinations of alleles of the vapBC and yacAB toxin–antitoxin (TA) systems, and the parAB partitioning system. This enables typing of S. flexneri pINV plasmids into distinct ‘virulence sequence types’ (vSTs). Furthermore, the phylogenies of vST alleles and bacterial host core genomes suggests an intimate co-evolution of pINV with the chromosome of its bacterial host, consistent with previous findings. This work demonstrates the potential of plasmid maintenance loci as genetic characteristics to study as well as to trace the molecular phylogenesis of S. flexneri pINV and the phylogenetic relationship of this plasmid with its bacterial host.

Phenotypic characterization and analysis of complete genomes of two distinct strains of the proposed species "L. swaminathanii"

Recently, a new Listeria species, “Listeria swaminathanii”, was proposed. Here, we phenotypically and genotypically characterize two additional strains that were previously obtained from soil samples and compare the results to the type strain. Complete genomes for both strains were assembled from hybrid Illumina and Nanopore sequencing reads and annotated. Further genomic analysis including average nucleotide identity (ANI) and detection of mobile genetic elements and genes of interest (e.g., virulence-associated) were conducted. The strains showed 98.7–98.8% ANI with the type strain. The UTK C1-0015 genome contained a partial monocin locus and a plasmid, while the UTK C1-0024 genome contained a full monocin locus and a prophage. Phenotypic characterization consistent with those performed on the proposed type strain was conducted to assess consistency of phenotypes across a greater diversity of the proposed species (n = 3 instead of n = 1). Only a few findings were notably different from those of the type strain, such as catalase activity, glycerol metabolism, starch metabolism, and growth at 41 °C. This study further expands our understanding of this newly proposed sensu stricto Listeria species.

CRISPR-Cas systems are widespread accessory elements across bacterial and archaeal plasmids

Many prokaryotes encode CRISPR-Cas systems as immune protection against mobile genetic elements (MGEs), yet a number of MGEs also harbor CRISPR-Cas components. With a few exceptions, CRISPR-Cas loci encoded on MGEs are uncharted and a comprehensive analysis of their distribution, prevalence, diversity, and function is lacking. Here, we systematically investigated CRISPR-Cas loci across the largest curated collection of natural bacterial and archaeal plasmids. CRISPR-Cas loci are widely but heterogeneously distributed across plasmids and, in comparison to host chromosomes, their mean prevalence per Mbp is higher and their distribution is distinct. Furthermore, the spacer content of plasmid CRISPRs exhibits a strong targeting bias towards other plasmids, while chromosomal arrays are enriched with virus-targeting spacers. These contrasting targeting preferences highlight the genetic independence of plasmids and suggest a major role for mediating plasmid-plasmid conflicts. Altogether, CRISPR-Cas are frequent accessory components of many plasmids, which is an overlooked phenomenon that possibly facilitates their dissemination across microbiomes.

Whole genome sequencing of Klebsiella pneumoniae clinical isolates sequence type 627 isolated from Egyptian patients

Klebsiella pneumoniae is considered a threat to public health especially due to multidrug resistance emergence. It is largely oligoclonal based on multi-locus sequence typing (MLST); in Egypt, ST 627 was recently detected. Despites the global dissemination of this ST, there is still paucity of information about it. Herein, we used 4 K. pneumoniae ST627 for whole genome sequencing utilizing an Illumina MiSeq platform. Genome sequences were examined for resistance and virulence determinants, capsular types, plasmids, insertion sequences, phage regions, and Clustered Regularly Interspaced Palindromic Repeats (CRISPR) regions using bioinformatic analysis. The molecular characterization revealed 15 and 65 antimicrobial resistance and virulence genes, respectively. Resistance genes such as tet(D), aph(3’’)-Ib, aph(6)-Id, blaTEM-234, fosA, and fosA6; were mainly responsible for tetracycline, aminoglycoside, and fosfomycin resistance; respectively. The capsular typing revealed that the four strains are KL-24 and O1v1. One plasmid was found in all samples known as pC17KP0052-1 and another plasmid with accession no. NZ_CP032191.1 was found only in K90. IncFIB(K) and IncFII(K) are two replicons found in all samples, while ColRNAI replicon was found only in K90. Entero P88, Salmon SEN5, and Klebsi phiKO2 intact phage regions were identified. All samples harbored CRISPR arrays including CRISPR1 and CRISPR2. Our results shed light on critical tasks of mobile genetic elements in ST 627 in antibiotic resistance spreading.

Emergence of Neonatal Sepsis Caused by MCR-9- and NDM-1-Co-Producing Enterobacter hormaechei in China

Mobile colistin resistance (mcr) genes represent an emerging threat to public health. Reports on the prevalence, antimicrobial profiles, and clonality of MCR-9-producing Enterobacter cloacae complex (ECC) isolates on a national scale in China are limited. We screened 3,373 samples from humans, animals, and the environment and identified eleven MCR-9-positive ECC isolates. We further investigated their susceptibility, epidemiology, plasmid profiles, genetic features, and virulence potential. Ten strains were isolated from severe bloodstream infection cases, especially three of them were recovered from neonatal sepsis. Enterobacter hormaechei was the most predominant species among the MCR-9-producing ECC population. Moreover, the co-existence of MCR-9, CTX-M, and SHV-12 encoding genes in MCR-9-positive isolates was globally observed. Notably, mcr-9 was mainly carried by IncHI2 plasmids, and we found a novel ~187 kb IncFII plasmid harboring mcr-9, with low similarity with known plasmids. In summary, our study presented genomic insights into genetic characteristics of MCR-9-producing ECC isolates retrieved from human, animal, and environment samples with one health perspective. This study is the first to reveal NDM-1- and MCR-9-co-producing ECC from neonatal sepsis in China. Our data highlights the risk for the hidden spread of the mcr-9 colistin resistance gene.

Metagenomics Insights Into the Microbial Diversity and Microbiome Network Analysis on the Heterogeneity of Influent to Effluent Water

Sanitizing the water sources of local communities is important to control the spread of microbial resistance genes, especially those for water-borne illnesses. The activities of antibiotic resistance gene (ARG)-host pathogens pose a threat to public health, and it has been estimated that the infection will lead up to 10 million deaths globally by the year 2050. Hence, in this study, we aim to analyze the efficiency of our municipal wastewater treatment plant (WWTP) process in producing pathogen-free water by investigating the microbial composition between influent and effluent water sites. Shotgun metagenomics sequencing using the Illumina platform was performed on the influent and effluent samples of six different WWTP sites located in Johore, Malaysia. After raw data pre-processing, the non-redundant contigs library was then aligned against BLASTP for taxonomy profiling and the Comprehensive Antibiotic Resistance Database for ARG annotation. Interestingly, the alpha-diversity result reported that effluent site samples showed higher abundance and diverse heterogeneity compared to the influent site. The principal component analysis (PCA) and non-metric multidimensional scaling (NMDS) plots also suggested that effluent sites showed high variation in the genetic material due to loosely clustered sample plots, as compared to the tightly clustered influent samples. This study has successfully identified the top three abundant phyla in influent-Proteobacteria, Firmicutes, and Bacteroidetes-and effluent-Proteobacteria, Actinobacteria, and Bacteroidetes-water. Despite the overlap within the top three abundant phyla in influent and effluent sites (Proteobacteria and Bacteroidetes), the ARG composition heat map and drug class phenotype plot bar exhibits a general trend of a downward shift, showing the efficiency of WWTP in reducing opportunistic pathogens. Overall, it was demonstrated that our municipal WWTP efficiently eliminated pathogenic microbes from the influent water before its total discharge to the environment, though not with the total elimination of microorganisms. This metagenomics study allowed for an examination of our water source and showed the potential interaction of species and ARGs residing in the influent and effluent environment. Both microbial profile structure and co-occurrence network analysis provide integrated understanding regarding the diversity of microorganisms and interactions for future advanced water sanitation treatments.

In-depth genome analysis of Bacillus sp. BH32, a salt stress-tolerant endophyte obtained from a halophyte in a semiarid region

Endophytic strains belonging to the Bacillus cereus group were isolated from the halophytes Atriplex halimus L. (Amaranthaceae) and Tamarix aphylla L. (Tamaricaceae) from costal and continental regions in Algeria. Based on their salt tolerance (up to 5%), the strains were tested for their ability to alleviate salt stress in tomato and wheat. Bacillus sp. strain BH32 showed the highest potential to reduce salinity stress (up to + 50% and + 58% of dry weight improvement, in tomato and wheat, respectively, compared to the control). To determine putative mechanisms involved in salt tolerance and plant growth promotion, the whole genome of Bacillus sp. BH32 was sequenced, annotated, and used for comparative genomics against the genomes of closely related strains. The pangenome of Bacillus sp. BH32 and its closest relative was further analyzed. The phylogenomic analyses confirmed its taxonomic position, a member of the Bacillus cereus group, with intergenomic distances (GBDP analysis) pinpointing to a new taxon (digital DNA-DNA hybridization, dDDH < 70%). Genome mining unveiled several genes involved in stress tolerance, production of anti-oxidants and genes involved in plant growth promotion as well as in the production of secondary metabolites. KEY POINTS : • Bacillus sp. BH32 and other bacterial endophytes were isolated from halophytes, to be tested on tomato and wheat and to limit salt stress adverse effects. • The strain with the highest potential was then studied at the genomic level to highlight numerous genes linked to plant growth promotion and stress tolerance. • Pangenome approaches suggest that the strain belongs to a new taxon within the Bacillus cereus group.

Use of Nanopore Sequencing to Characterise the Genomic Architecture of Mobile Genetic Elements Encoding bla CTX-M-15 in Escherichia coli Causing Travellers' Diarrhoea

Increasing levels of antimicrobial resistance (AMR) have been documented in Escherichia coli causing travellers’ diarrhoea, particularly to the third-generation cephalosporins. Diarrhoeagenic E. coli (DEC) can act as a reservoir for the exchange of AMR genes between bacteria residing in the human gut, enabling them to survive and flourish through the selective pressures of antibiotic treatments. Using Oxford Nanopore Technology (ONT), we sequenced eight isolates of DEC from four patients’ specimens who had all recently returned to the United Kingdome from Pakistan. Sequencing yielded two DEC harbouring bla_CTX-M-15 per patient, all with different sequence types (ST) and belonging to five different pathotypes. The study aimed to determine whether bla_CTX-M-15 was located on the chromosome or plasmid and to characterise the drug-resistant regions to better understand the mechanisms of onward transmission of AMR determinants. Patients A and C both had one isolate where bla_CTX-M-15 was located on the plasmid (899037 & 623213, respectively) and one chromosomally encoded (899091 & 623214, respectively). In patient B, bla_CTX-M-15 was plasmid-encoded in both DEC isolates (786605 & 7883090), whereas in patient D, bla_CTX-M-15 was located on the chromosome in both DEC isolates (542093 & 542099). The two bla_CTX-M-15-encoding plasmids associated with patient B were different although the bla_CTX-M-15-encoding plasmid isolated from 788309 (IncFIB) exhibited high nucleotide similarity to the bla_CTX-M-15-encoding plasmid isolated from 899037 (patient A). In the four isolates where bla_CTX-M-15 was chromosomally encoded, two isolates (899091 & 542099) shared the same insertion site. The bla_CTX-M-15 insertion site in isolate 623214 was described previously, whereas that of isolate 542093 was unique to this study. Analysis of Nanopore sequencing data enables us to characterise the genomic architecture of mobile genetic elements encoding AMR determinants. These data may contribute to a better understanding of persistence and onward transmission of AMR determinants in multidrug-resistant (MDR) E. coli causing gastrointestinal and extra-intestinal infections.

Precise genotyping of circular mobile elements from metagenomic data uncovers human-associated plasmids with recent common ancestors

Mobile genetic elements with circular genomes play a key role in the evolution of microbial communities. Their circular genomes correspond to circular walks in metagenome graphs, and yet, assemblies derived from natural microbial communities produce graphs riddled with spurious cycles, complicating the accurate reconstruction of circular genomes. We present DomCycle, an algorithm that reconstructs likely circular genomes based on the identification of so-called 'dominant' graph cycles. In the implementation we leverage paired reads to bridge assembly gaps and scrutinize cycles through a nucleotide-level analysis, making the approach robust to misassembly artifacts. We validated the approach using simulated and real sequencing data. Application of DomCycle to 32 publicly available DNA shotgun sequence data sets from diverse natural environments led to the reconstruction of hundreds of circular mobile genomes. Clustering revealed 20 highly prevalent and cryptic plasmids that have clonal population structures with recent common ancestors. This method facilitates the study of microbial communities that evolve through horizontal gene transfer.

Comparative Genomic Analysis of a Panton–Valentine Leukocidin-Positive ST22 Community-Acquired Methicillin-Resistant Staphylococcus aureus from Pakistan

Staphylococcus aureus (S. aureus) ST22 is considered a clinically important clone because an epidemic strain EMRSA-15 belongs to ST22, and several outbreaks of this clone have been documented worldwide. We performed genomic analysis of an S. aureus strain Lr2 ST22 from Pakistan and determined comparative analysis with other ST22 strains. The genomic data show that Lr2 belongs to spa-type t2986 and harbors staphylococcal cassette chromosome mec (SCCmec) type IVa(2B), one complete plasmid, and seven prophages or prophage-like elements. The strain harbors several prophage-associated virulence factors, including Panton–Valentine leukocidin (PVL) and toxic shock syndrome toxin (TSST). The single nucleotide polymorphism (SNPs)-based phylogenetic relationship inferred from whole genome and core genome revealed that strain Lr2 exhibits the nearest identities to a South African community-acquired methicillin-resistant S. aureus (CA-MRSA) ST22 strain and makes a separate clade with an Indian CA-MRSA ST22 strain. Although most ST22 strains carry blaZ, mecA, and mutations in gyrA, the Lr2 strain does not have the blaZ gene but, unlike other ST22 strains, carries the antibiotic resistance genes erm(C) and aac(6′)-Ie-aph(2″)-Ia. Among ST22 strains analyzed, only the strain Lr2 possesses both PVL and TSST genes. The functional annotation of genes unique to Lr2 revealed that mobilome is the third-largest Cluster of Orthologous Genes (COGs) category, which encodes genes associated with prophages and transposons. This possibly makes methicillin-resistant S. aureus (MRSA) Lr2 ST22 strain highly virulent, and this study would improve the knowledge of MRSA ST22 strains in Pakistan. However, further studies are needed on a large collection of MRSA to comprehend the genomic epidemiology and evolution of this clone in Pakistan.

Ultrafast and Cost-Effective Pathogen Identification and Resistance Gene Detection in a Clinical Setting Using Nanopore Flongle Sequencing

Rapid bacterial identification and antimicrobial resistance gene (ARG) detection are crucial for fast optimization of antibiotic treatment, especially for septic patients where each hour of delayed antibiotic prescription might have lethal consequences. This work investigates whether the Oxford Nanopore Technology’s (ONT) Flongle sequencing platform is suitable for real-time sequencing directly from blood cultures to identify bacteria and detect resistance-encoding genes. For the analysis, we used pure bacterial cultures of four clinical isolates of Escherichia coli and Klebsiella pneumoniae and two blood samples spiked with either E. coli or K. pneumoniae that had been cultured overnight. We sequenced both the whole genome and plasmids isolated from these bacteria using two different sequencing kits. Generally, Flongle data allow rapid bacterial ID and resistome detection based on the first 1,000–3,000 generated sequences (10 min to 3 h from the sequencing start), albeit ARG variant identification did not always correspond to ONT MinION and Illumina sequencing-based data. Flongle data are sufficient for 99.9% genome coverage within at most 20,000 (clinical isolates) or 50,000 (positive blood cultures) sequences generated. The SQK-LSK110 Ligation kit resulted in higher genome coverage and more accurate bacterial identification than the SQK-RBK004 Rapid Barcode kit.

The Dynamics of the Antimicrobial Resistance Mobilome of Salmonella enterica and Related Enteric Bacteria

The foodborne pathogen Salmonella enterica is considered a global public health risk. Salmonella enterica isolates can develop resistance to several antimicrobial drugs due to the rapid spread of antimicrobial resistance (AMR) genes, thus increasing the impact on hospitalization and treatment costs, as well as the healthcare system. Mobile genetic elements (MGEs) play key roles in the dissemination of AMR genes in S. enterica isolates. Multiple phenotypic and molecular techniques have been utilized to better understand the biology and epidemiology of plasmids including DNA sequence analyses, whole genome sequencing (WGS), incompatibility typing, and conjugation studies of plasmids from S. enterica and related species. Focusing on the dynamics of AMR genes is critical for identification and verification of emerging multidrug resistance. The aim of this review is to highlight the updated knowledge of AMR genes in the mobilome of Salmonella and related enteric bacteria. The mobilome is a term defined as all MGEs, including plasmids, transposons, insertion sequences (ISs), gene cassettes, integrons, and resistance islands, that contribute to the potential spread of genes in an organism, including S. enterica isolates and related species, which are the focus of this review.

Clustering biological sequences with dynamic sequence similarity threshold

BACKGROUND

Biological sequence clustering is a complicated data clustering problem owing to the high computation costs incurred for pairwise sequence distance calculations through sequence alignments, as well as difficulties in determining parameters for deriving robust clusters. While current approaches are successful in reducing the number of sequence alignments performed, the generated clusters are based on a single sequence identity threshold applied to every cluster. Poor choices of this identity threshold would thus lead to low quality clusters. There is however little support provided to users in selecting thresholds that are well matched with the input sequences.

RESULTS

We present a novel sequence clustering approach called ALFATClust that exploits rapid pairwise alignment-free sequence distance calculations and community detection in graph for clusters generation. Instead of a single threshold applied to every generated cluster, ALFATClust is capable of dynamically determining the cut-off threshold for each individual cluster by considering both cluster separation and intra-cluster sequence similarity. Benchmarking analysis shows that ALFATClust generally outperforms existing approaches by simultaneously maintaining cluster robustness and substantial cluster separation for the benchmark datasets. The software also provides an evaluation report for verifying the quality of the non-singleton clusters obtained.

CONCLUSIONS

ALFATClust is able to generate sequence clusters having high intra-cluster sequence similarity and substantial separation between clusters without having users to decide precise similarity cut-off thresholds.

Plasmid-Mediated Stabilization of Prophages

Mobile genetic elements (MGEs) drive bacterial evolution, alter gene availability within microbial communities, and facilitate adaptation to ecological niches. In natural systems, bacteria simultaneously possess or encounter multiple MGEs, yet their combined influences on microbial communities are poorly understood. Here, we investigate interactions among MGEs in the marine bacterium Sulfitobacter pontiacus. Two related strains, CB-D and CB-A, each harbor a single prophage. These prophages share high sequence identity with one another and an integration site within the host genome, yet these strains exhibit differences in “spontaneous” prophage induction (SPI) and consequent fitness. To better understand mechanisms underlying variation in SPI between these lysogens, we closed their genomes, which revealed that in addition to harboring different prophage genotypes, CB-A lacks two of the four large, low-copy-number plasmids possessed by CB-D. To assess the relative roles of plasmid content versus prophage genotype on host physiology, a panel of derivative strains varying in MGE content were generated. Characterization of these derivatives revealed a robust link between plasmid content and SPI, regardless of prophage genotype. Strains possessing all four plasmids had undetectable phage in cell-free lysates, while strains lacking either one plasmid (pSpoCB-1) or a combination of two plasmids (pSpoCB-2 and pSpoCB-4) produced high (>10⁵ PFU/mL) phage titers. Homologous plasmid sequences were identified in related bacteria, and plasmid and phage genes were found to be widespread in Tara Oceans metagenomic data sets. This suggests that plasmid-dependent stabilization of prophages may be commonplace throughout the oceans.

IMPORTANCE The consequences of prophage induction on the physiology of microbial populations are varied and include enhanced biofilm formation, conferral of virulence, and increased opportunity for horizontal gene transfer. These traits lead to competitive advantages for lysogenized bacteria and influence bacterial lifestyles in a variety of niches. However, biological controls of “spontaneous” prophage induction, the initiation of phage replication and phage-mediated cell lysis without an overt stressor, are not well understood. In this study, we observed a novel interaction between plasmids and prophages in the marine bacterium Sulfitobacter pontiacus. We found that loss of one or more distinct plasmids—which we show carry genes ubiquitous in the world’s oceans—resulted in a marked increase in prophage induction within lysogenized strains. These results demonstrate cross talk between different mobile genetic elements and have implications for our understanding of the lysogenic-lytic switches of prophages found not only in marine environments, but throughout all ecosystems.

Strain-level characterization of broad host range mobile genetic elements transferring antibiotic resistance from the human microbiome

Mobile genetic elements (MGEs) carrying antibiotic resistance genes (ARGs) disseminate ARGs when they mobilise into new bacterial hosts. The nature of such horizontal gene transfer (HGT) events between human gut commensals and pathogens remain poorly characterised. Here, we compare 1354 cultured commensal strains (540 species) to 45,403 pathogen strains (12 species) and find 64,188 MGE-mediated ARG transfer events between the two groups using established methods. Among the 5931 MGEs, we find 15 broad host range elements predicted to have crossed different bacterial phyla while also occurring in animal and environmental microbiomes. We experimentally demonstrate that predicted broad host range MGEs can mobilise from commensals Dorea longicatena and Hungatella hathewayi to pathogen Klebsiella oxytoca, crossing phyla simultaneously. Our work establishes the MGE-mediated ARG dissemination network between human gut commensals and pathogens and highlights broad host range MGEs as targets for future ARG dissemination management.

Healthcare-associated infections caused by chlorhexidine-tolerant Serratia marcescens carrying a promiscuous IncHI2 multi-drug resistance plasmid in a veterinary hospital

The bacterium Serratia marcescens can cause opportunistic infections in humans and in animals. In veterinary settings, the diversity, reservoirs and modes of transmission of this pathogen are poorly understood. The phenotypes and genotypes of Serratia spp. isolated from dogs, cats, horses, a bird and a rabbit examined at an Australian veterinary hospital between 2008 and 2019 were characterised. The isolates were identified as S. marcescens (n = 15) or S. ureilytica (n = 3) and were placed into four distinct phylogenetic groups. Nine quasi-clonal isolates associated with post-surgical complications in different patients displayed high levels of resistance to the antimicrobials fluoroquinolones, cephalosporins, aminoglycosides, and to the disinfectant chlorhexidine. A Serratia sp. with a similar resistance profile was also isolated from chlorhexidine solutions used across the Hospital, suggesting that these infections had a nosocomial origin. A genomic island encoding a homolog of the Pseudomonas MexCD-OprJ biocide efflux system was detected in the chlorhexidine-tolerant Serratia. The nine multi-drug resistant Serratia isolates also possessed a Ser-83-Ile mutation in GyrA conferring fluoroquinolone resistance, and carried a large IncHI2 conjugative plasmid encoding antimicrobial and heavy metal resistances. This replicon was highly similar to a plasmid previously detected in a strain of Enterobacter hormaechei recovered from the Hospital environment. IncHI2 plasmids are commonly found in Enterobacteriaceae, but are rarely present in Serratia spp., suggesting that this plasmid was acquired from another organism. A chlorhexidine-tolerant Serratia isolate which lacked the IncHI2 plasmid was used in mating experiments to demonstrate the transfer of multi-drug resistance from a E. hormaechei donor. This study illustrates the importance of environmental surveillance of biocide-resistance in veterinary hospitals.

The Dissemination of NDM-1 in Acinetobacter baumannii Strains

Acinetobacter baumannii is a common pathogen in hospitals and usually causes bacteremia, pneumonia, meningitis, peritonitis and other diseases. Isolates carried NDM-1 gene can make several antibiotics such as carbapenems and other beta-lactams ineffective. Nowadays, the number of A. baumannii strains carrying NDM-1 has been climbing year by year in recent years. To characterise the transmission of NDM-1 in A. baumannii, we collected 2576 human-derived genomes of A. baumannii strains from NCBI database and found that 186 strains contained NDM-1 gene. The multi-locus sequence typing, phylogenetic tree, NDM-1 gene organization and the single nucleotide polymorphisms of NDM-1 were investigated. We hope that our work will provide a theoretical basis for the prevention of dissemination of NDM-1 in A. baumannii.

Phylogenetic Comparison and Characterization of an mcr-1-Harboring Complete Plasmid Genome Isolated from Enterobacteriaceae

Global dissemination of mobilized colistin resistance (mcr)-1-carrying plasmids has been reported. This study aimed to investigate the global dissemination of these plasmids using whole genome sequencing to provide better understanding on genetic characteristics. Sixty-seven complete plasmid genomes harboring mcr-1 were obtained. Phylogeny was built against full plasmid genomes. Different replicon types of plasmid were compared in terms of antimicrobial resistance genes (ARGs), insertion sequence, and other functional genes. Five different replicon types of plasmid (IncX4, IncI2, IncP1, IncHIA, and IncFIB) were found to harbor mcr-1. IncX4 and IncI2 types of plasmid were well clustered in accordance with the country where they were isolated (and not as IncHIA and IncFIB). Three insertion sequences (ISApl1, ISKpn26, and IS1294) were identified in up- and/or downstream of mcr-1. Plasmids IncX4 and IncI2 were observed across the sample origin. Plasmids IncX4 showed high uniformity regardless of the origin of isolates and harbored H–NS coding genes, a facilitator for successful plasmid transfer. All three insertion sequences were observed in IncI2 plasmids. IncHI2 plasmids harbored various ARGs in addition to mcr-1. Our results elucidate the characteristics and phylogenetic relationships of complete mcr-1-harboring plasmids, indicating that global dissemination of mcr-1 is primarily owing to plasmid transfer rather than clonal spread.

Genome binning of viral entities from bulk metagenomics data

Despite the accelerating number of uncultivated virus sequences discovered in metagenomics and their apparent importance for health and disease, the human gut virome and its interactions with bacteria in the gastrointestinal tract are not well understood. This is partly due to a paucity of whole-virome datasets and limitations in current approaches for identifying viral sequences in metagenomics data. Here, combining a deep-learning based metagenomics binning algorithm with paired metagenome and metavirome datasets, we develop Phages from Metagenomics Binning (PHAMB), an approach that allows the binning of thousands of viral genomes directly from bulk metagenomics data, while simultaneously enabling clustering of viral genomes into accurate taxonomic viral populations. When applied on the Human Microbiome Project 2 (HMP2) dataset, PHAMB recovered 6,077 high-quality genomes from 1,024 viral populations, and identified viral-microbial host interactions. PHAMB can be advantageously applied to existing and future metagenomes to illuminate viral ecological dynamics with other microbiome constituents.

Ultra-deep long-read sequencing detects IS-mediated gene duplications as a potential trigger to generate arrays of resistance genes and a mechanism to induce novel gene variants such as blaCTX-M-243

BACKGROUND

Extended-spectrum β-lactamases (ESBLs) are enzymes that can render their hosts resistant to various β-lactam antibiotics. CTX-M-type enzymes are the most prevalent ESBLs and the main cause of resistance to third-generation cephalosporins in Enterobacteriaceae. The number of described CTX-M types is continuously rising, currently comprising over 240 variants. During routine screening we identified a novel blaCTX-M gene.

OBJECTIVES

To characterize a novel blaCTX-M variant harboured by a multidrug-resistant Escherichia coli isolate of sequence type ST354.

METHODS

Antibiotic susceptibilities were determined using broth microdilution. Genome and plasmid sequences were reconstructed using short- and long-read sequencing. The novel blaCTX-M locus was analysed using long-read and Sanger sequencing. Plasmid polymorphisms were determined in silico on a single plasmid molecule level.

RESULTS

The novel blaCTX-M-243 allele was discovered alongside a nearly identical blaCTX-M-104-containing gene array on a 219 kbp IncHI2A plasmid. CTX-M-243 differed from CTX-M-104 by only one amino acid substitution (N109S). Ultra-deep (2300-fold coverage) long-read sequencing revealed dynamic scaling of the blaCTX-M genetic contexts from one to five copies. Further antibiotic resistance genes such as blaTEM-1 also exhibited sequence heterogeneity but were stable in copy number.

CONCLUSIONS

We identified the novel ESBL gene blaCTX-M-243 and illustrate a dynamic system of varying blaCTX-M copy numbers. Our results highlight the constant emergence of new CTX-M family enzymes and demonstrate a potential evolutionary platform to generate novel ESBL variants and possibly other antibiotic resistance genes.

Comparative genomics of the black rot pathogen Xanthomonas campestris pv. campestris and non-pathogenic co-inhabitant Xanthomonas melonis from Trinidad reveal unique pathogenicity determinants and secretion system profiles

Black-rot disease caused by the phytopathogen Xanthomonas campestris pv. campestris (Xcc) continues to have considerable impacts on the productivity of cruciferous crops in Trinidad and Tobago and the wider Caribbean region. While the widespread occurrence of resistance of Xcc against bactericidal agrochemicals can contribute to the high disease burdens, the role of virulence and pathogenicity features of local strains on disease prevalence and severity has not been investigated yet. In the present study, a comparative genomic analysis was performed on 6 pathogenic Xcc and 4 co-isolated non-pathogenic Xanthomonas melonis (Xmel) strains from diseased crucifer plants grown in fields with heavy chemical use in Trinidad. Native isolates were grouped into two known and four newly assigned ribosomal sequence types (rST). Mobile genetic elements were identified which belonged to the IS3, IS5 family, Tn3 transposon, resolvases, and tra T4SS gene clusters. Additionally, exogenous plasmid derived sequences with origins from other bacterial species were characterised. Although several instances of genomic rearrangements were observed, native Xcc and Xmel isolates shared a significant level of structural homology with reference genomes, Xcc ATCC 33913 and Xmel CFBP4644, respectively. Complete T1SS hlyDB, T2SS, T4SS vir and T5SS xadA, yapH and estA gene clusters were identified in both species. Only Xmel strains contained a complete T6SS but no T3SS. Both species contained a complex repertoire of extracellular cell wall degrading enzymes. Native Xcc strains contained 37 T3SS and effector genes but a variable and unique profile of 8 avr, 4 xop and 1 hpa genes. Interestingly, Xmel strains contained several T3SS effectors with low similarity to references including avrXccA1 (~89%), hrpG (~73%), hrpX (~90%) and xopAZ (~87%). Furthermore, only Xmel genomes contained a CRISPR-Cas I-F array, but no lipopolysaccharide wxc gene cluster. Xmel strains were confirmed to be non-pathogenic by pathogenicity assays. The results of this study will be useful to guide future research into virulence mechanisms, agrochemical resistance, pathogenomics and the potential role of the co-isolated non-pathogenic Xanthomonas strains on Xcc infections.

IncHI1A plasmids potentially facilitate a horizontal flow of antibiotic resistance genes to pathogens in microbial communities of urban residential sewage

Horizontal gene transfer via plasmids is important for the dissemination of antibiotic resistance genes among medically relevant pathogens. Specifically, the transfer of IncHI1A plasmids is believed to facilitate the spread of antibiotic resistance genes, such as carbapenemases, within the clinically important family Enterobacteriaceae. The microbial community of urban wastewater treatment plants has been shown to be highly permissive towards conjugal transfer of IncP1 plasmids. Here, we tracked the transfer of the P1 plasmid pB10 and the clinically relevant HI1A plasmid R27 in the microbial communities present in urban residential sewage entering full-scale wastewater treatment plants. We found that both plasmids readily transferred to these communities and that strains in the sewage were able to further disseminate them. Furthermore, that R27 has a broad potential host range, but a low host divergence. Interestingly, although the majority of R27 transfer events were to members of Enterobacteriaceae, we found a subset of transfer to other families, even other phyla. Indicating, that HI1A plasmids facilitate horizontal gene transfer both within Enterobacteriaceae, but also across families of especially Gammaproteobacteria, such as Moraxellaceae, Pseudomonadaceae and Shewanellaceae. pB10 displayed a similar potential host range as R27. In contrast to R27, pB10 had a high host divergence. By culture enrichment of the transconjugant communities, we show that sewage strains of Enterobacteriaceae and Aeromonadaceae can stably maintain R27 and pB10, respectively. Our results suggest that dissemination in the urban residual water system of HI1A plasmids may result in an accelerated acquisition of antibiotic resistance genes among pathogens.

OUP accepted manuscript

Motivation

Phage–host associations play important roles in microbial communities. But in natural communities, as opposed to culture-based lab studies where phages are discovered and characterized metagenomically, their hosts are generally not known. Several programs have been developed for predicting which phage infects which host based on various sequence similarity measures or machine learning approaches. These are often based on whole viral and host genomes, but in metagenomics-based studies, we rarely have whole genomes but rather must rely on contigs that are sometimes as short as hundreds of bp long. Therefore, we need programs that predict hosts of phage contigs on the basis of these short contigs. Although most existing programs can be applied to metagenomic datasets for these predictions, their accuracies are generally low. Here, we develop ContigNet, a convolutional neural network-based model capable of predicting phage–host matches based on relatively short contigs, and compare it to previously published VirHostMatcher (VHM) and WIsH.

Results

On the validation set, ContigNet achieves 72–85% area under the receiver operating characteristic curve (AUROC) scores, compared to the maximum of 68% by VHM or WIsH for contigs of lengths between 200 bps to 50 kbps. We also apply the model to the Metagenomic Gut Virus (MGV) catalogue, a dataset containing a wide range of draft genomes from metagenomic samples and achieve 60–70% AUROC scores compared to that of VHM and WIsH of 52%. Surprisingly, ContigNet can also be used to predict plasmid-host contig associations with high accuracy, indicating a similar genetic exchange between mobile genetic elements and their hosts.

Availability and implementation

The source code of ContigNet and related datasets can be downloaded from https://github.com/tianqitang1/ContigNet.

Complete sequence of classic F-type plasmid pRK100 shows unique conservation over time and geographic location

Plasmids exhibit great diversity of gene content and host ranges and are famous for quick adaptation to the genetic background of the bacterial host cell. In addition to observing ever evolving plasmids, some plasmids have conserved backbones: a stable core composition and arrangement of genes in addition to variable regions. There are a few reports of extremely conserved plasmids. Here we report the complete sequence of pRK100 plasmid - a large, well-characterized conjugative F-like plasmid found in an Escherichia coli strain isolated from a urinary tract infection patient in 1990. The sequence shows that the 142 kb-long pRK100 plasmid is nearly identical to plasmids circulating in distant geographical locations and found in different host E. coli strains between 2007 and 2017. We also performed additional functional characterization of pRK100. Our results showed that pRK100 does not have a strong pathogenicity phenotype in porcine primary bladder epithelial cell culture. Moreover, the conjugation of pRK100 seems to strongly depend on recipient characteristics. These observations and identification of the pRK100 plasmid in different strain genotypes leave the extreme sequence conservation and broad distribution of this plasmid unexplained.

Metagenomic evidence for co-occurrence of antibiotic, biocide and metal resistance genes in pigs

Antibiotic-resistant pathogens constitute an escalating public health concern. Hence a better understanding of the underlying processes responsible for this expansion is urgently needed. Co-selection of heavy metal/biocide and antibiotic resistance genes (ARGs) has been suggested as one potential mechanism promoting the proliferation of antimicrobial resistance (AMR). This paper aims to elucidate this interplay and exploit differences in antibiotic usage to infer patterns of co-selection by the non-antibiotic factors metals and biocides in the context of pig farming. We examined 278 gut metagenomes from pigs with continuous antibiotic exposure, only at weaning and at no exposure. Metals as growth promoters and biocides as disinfectants are currently used with little restrictions in stock farming. The pigs under continuous antibiotic exposure displayed the highest co-occurrence of ARGs and other genetic elements while the pigs under limited use of antibiotics still showed abundant co-occurrences. Pathogens belonging to Enterobacteriaceae displayed increased co-occurrence phenomena, suggesting that this maintenance is not a random selection process from a mobilized pool but pertains to specific phylogenetic clades. These results suggest that metals and biocides displayed strong selective pressures on ARGs exerted by intensive farming, regardless of the current use of antibiotics.

Plant Growth Promotion Diversity in Switchgrass-Colonizing, Diazotrophic Endophytes

Endophytic nitrogen-fixing (diazotrophic) bacteria are essential members of the microbiome of switchgrass (Panicum virgatum), considered to be an important commodity crop in bioenergy production. While endophytic diazotrophs are known to provide fixed atmospheric nitrogen to their host plant, there are many other plant growth-promoting (PGP) capabilities of these organisms to be demonstrated. The diversity of PGP traits across different taxa of switchgrass-colonizing endophytes is understudied, yet critical for understanding endophytic function and improving cultivation methods of important commodity crops. Here, we present the isolation and characterization of three diazotrophic endophytes: Azospirillum agricola R1C, Klebsiella variicola F10Cl, and Raoultella terrigena R1Gly. Strains R1C and F10Cl were isolated from switchgrass and strain R1Gly, while isolated from tobacco, is demonstrated herein to colonize switchgrass. Each strain exhibited highly diverse genomic and phenotypic PGP capabilities. Strain F10Cl and R1Gly demonstrated the highest functional similarity, suggesting that, while endophyte community structure may vary widely based on host species, differences in functional diversity are not a clearly delineated. The results of this study advance our understanding of diazotrophic endophyte diversity, which will allow us to design robust strategies to improve cultivation methods of many economically important commodity crops.

Deeplasmid: deep learning accurately separates plasmids from bacterial chromosomes

Plasmids are mobile genetic elements that play a key role in microbial ecology and evolution by mediating horizontal transfer of important genes, such as antimicrobial resistance genes. Many microbial genomes have been sequenced by short read sequencers and have resulted in a mix of contigs that derive from plasmids or chromosomes. New tools that accurately identify plasmids are needed to elucidate new plasmid-borne genes of high biological importance. We have developed Deeplasmid, a deep learning tool for distinguishing plasmids from bacterial chromosomes based on the DNA sequence and its encoded biological data. It requires as input only assembled sequences generated by any sequencing platform and assembly algorithm and its runtime scales linearly with the number of assembled sequences. Deeplasmid achieves an AUC–ROC of over 89%, and it was more accurate than five other plasmid classification methods. Finally, as a proof of concept, we used Deeplasmid to predict new plasmids in the fish pathogen Yersinia ruckeri ATCC 29473 that has no annotated plasmids. Deeplasmid predicted with high reliability that a long assembled contig is part of a plasmid. Using long read sequencing we indeed validated the existence of a 102 kb long plasmid, demonstrating Deeplasmid's ability to detect novel plasmids.