PlaScope: a targeted approach to assess the plasmidome from genome assemblies at the species level

IncC plasmid genome rearrangements influence the vertical and horizontal transmission tradeoff in Escherichia coli

It has been shown that an evolutionary tradeoff between vertical (host growth rate) and horizontal (plasmid conjugation) transmissions contributes to global plasmid fitness. As conjugative IncC plasmids are important for the spread of multidrug resistance (MDR), in a broad range of bacterial hosts, we investigated vertical and horizontal transmissions of two multidrug-resistant IncC plasmids according to their backbones and MDR-region rearrangements, upon plasmid entry into a new host. We observed plasmid genome deletions after conjugation in three diverse natural Escherichia coli clinical strains, varying from null to high number depending on the plasmid, all occurring in the MDR region. The plasmid burden on bacterial fitness depended more on the strain background than on the structure of the MDR region, with deletions appearing to have no impact. Besides, we observed an increase in plasmid transfer rate, from ancestral host to new clinical recipient strains, when the IncC plasmid was rearranged. Finally, using a second set of conjugation experiments, we investigated the evolutionary tradeoff of the IncC plasmid during the critical period of plasmid establishment in E. coli K-12, by correlating the transfer rates of deleted or non-deleted IncC plasmids and their costs on the recipient strain. Plasmid deletions strongly improved conjugation efficiency with no negative growth effect. Our findings indicate that the flexibility of the MDR-region of the IncC plasmids can promote their dissemination, and provide diverse opportunities to capture new resistance genes. In a broader view, they suggest that the vertical-horizontal transmission tradeoff can be manipulated by the plasmid to improve its fitness.

Genomic evidence of Escherichia coli gut population diversity translocation in leukemia patients

Escherichia coli, a commensal species of the human gut, is an opportunistic pathogen that can reach extra-intestinal compartments, including the bloodstream and the bladder, among others. In non-immunosuppressed patients, purifying or neutral evolution of E. coli populations has been reported in the gut. Conversely, it has been suggested that when migrating to extra-intestinal compartments, E. coli genomes undergo diversifying selection as supported by strong evidence for adaptation. The level of genomic polymorphism and the size of the populations translocating from gut to extra-intestinal compartments is largely unknown. To gain insights into the pathophysiology of these translocations, we investigated the level of polymorphism and the evolutionary forces acting on the genomes of 77 E. coli isolated from various compartments in three immunosuppressed patients. Each patient had a unique strain, which was a mutator in one case. In all instances, we observed that translocation encompasses much of the genomic diversity present in the gut. The same signature of selection, whether purifying or diversifying, and as anticipated, neutral for mutator isolates, was observed in both the gut and bloodstream. Additionally, we found a limited number of non-specific mutations among compartments for non-mutator isolates. In all cases, urine isolates were dominated by neutral selection. These findings indicate that substantial proportions of populations are undergoing translocation and that they present a complex compartment-specific pattern of selection at the patient level.IMPORTANCEIt has been suggested that intra and extra-intestinal compartments differentially constrain the evolution of E. coli strains. Whether host particular conditions, such as immunosuppression, could affect the strain evolutionary trajectories remains understudied. We found that, in immunosuppressed patients, large fractions of E. coli gut populations are translocating with variable modifications of the signature of selection for commensal and pathogenic isolates according to the compartment and/or the patient. Such multiple site sampling should be performed in large cohorts of patients to gain a better understanding of E. coli extra-intestinal diseases.

PlasmidEC and gplas2: an optimized short-read approach to predict and reconstruct antibiotic resistance plasmids in Escherichia coli

Accurate reconstruction of Escherichia coli antibiotic resistance gene (ARG) plasmids from Illumina sequencing data has proven to be a challenge with current bioinformatic tools. In this work, we present an improved method to reconstruct E. coli plasmids using short reads. We developed plasmidEC, an ensemble classifier that identifies plasmid-derived contigs by combining the output of three different binary classification tools. We showed that plasmidEC is especially suited to classify contigs derived from ARG plasmids with a high recall of 0.941. Additionally, we optimized gplas, a graph-based tool that bins plasmid-predicted contigs into distinct plasmid predictions. Gplas2 is more effective at recovering plasmids with large sequencing coverage variations and can be combined with the output of any binary classifier. The combination of plasmidEC with gplas2 showed a high completeness (median=0.818) and F1-Score (median=0.812) when reconstructing ARG plasmids and exceeded the binning capacity of the reference-based method MOB-suite. In the absence of long-read data, our method offers an excellent alternative to reconstruct ARG plasmids in E. coli.

Dissemination of IncI plasmid encoding bla CTX-M-1 is not hampered by its fitness cost in the pig's gut

Multiresistance plasmids belonging to the IncI incompatibility group have become one of the most pervasive plasmid types in extended-spectrum beta-lactamase-producing Escherichia coli of animal origin. The extent of the burden imposed on the bacterial cell by these plasmids seems to modulate the emergence of "epidemic" plasmids. However, in vivo data in the natural environment of the strains are scarce. Here, we investigated the cost of a bla CTX-M-1-IncI1 epidemic plasmid in a commensal E. coli animal strain, UB12-RC, before and after oral inoculation of 15 6- to 8-week- old specific-pathogen-free pigs. Growth rate in rich medium was determined on (i) UB12-RC and derivatives, with or without plasmid, in vivo and/or in vitro evolved, and (ii) strains that acquired the plasmid in the gut during the experiment. Although bla CTX-M-1-IncI1 plasmid imposed no measurable burden on the recipient strain after conjugation and during the longitudinal carriage in the pig's gut, we observed a significant difference in the bacterial growth rate between IncI1 plasmid-carrying and plasmid-free isolates collected during in vivo carriage. Only a few mutations on the chromosome of the UB12-RC derivatives were detected by whole-genome sequencing. RNA-Seq analysis of a selected set of these strains showed that transcriptional responses to the bla CTX-M-1-IncI1 acquisition were limited, affecting metabolism, stress response, and motility functions. Our data suggest that the effect of IncI plasmid on host cells is limited, fitness cost being insufficient to act as a barrier to IncI plasmid spread among natural population of E. coli in the gut niche.

Faecal carriage of extended-spectrum β-lactamase-producing Escherichia coli in a remote region of Niger

OBJECTIVE

Whole genome sequencing (WGS) of extended-spectrum β-lactamase-producing Escherichia coli (ESBL-E. coli) in developing countries is lacking. Here we describe the population structure and molecular characteristics of ESBL-E. coli faecal isolates in rural Southern Niger.

METHODS

Stools of 383 healthy participants were collected among which 92.4% were ESBL-Enterobacterales carriers. A subset of 90 ESBL-E. coli containing stools (109 ESBL-E. coli isolates) were further analysed by WGS, using short- and long-reads.

RESULTS

Most isolates belonged to the commensalism-adapted phylogroup A (83.5%), with high clonal diversity. The blaCTX-M-15 gene was the major ESBL determinant (98.1%), chromosome-integrated in approximately 50% of cases, in multiple integration sites. When plasmid-borne, blaCTX-M-15 was found in IncF (57.4%) and IncY plasmids (26.2%). Closely related plasmids were found in different genetic backgrounds. Genomic environment analysis of blaCTX-M-15 in closely related strains argued for mobilisation between plasmids or from plasmid to chromosome.

CONCLUSIONS

Massive prevalence of community faecal carriage of CTX-M-15-producing E. coli was observed in a rural region of Niger due to the spread of highly diverse A phylogroup commensalism-adapted clones, with frequent chromosomal integration of blaCTX-M-15. Plasmid spread was also observed. These data suggest a risk of sustainable implementation of ESBL in community faecal carriage.

Characterization of the Environmental Plasmidome of the Red Sea

Plasmids contribute to microbial diversity and adaptation, providing microorganisms with the ability to thrive in a wide range of conditions in extreme environments. However, while the number of marine microbiome studies is constantly increasing, very little is known about marine plasmids, and they are very poorly represented in public databases. To extend the repertoire of environmental marine plasmids, we established a pipeline for the de novo assembly of plasmids in the marine environment by analyzing available microbiome metagenomic sequencing data. By applying the pipeline to data from the Red Sea, we identified 362 plasmid candidates. We showed that the distribution of plasmids corresponds to environmental conditions, particularly, depth, temperature, and physical location. At least 7 of the 362 candidates are most probably real plasmids, based on a functional analysis of their open reading frames (ORFs). Only one of the seven has been described previously. Three plasmids were identified in other public marine metagenomic data from different locations all over the world; these plasmids contained different cassettes of functional genes at each location. Analysis of antibiotic and metal resistance genes revealed that the same positions that were enriched with genes encoding resistance to antibiotics were also enriched with resistance to metals, suggesting that plasmids contribute site-dependent phenotypic modules to their ecological niches. Finally, half of the ORFs (50.8%) could not be assigned to a function, emphasizing the untapped potential of the unique marine plasmids to provide proteins with multiple novel functions. IMPORTANCE Marine plasmids are understudied and hence underrepresented in databases. Plasmid functional annotation and characterization is complicated but, if successful, may provide a pool of novel genes and unknown functions. Newly discovered plasmids and their functional repertoire are potentially valuable tools for predicting the dissemination of antimicrobial resistance, providing vectors for molecular cloning and an understanding of plasmid-bacterial interactions in various environments.

Epistatic interactions between the high pathogenicity island and other iron uptake systems shape Escherichia coli extra-intestinal virulence

The intrinsic virulence of extra-intestinal pathogenic Escherichia coli is associated with numerous chromosomal and/or plasmid-borne genes, encoding diverse functions such as adhesins, toxins, and iron capture systems. However, the respective contribution to virulence of those genes seems to depend on the genetic background and is poorly understood. Here, we analyze genomes of 232 strains of sequence type complex STc58 and show that virulence (quantified in a mouse model of sepsis) emerged in a sub-group of STc58 due to the presence of the siderophore-encoding high-pathogenicity island (HPI). When extending our genome-wide association study to 370 Escherichia strains, we show that full virulence is associated with the presence of the aer or sit operons, in addition to the HPI. The prevalence of these operons, their co-occurrence and their genomic location depend on strain phylogeny. Thus, selection of lineage-dependent specific associations of virulence-associated genes argues for strong epistatic interactions shaping the emergence of virulence in E. coli.

Plasmer: an Accurate and Sensitive Bacterial Plasmid Prediction Tool Based on Machine Learning of Shared k-mers and Genomic Features

Identification of plasmids in bacterial genomes is critical for many factors, including horizontal gene transfer, antibiotic resistance genes, host-microbe interactions, cloning vectors, and industrial production. There are several in silico methods to predict plasmid sequences in assembled genomes. However, existing methods have evident shortcomings, such as unbalance in sensitivity and specificity, dependency on species-specific models, and performance reduction in sequences shorter than 10 kb, which has limited their scope of applicability. In this work, we proposed Plasmer, a novel plasmid predictor based on machine-learning of shared k-mers and genomic features. Unlike existing k-mer or genomic-feature based methods, Plasmer employs the random forest algorithm to make predictions using the percent of shared k-mers with plasmid and chromosome databases combined with other genomic features, including alignment E value and replicon distribution scores (RDS). Plasmer can predict on multiple species and has achieved an average the area under the curve (AUC) of 0.996 with accuracy of 98.4%. Compared to existing methods, tests of both sliding sequences and simulated and de novo assemblies have consistently shown that Plasmer has outperforming accuracy and stable performance across long and short contigs above 500 bp, demonstrating its applicability for fragmented assemblies. Plasmer also has excellent and balanced performance on both sensitivity and specificity (both >0.95 above 500 bp) with the highest F1-score, which has eliminated the bias on sensitivity or specificity that was common in existing methods. Plasmer also provides taxonomy classification to help identify the origin of plasmids. IMPORTANCE In this study, we proposed a novel plasmid prediction tool named Plasmer. Technically, unlike existing k-mer or genomic features-based methods, Plasmer is the first tool to combine the advantages of the percent of shared k-mers and the alignment score of genomic features. This has given Plasmer (i) evident improvement in performance compared to other methods, with the best F1-score and accuracy on sliding sequences, simulated contigs, and de novo assemblies; (ii) applicability for contigs above 500 bp with highest accuracy, enabling plasmid prediction in fragmented short-read assemblies; (iii) excellent and balanced performance between sensitivity and specificity (both >0.95 above 500 bp) with the highest F1-score, which eliminated the bias on sensitivity or specificity that commonly existed in other methods; and (iv) no dependency of species-specific training models. We believe that Plasmer provides a more reliable alternative for plasmid prediction in bacterial genome assemblies.

First report of enterotoxigenic Staphylococcus argenteus as a foodborne pathogen

Staphylococcal enterotoxins preformed in food are the causative agents of staphylococcal food poisoning outbreaks (SFPO). In this study we characterised in depth two coagulase-positive non-pigmented staphylococci involved in two independent outbreaks that occurred in France. While indistinguishable from Staphylococcus aureus using PCR methods and growth phenotype comparisons, both isolates were identified as Staphylococcus argenteus by whole genome sequencing. The genomes were analysed for the presence of enterotoxin genes, whose expression was determined in laboratory medium and, for the first time, in artificially-contaminated milk samples by using liquid chromatography-mass spectrometry and ELISA methods. The concentration measured for the SEB toxin in milk (0.67 ng/ml) was comparable to concentrations reported for other types of enterotoxins behind SFPO. From a collection of publicly available genomes, we performed an unprecedented systematic investigation of the enterotoxin gene set of S. argenteus, including variants and pseudogenes. The most prevalent genes were sex, followed by sel26, sel27 and sey. The egc cluster was less frequent and most of the time carried a dysfunctional seg gene. Our results shed light on the enterotoxigenic properties of S. argenteus, and emphasize the importance in monitoring of S. argenteus as an emerging foodborne pathogen.

Evolution of fluoroquinolone-resistant Escherichia coli in the gut after ciprofloxacin treatment

BACKGROUND

Three healthy volunteers carried similar quinolone-resistant E. coli (QREC) (pulsed field gel electrophoresis profiles) in their gut before and after 14 days ciprofloxacin treatment. Given the intensity of the selective pressure and the mutagenic properties of quinolones, we determined whether these strains had evolved at the phenotypic and/or genomic levels.

MATERIAL AND METHODS

Commensal QREC from before day-0 (D0), and a month after 14 days of ciprofloxacin (D42) were compared in 3 volunteers. Growth experiments were performed; acetate levels, mutation frequencies, quinolone MICs and antibiotic tolerance were measured at D0 and D42. Genomes were sequenced and single nucleotide polymorphisms (SNPs) between D0 and D42 were analyzed using DiscoSNP and breseq methods. Cytoplasmic proteins were extracted, HPLC performed and proteins identified using X!tandem software; abundances were measured by mass spectrometry using the Spectral Counting (SC) and eXtraction Ion Chromatograms (XIC) integration methods.

RESULTS

No difference was found in MICs, growth characteristics, acetate concentrations, mutation frequencies, tolerance profiles, phylogroups, O-and H-types, fimH alleles and sequence types between D0 and D42. No SNP variation was evidenced between D0 and D42 isolates for 2/3 subjects; 2 SNP variations were evidenced in one. At the protein level, very few significant protein abundance differences were identified between D0 and D42.

CONCLUSION

No fitness, tolerance, metabolic or genomic evolution of commensal QREC was observed overtime, despite massive exposure to ciprofloxacin in the gut. The three strains behaved as if they had been unaffected by ciprofloxacin, suggesting that gut may act as a sanctuary where bacteria would be protected from the effect of antibiotics and survive without any detrimental effect of stress.

Prevalence, Risk Factors, and Genetic Characterization of Extended-Spectrum Beta-Lactamase Escherichia coli Isolated From Healthy Pregnant Women in Madagascar

Antimicrobial resistance is a major public health concern worldwide affecting humans, animals and the environment. However, data is lacking especially in developing countries. Thus, the World Health Organization developed a One-Health surveillance project called Tricycle focusing on the prevalence of ESBL-producing Escherichia coli in humans, animals, and the environment. Here we present the first results of the human community component of Tricycle in Madagascar. From July 2018 to April 2019, rectal swabs from 492 pregnant women from Antananarivo, Mahajanga, Ambatondrazaka, and Toamasina were tested for ESBL-E. coli carriage. Demographic, sociological and environmental risk factors were investigated, and E. coli isolates were characterized (antibiotic susceptibility, resistance and virulence genes, plasmids, and genomic diversity). ESBL-E. coli prevalence carriage in pregnant women was 34% varying from 12% (Toamasina) to 65% (Ambatondrazaka). The main risk factor associated with ESBL-E. coli carriage was the rainy season (OR = 2.9, 95% CI 1.3-5.6, p = 0.009). Whole genome sequencing was performed on 168 isolates from 144 participants. bla CTX-M-15 was the most frequent ESBL gene (86%). One isolate was resistant to carbapenems and carried the bla NDM-5 gene. Most isolates belonged to commensalism associated phylogenetic groups A, B1, and C (90%) and marginally to extra-intestinal virulence associated phylogenetic groups B2, D and F (10%). Multi locus sequence typing showed 67 different sequence types gathered in 17 clonal complexes (STc), the most frequent being STc10/phylogroup A (35%), followed distantly by the emerging STc155/phylogroup B1 (7%), STc38/phylogroup D (4%) and STc131/phylogroup B2 (3%). While a wide diversity of clones has been observed, SNP analysis revealed several genetically close isolates (n = 34/168) which suggests human-to-human transmissions. IncY plasmids were found with an unusual prevalence (23%), all carrying a bla CTX-M-15. Most of them (85%) showed substantial homology (≥85%) suggesting a dissemination of IncY ESBL plasmids in Madagascar. This large-scale study reveals a high prevalence of ESBL-E. coli among pregnant women in four cities in Madagascar associated with warmth and rainfall. It shows the great diversity of E. coli disseminating throughout the country but also transmission of specific clones and spread of plasmids. This highlights the urgent need of public-health interventions to control antibiotic resistance in the country.

RFPlasmid: predicting plasmid sequences from short-read assembly data using machine learning

Antimicrobial-resistance (AMR) genes in bacteria are often carried on plasmids and these plasmids can transfer AMR genes between bacteria. For molecular epidemiology purposes and risk assessment, it is important to know whether the genes are located on highly transferable plasmids or in the more stable chromosomes. However, draft whole-genome sequences are fragmented, making it difficult to discriminate plasmid and chromosomal contigs. Current methods that predict plasmid sequences from draft genome sequences rely on single features, like k-mer composition, circularity of the DNA molecule, copy number or sequence identity to plasmid replication genes, all of which have their drawbacks, especially when faced with large single-copy plasmids, which often carry resistance genes. With our newly developed prediction tool RFPlasmid, we use a combination of multiple features, including k-mer composition and databases with plasmid and chromosomal marker proteins, to predict whether the likely source of a contig is plasmid or chromosomal. The tool RFPlasmid supports models for 17 different bacterial taxa, including Campylobacter, Escherichia coli and Salmonella, and has a taxon agnostic model for metagenomic assemblies or unsupported organisms. RFPlasmid is available both as a standalone tool and via a web interface.

Interplay between Bacterial Clones and Plasmids in the Spread of Antibiotic Resistance Genes in the Gut: Lessons from a Temporal Study in Veal Calves

Intestinal carriage of extended spectrum β-lactamase (ESBL)-producing Escherichia coli is a frequent, increasing, and worrying phenomenon, but little is known about the molecular scenario and the evolutionary forces at play. We screened 45 veal calves, known to have high prevalence of carriage, for ESBL-producing E. coli on 514 rectal swabs (one randomly selected colony per sample) collected over 6 months. We characterized the bacterial clones and plasmids carrying bla_ESBL genes with a combination of genotyping methods, whole genome sequencing, and conjugation assays. One hundred and seventy-three ESBL-producing E. coli isolates [bla_CTX-M-1 (64.7%), bla_CTX-M-14 (33.5%), or bla_CTX-M-15 (1.8%)] were detected, belonging to 32 bacterial clones, mostly of phylogroup A. Calves were colonized successively by different clones with a trend in decreasing carriage. The persistence of a clone in a farm was significantly associated with the number of calves colonized. Despite a high diversity of E. coli clones and bla_CTX-M-carrying plasmids, few bla_CTX-M gene/plasmid/chromosomal background combinations dominated, due to (i) efficient colonization of bacterial clones and/or (ii) successful plasmid spread in various bacterial clones. The scenario "clone versus plasmid spread" depended on the farm. Thus, epistatic interactions between resistance genes, plasmids, and bacterial clones contribute to optimize fitness in specific environments. IMPORTANCE The gut microbiota is the epicenter of the emergence of resistance. Considerable amount of knowledge on the molecular mechanisms of resistance has been accumulated, but the ecological and evolutionary forces at play in nature are less studied. In this context, we performed a field work on temporal intestinal carriage of extended spectrum β-lactamase (ESBL)-producing Escherichia coli in veal farms. Veal calves are animals with one of the highest levels of ESBL producing E. coli fecal carriage, due to early high antibiotic exposure. We were able to show that calves were colonized successively by different ESBL-producing E. coli clones, and that two main scenarios were at play in the spread of bla_CTX-M genes among calves: efficient colonization of several calves by a few bacterial clones and successful plasmid spread in various bacterial clones. Such knowledge should help develop new strategies to fight the emergence of antibiotic-resistance.

Recovering Escherichia coli Plasmids in the Absence of Long-Read Sequencing Data

The incidence of infections caused by multidrug-resistant E. coli strains has risen in the past years. Antibiotic resistance in E. coli is often mediated by acquisition and maintenance of plasmids. The study of E. coli plasmid epidemiology and genomics often requires long-read sequencing information, but recently a number of tools that allow plasmid prediction from short-read data have been developed. Here, we reviewed 25 available plasmid prediction tools and categorized them into binary plasmid/chromosome classification tools and plasmid reconstruction tools. We benchmarked six tools (MOB-suite, plasmidSPAdes, gplas, FishingForPlasmids, HyAsP and SCAPP) that aim to reliably reconstruct distinct plasmids, with a special focus on plasmids carrying antibiotic resistance genes (ARGs) such as extended-spectrum beta-lactamase genes. We found that two thirds (n = 425, 66.3%) of all plasmids were correctly reconstructed by at least one of the six tools, with a range of 92 (14.58%) to 317 (50.23%) correctly predicted plasmids. However, the majority of plasmids that carried antibiotic resistance genes (n = 85, 57.8%) could not be completely recovered as distinct plasmids by any of the tools. MOB-suite was the only tool that was able to correctly reconstruct the majority of plasmids (n = 317, 50.23%), and performed best at reconstructing large plasmids (n = 166, 46.37%) and ARG-plasmids (n = 41, 27.9%), but predictions frequently contained chromosome contamination (40%). In contrast, plasmidSPAdes reconstructed the highest fraction of plasmids smaller than 18 kbp (n = 168, 61.54%). Large ARG-plasmids, however, were frequently merged with sequences derived from distinct replicons. Available bioinformatic tools can provide valuable insight into E. coli plasmids, but also have important limitations. This work will serve as a guideline for selecting the most appropriate plasmid reconstruction tool for studies focusing on E. coli plasmids in the absence of long-read sequencing data.

Dynamics of extended-spectrum beta-lactamase-producing Enterobacterales colonization in long-term carriers following travel abroad

Travel to tropical regions is associated with high risk of acquiring extended-spectrum beta-lactamase-producing Enterobacterales (ESBL-E) that are typically cleared in less than 3 months following return. The conditions leading to persistent carriage that exceeds 3 months in some travellers require investigation. Whole-genome sequencing (Illumina MiSeq) was performed on the 82 ESBL-E isolates detected upon return and 1, 2, 3, 6 and 12 months later from the stools of 11 long-term (>3 months) ESBL-E carriers following travel abroad. One to five different ESBL Escherichia coli strains were detected per traveller upon return, and this diminished to one after 3 months. Long-term carriage was due to the presence of the same ESBL E. coli strain, for more than 3 months, in 9 out of 11 travellers, belonging to epidemic sequence type complexes (STc 10, 14, 38, 69, 131 and 648). The mean carriage duration of strains belonging to phylogroups B2/D/F, associated with extra-intestinal virulence, was higher than that for commensal-associated A/B1/E phylogroups (3.5 vs 0.5 months, P=0.021). Genes encoding iron capture systems (fyuA, irp), toxins (senB, sat), adhesins (flu, daaF, afa/nfaE, pap, ecpA) and colicin (cjrA) were more often present in persistent strains than in transient ones. Single-nucleotide polymorphism (SNP) analysis in persistent strains showed a maximum divergence of eight SNPs over 12 months without signs of adaptation. Genomic plasticity was observed during the follow-up with the loss or gain of mobile genetic elements such as plasmids, integrons and/or transposons that may contain resistance genes at different points in the follow-up. Long-term colonization of ESBL-E following travel is primarily due to the acquisition of E. coli strains belonging to epidemic clones and harbouring ‘virulence genes’, allowing good adaptation to the intestinal microbiota.

Phylogroup stability contrasts with high within sequence type complex dynamics of Escherichia coli bloodstream infection isolates over a 12-year period

Background

Escherichia coli is the leading cause of bloodstream infections, associated with a significant mortality. Recent genomic analyses revealed that few clonal lineages are involved in bloodstream infections and captured the emergence of some of them. However, data on within sequence type (ST) population genetic structure evolution are rare.

Methods

We compared whole genome sequences of 912 E. coli isolates responsible for bloodstream infections from two multicenter clinical trials that were conducted in the Paris area, France, 12 years apart, in teaching hospitals belonging to the same institution (“Assistance Publique-Hôpitaux de Paris”). We analyzed the strains at different levels of granularity, i.e., the phylogroup, the ST complex (STc), and the within STc clone taking into consideration the evolutionary history, the resistance, and virulence gene content as well as the antigenic diversity of the strains.

Results

We found a mix of stability and changes overtime, depending on the level of comparison. Overall, we observed an increase in antibiotic resistance associated to a restricted number of genetic determinants and in strain plasmidic content, whereas phylogroup distribution and virulence gene content remained constant. Focusing on STcs highlighted the pauci-clonality of the populations, with only 11 STcs responsible for more than 73% of the cases, dominated by five STcs (STc73, STc131, STc95, STc69, STc10). However, some STcs underwent dramatic variations, such as the global pandemic STc131, which replaced the previously predominant STc95. Moreover, within STc131, 95 and 69 genomic diversity analysis revealed a highly dynamic pattern, with reshuffling of the population linked to clonal replacement sometimes coupled with independent acquisitions of virulence factors such as the pap gene cluster bearing a papGII allele located on various pathogenicity islands. Additionally, STc10 exhibited huge antigenic diversity evidenced by numerous O:H serotype/fimH allele combinations, whichever the year of isolation.

Conclusions

Altogether, these data suggest that the bloodstream niche is occupied by a wide but specific phylogenetic diversity and that highly specialized extra-intestinal clones undergo frequent turnover at the within ST level. Additional worldwide epidemiological studies overtime are needed in different geographical and ecological contexts to assess how generalizable these data are.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-00892-0.

Detection of plasmid contigs in draft genome assemblies using customized Kraken databases

Plasmids play an important role in bacterial evolution and mediate horizontal transfer of genes including virulence and antimicrobial resistance genes. Although short-read sequencing technologies have enabled large-scale bacterial genomics, the resulting draft genome assemblies are often fragmented into hundreds of discrete contigs. Several tools and approaches have been developed to identify plasmid sequences in such assemblies, but require trade-off between sensitivity and specificity. Here we propose using the Kraken classifier, together with a custom Kraken database comprising known chromosomal and plasmid sequences of Klebsiella pneumoniae species complex (KpSC), to identify plasmid-derived contigs in draft assemblies. We assessed performance using Illumina-based draft genome assemblies for 82 KpSC isolates, for which complete genomes were available to supply ground truth. When benchmarked against five other classifiers (Centrifuge, RFPlasmid, mlplasmids, PlaScope and Platon), Kraken showed balanced performance in terms of overall sensitivity and specificity (90.8 and 99.4 %, respectively, for contig count; 96.5 and >99.9 %, respectively, for cumulative contig length), and the highest accuracy (96.8% vs 91.8-96.6% for contig count; 99.8% vs 99.0-99.7 % for cumulative contig length), and F1-score (94.5 % vs 84.5-94.1 %, for contig count; 98.0 % vs 88.9-96.7 % for cumulative contig length). Kraken also achieved consistent performance across our genome collection. Furthermore, we demonstrate that expanding the Kraken database with additional known chromosomal and plasmid sequences can further improve classification performance. Although we have focused here on the KpSC, this methodology could easily be applied to other species with a sufficient number of completed genomes.

The role of mobile genetic elements in organic micropollutant degradation during biological wastewater treatment

Wastewater treatment plants (WWTPs) are crucial for producing clean effluents from polluting sources such as hospitals, industries, and municipalities. In recent decades, many new organic compounds have ended up in surface waters in concentrations that, while very low, cause (chronic) toxicity to countless organisms. These organic micropollutants (OMPs) are usually quite recalcitrant and not sufficiently removed during wastewater treatment. Microbial degradation plays a pivotal role in OMP conversion. Microorganisms can adapt their metabolism to the use of novel molecules via mutations and rearrangements of existing genes in new clusters. Many catabolic genes have been found adjacent to mobile genetic elements (MGEs), which provide a stable scaffold to host new catabolic pathways and spread these genes in the microbial community. These mobile systems could be engineered to enhance OMP degradation in WWTPs, and this review aims to summarize and better understand the role that MGEs might play in the degradation and wastewater treatment process. Available data about the presence of catabolic MGEs in WWTPs are reviewed, and current methods used to identify and measure MGEs in environmental samples are critically evaluated. Finally, examples of how these MGEs could be used to improve micropollutant degradation in WWTPs are outlined. In the near future, advances in the use of MGEs will hopefully enable us to apply selective augmentation strategies to improve OMP conversion in WWTPs.

Platon: identification and characterization of bacterial plasmid contigs in short-read draft assemblies exploiting protein sequence-based replicon distribution scores

Plasmids are extrachromosomal genetic elements that replicate independently of the chromosome and play a vital role in the environmental adaptation of bacteria. Due to potential mobilization or conjugation capabilities, plasmids are important genetic vehicles for antimicrobial resistance genes and virulence factors with huge and increasing clinical implications. They are therefore subject to large genomic studies within the scientific community worldwide. As a result of rapidly improving next-generation sequencing methods, the quantity of sequenced bacterial genomes is constantly increasing, in turn raising the need for specialized tools to (i) extract plasmid sequences from draft assemblies, (ii) derive their origin and distribution, and (iii) further investigate their genetic repertoire. Recently, several bioinformatic methods and tools have emerged to tackle this issue; however, a combination of high sensitivity and specificity in plasmid sequence identification is rarely achieved in a taxon-independent manner. In addition, many software tools are not appropriate for large high-throughput analyses or cannot be included in existing software pipelines due to their technical design or software implementation. In this study, we investigated differences in the replicon distributions of protein-coding genes on a large scale as a new approach to distinguish plasmid-borne from chromosome-borne contigs. We defined and computed statistical discrimination thresholds for a new metric: the replicon distribution score (RDS), which achieved an accuracy of 96.6 %. The final performance was further improved by the combination of the RDS metric with heuristics exploiting several plasmid-specific higher-level contig characterizations. We implemented this workflow in a new high-throughput taxon-independent bioinformatics software tool called Platon for the recruitment and characterization of plasmid-borne contigs from short-read draft assemblies. Compared to PlasFlow, Platon achieved a higher accuracy (97.5 %) and more balanced predictions (F1=82.6 %) tested on a broad range of bacterial taxa and better or equal performance against the targeted tools PlasmidFinder and PlaScope on sequenced Escherichia coli isolates. Platon is available at: http://platon.computational.bio/.

Phylogenetic background and habitat drive the genetic diversification of Escherichia coli

Escherichia coli is mostly a commensal of birds and mammals, including humans, where it can act as an opportunistic pathogen. It is also found in water and sediments. We investigated the phylogeny, genetic diversification, and habitat-association of 1,294 isolates representative of the phylogenetic diversity of more than 5,000 isolates from the Australian continent. Since many previous studies focused on clinical isolates, we investigated mostly other isolates originating from humans, poultry, wild animals and water. These strains represent the species genetic diversity and reveal widespread associations between phylogroups and isolation sources. The analysis of strains from the same sequence types revealed very rapid change of gene repertoires in the very early stages of divergence, driven by the acquisition of many different types of mobile genetic elements. These elements also lead to rapid variations in genome size, even if few of their genes rise to high frequency in the species. Variations in genome size are associated with phylogroup and isolation sources, but the latter determine the number of MGEs, a marker of recent transfer, suggesting that gene flow reinforces the association of certain genetic backgrounds with specific habitats. After a while, the divergence of gene repertoires becomes linear with phylogenetic distance, presumably reflecting the continuous turnover of mobile element and the occasional acquisition of adaptive genes. Surprisingly, the phylogroups with smallest genomes have the highest rates of gene repertoire diversification and fewer but more diverse mobile genetic elements. This suggests that smaller genomes are associated with higher, not lower, turnover of genetic information. Many of these genomes are from freshwater isolates and have peculiar traits, including a specific capsule, suggesting adaptation to this environment. Altogether, these data contribute to explain why epidemiological clones tend to emerge from specific phylogenetic groups in the presence of pervasive horizontal gene transfer across the species.

Previous large scale studies on the evolution of E. coli focused on clinical isolates emphasizing virulence and antibiotic resistance in medically important lineages. Yet, most E. coli strains are either human commensals or not associated with humans at all. Here, we analyzed a large collection of non-clinical isolates of the species to assess the mechanisms of gene repertoire diversification in the light of isolation sources and phylogeny. We show that gene repertoires evolve so rapidly by the high turnover of mobile genetic elements that epidemiologically indistinguishable strains can be phenotypically extremely heterogeneous, illustrating the velocity of bacterial adaptation and the importance of accounting for the information on the whole genome at the epidemiological scale. Phylogeny and habitat shape the genetic diversification of E. coli to similar extents. Surprisingly, freshwater strains seem specifically adapted to this environment, breaking the paradigm that E. coli environmental isolates are systematically fecal contaminations. As a consequence, the evolution of this species is also shaped by environmental habitats, and it may diversify by acquiring genes and mobile elements from environmental bacteria (and not just from gut bacteria). This may facilitate the acquisition of virulence factors and antibiotic resistance in the strains that become pathogenic.

Differential dynamics and impacts of prophages and plasmids on the pangenome and virulence factor repertoires of Shiga toxin-producing Escherichia coli O145:H28

Phages and plasmids play important roles in bacterial evolution and diversification. Although many draft genomes have been generated, phage and plasmid genomes are usually fragmented, limiting our understanding of their dynamics. Here, we performed a systematic analysis of 239 draft genomes and 7 complete genomes of Shiga toxin (Stx)-producing Escherichia coli O145:H28, the major virulence factors of which are encoded by prophages (PPs) or plasmids. The results indicated that PPs are more stably maintained than plasmids. A set of ancestrally acquired PPs was well conserved, while various PPs, including Stx phages, were acquired by multiple sublineages. In contrast, gains and losses of a wide range of plasmids have frequently occurred across the O145:H28 lineage, and only the virulence plasmid was well conserved. The different dynamics of PPs and plasmids have differentially impacted the pangenome of O145:H28, with high proportions of PP- and plasmid-associated genes in the variably present and rare gene fractions, respectively. The dynamics of PPs and plasmids have also strongly impacted virulence gene repertoires, such as the highly variable distribution of stx genes and the high conservation of a set of type III secretion effectors, which probably represents the core effectors of O145:H28 and the genes on the virulence plasmid in the entire O145:H28 population. These results provide detailed insights into the dynamics of PPs and plasmids, and show the application of genomic analyses using a large set of draft genomes and appropriately selected complete genomes.

mcr-9, an Inducible Gene Encoding an Acquired Phosphoethanolamine Transferase in Escherichia coli, and Its Origin

The plasmid-located mcr-9 gene, encoding a putative phosphoethanolamine transferase, was identified in a colistin-resistant human fecal Escherichia coli strain belonging to a very rare phylogroup, the D-ST69-O15:H6 clone. This MCR-9 protein shares 33% to 65% identity with the other plasmid-encoded MCR-type enzymes identified (MCR-1 to -8) that have been found as sources of acquired resistance to polymyxins in Enterobacteriaceae Analysis of the lipopolysaccharide of the MCR-9-producing isolate revealed a function similar to that of MCR-1 by adding a phosphoethanolamine group to lipid A and subsequently modifying the structure of the lipopolysaccharide. However, a minor impact on susceptibility to polymyxins was noticed once the mcr-9 gene was cloned and produced in an E. coli K-12-derived strain. Nevertheless, we showed here that subinhibitory concentrations of colistin induced the expression of the mcr-9 gene, leading to increased MIC levels. This inducible expression was mediated by a two-component regulatory system encoded by the qseC and qseB genes located downstream of mcr-9 Genetic analysis showed that the mcr-9 gene was carried by an IncHI2 plasmid. In silico analysis revealed that the plasmid-encoded MCR-9 shared significant amino acid identity (ca. 80%) with the chromosomally encoded MCR-like proteins from Buttiauxella spp. In particular, Buttiauxella gaviniae was found to harbor a gene encoding MCR-BG, sharing 84% identity with MCR-9. That gene was neither expressed nor inducible in its original host, which was fully susceptible to polymyxins. This work showed that mcr genes may circulate silently and remain undetected unless induced by colistin.

The Arginine Deiminase Operon Is Responsible for a Fitness Trade-Off in Extended-Spectrum-β-Lactamase-Producing Strains of Escherichia coli

We previously identified an operon involved in an arginine deiminase (ADI) pathway (arc operon) on a CTX-M-producing plasmid from an O102-ST405 strain of Escherichia coli As the ADI pathway was shown to be involved in the virulence of various Gram-positive bacteria, we tested whether the ADI pathway could be involved in the epidemiological success of extended-spectrum-β-lactamase (ESBL)-producing E. coli strains. We studied two collections of human E. coli isolated in France (n = 493) and England (n = 1,509) and show that the prevalence of the arc operon (i) is higher in ESBL-producing strains (12.1%) than in nonproducers (2.5%), (ii) is higher in CTX-M-producing strains (16%) than in other ESBL producers (3.5%), and (iii) increased over time in ESBL-producing strains from 0% before 2000 to 43.3% in 2011 to 2012. The arc operon, found in strains from various phylogenetic backgrounds, is carried by IncF plasmids (85%) or chromosomes (15%) in regions framed by numerous insertion sequences, indicating multiple arrivals. Competition experiments showed that the arc operon enhances fitness of the strain in vitro in lysogeny broth with arginine. In vivo competition experiments showed that the arc operon is advantageous for the strain in a mouse model of urinary tract infection (UTI), whereas it is a burden in a mouse model of intestinal colonization. In summary, we have identified a trait linked to CTX-M-producing strains that is responsible for a trade-off between two main E. coli lifestyles, UTI and gut commensalism. This trait alone cannot explain the wide spread of ESBLs in E. coli but merits epidemiological surveillance.