PubMed Health⌕ Search

Biomedical subjects

Edward J Feil

Publications and source records attributed to Edward J Feil.

At least 19 recordsLinked to original sources

Convergence and global molecular epidemiology of Klebsiella pneumoniae plasmids harbouring the iuc3 virulence locus: a population genomic analysis.

BACKGROUND: Klebsiella pneumoniae is an important pathogen of humans and animals. In the past five years, increasing reports of convergent strains that carry both virulence factors and antimicrobial resistance genes (ARGs) have raised serious public health concerns. The aim of this study is to describe the global diversity of plasmids carrying iuc3 (a key virulence factor in K pneumoniae associated with pigs and clinical isolates) from diverse settings, and their role in the emergence of convergent strains through hybridisation with plasmids carrying ARGs. METHODS: This population genomic analysis study was designed to describe both the global and local diversity of iuc3-carrying plasmids from diverse sources, and the co-occurrence of iuc3 with ARGs. We used all 4148 Klebsiella spp isolates from two large One-Health studies (SpARK, Italy, and OH-DART, Thailand), including 191 Klebsiella isolates from pigs, 635 from clinical isolates, 1040 from hospital and community carriage, and 2282 from other sources. Short-read sequencing of Klebsiella isolates was performed as part of the SpARK study. We sequenced Klebsiella isolates from the OH-DART (MicrobesNG, Birmingham, UK; HiSeq and NovaSeq, Illumina San Diego, CA, USA; GridION, Oxford Nanopore Technologies, Oxford, UK) and SpARK (MinION or GridION, Oxford Nanopore Technologies, Oxford, UK) studies. We also retrieved plasmid sequences carrying iuc3 from the National Centre for Biotechnology Information (NCBI). To ascertain the degree of diversity, evolutionary dynamics, and structuring across ecological and geographical axes, we detected ARGs and virulence loci, analysed clustering patterns and generated approximate maximum-likelihood phylogenetic trees. FINDINGS: We identified 48 K pneumoniae isolates with iuc3 in the SpARK data and 79 in the OH-DART data. Three (2·4%) of these 127 isolates were from clinical sources, 73 (57·5%) were from pig or pork meat. iuc3 isolates corresponded to multiple (n=47) host sequence types (STs), with ST35, ST45, ST881, ST25, and ST967 harbouring iuc3 in both datasets. We generated hybrid assemblies for 44 (SpARK) and 36 (OH-DART) isolates, plus a single iuc3 isolate from Germany. 53 (65·4%) of these isolates were from pigs, three (3·7%) from clinical sources, and 25 (30·9%) from other sources. There were an additional 48 iuc3 positive isolates from our collections for which only short read data was available. A single iuc3-positive Klebsiella oxytoca isolate from a pig farm was detected in the SpARK data, which was also sequenced. We identified 330 iuc3-positive isolates and 58 iuc3-carrying plasmid assemblies from NCBI, of which 83 (21·4%) were from clinical sources, 120 from pigs (30·9%), and 185 (47·7%) from other sources or of unknown provenance. These isolates were from K pneumoniae except two isolates of Klebsiella quasipneumoniae subsp similipneumoniae and one of Enterobacter hormaechei. The combined dataset of 517 iuc3 plasmids ranged in size from 110 375 bp to 365 580 bp and mostly corresponded to multiple IncFIB(K) and IncFII replicon types. We found seven convergent K pneumoniae plasmids in the Thai data: six from fresh markets and one from a neighbouring hospital. These plasmids emerged through the hybridisation of cocirculating iuc3 plasmids and plasmids encoding extended-spectrum β-lactamases (ESBLs), although none of these seven plasmids carried genes encoding carbapenemases. We also identified putative cocirculating parental plasmids carrying iuc3 and ESBL-encoding genes. Clustering and phylogenetic analysis resolved the iuc3 plasmid sequences into three groups, which were consistent using both complete plasmid sequences (n=139) and short-read data (n=517). In the complete plasmid sequence data, 66 strains contained group 1 plasmids, 38 strains contained group 2 plasmids, and 35 strains contained group 3 plasmids. Group 3 plasmids are mostly carried by isolates circulating in hospitals throughout Asia, with occasional examples in Europe and elsewhere, and carry multiple ARGs and potential virulence factors. By contrast, group 1 plasmids are commonly carried by porcine isolates in Europe, and group 2 are a heterogeneous mixture of geographical and ecological sources. INTERPRETATION: Plasmid hybridisation occurs frequently outside of the health-care environment and can lead to the convergence of resistance and virulence traits. Generating complete plasmid sequences from regional population-scale samples facilitates the identification of convergent plasmids and their putative parental plasmids. Three robust groups of iuc3 plasmids were resolved, which show both epidemiological and geographical differences; one of these groups was associated with clinical isolates in Asia and warrants targeted plasmid surveillance. FUNDING: UKRI, JPIAMR, Evolution Education Trust, and a Schlumberger Foundation Fellowship.

Plasmids↗

Global emergence and transmission dynamics of carbapenemase-producing Citrobacter freundii sequence type 22 high-risk international clone: a retrospective, genomic, epidemiological study.

BACKGROUND: Carbapenemase-producing Citrobacter (CPC) species have recently been recognised as emerging pathogens associated with nosocomial infections in humans. The increased rate of Citrobacter freundii infections is a public health concern and there is a paucity of genomic data regarding its global transmission dynamics. We aimed to characterise the genetic features of CPC species, and their associated carbapenemase-encoding plasmids, obtained from hospitalised patients in China and from publicly available global data, with a particular focus on high-risk clones. METHODS: This was a retrospective, genomic epidemiological study of CPC species obtained from a tertiary hospital in Zhejiang Province, China, from March 5, 2013, to March 5, 2023. We used antimicrobial susceptibility testing, short-read and long-read whole-genome sequencing, phylogenomic analysis, and plasmid structure analysis. A global dataset of complete plasmid sequences encoding blaKPC, blaNDM, and blaIMP was constructed from the National Center for Biotechnology Information (NCBI) RefSeq database to provide insights into their diversity and distribution. All carbapenemase-producing Citrobacter freundii genomes from the NCBI GenBank database were incorporated in the comparative genomic analyses. Bayesian phylogeographical analysis and growth rate assays were carried out to characterise the high-risk C freundii sequence type (ST) 22 clone. FINDINGS: 1724 Citrobacter species isolates were collected from diverse clinical specimens, with 48 identified as CPC species. Citrobacter koseri (22 [46%] of 48) and C freundii (20 [42%]) were the predominant CPC species. Comparative analysis found C freundii carried significantly higher median numbers of plasmid replicons (5&#xb7;0 [IQR 3&#xb7;3-6&#xb7;0] vs 2&#xb7;0 [2&#xb7;0-3&#xb7;0]; p<0&#xb7;0001) and acquired antimicrobial resistance genes (12&#xb7;0 [7&#xb7;3-15&#xb7;8] vs 3&#xb7;0 [3&#xb7;0-5&#xb7;3]; p<0&#xb7;0001) than did C koseri. Molecular characterisation identified Inc-type plasmids, In823::Kl.pn.I3/In1589-like/In837-like integrons, Tn6296/Tn125/Tn5060 transposons, and insertion sequences (eg, IS26, IS3000, IS5, ISAba125, ISCR1), collectively facilitating the dissemination of carbapenemase genes. Global analysis of 3126 carbapenemase-encoding plasmids found epidemic plasmids with broad host ranges and global diversity. Phylogenetic investigation of predominant carbapenemase-encoding plasmids showed their persistence across geographical regions, temporal spans, and Enterobacterales species, exhibiting high genetic similarity to our clinical plasmids. A phylogenetic tree of 726 global carbapenemase-producing C freundii genomes showed that ST22 (227 [31&#xb7;3%]) represents the predominant multidrug-resistant clone across community, health-care, and environmental niches. Transmission across continents contributes to the global predominance of the ST22 clone, which carries a high load of resistance genes (median 15&#xb7;0 [IQR 11&#xb7;0-17&#xb7;0] vs 12&#xb7;0 [3&#xb7;0-16&#xb7;0]; p<0&#xb7;0001) and enhanced plasmid maintenance capacity (median replicons 5&#xb7;0 [IQR 4&#xb7;0-7&#xb7;0] vs 4&#xb7;0 [3&#xb7;0-6&#xb7;0]; p<0&#xb7;0001) relative to non-ST22 clones. INTERPRETATION: Our study provides evidence to suggest that Citrobacter species are emerging carriers of carbapenem-resistance genes. These findings provide insight into the population structure of CPC species and highlight C freundii ST22 as a prominent high-risk international clone. FUNDING: National Natural Science Foundation of China, National Health Commission Scientific Research Fund-Zhejiang Provincial Major Health Science and Technology Plan Project, Zhejiang Province Natural Science Foundation Project, Outstanding Youth Foundation of Jiangsu Province of China, the Priority Academic Program Development of Jiangsu Higher Education Institutions, and Postgraduate Research and Practice Innovation Program of Jiangsu Province.

Citrobacter freundii↗

The secret life of the multilocus sequence type.

Such are the challenges, and the potential, presented by complete genome sequences that the eventual erosion of the boundaries between biochemistry, ecology, bioinformatics, population biology, epidemiology and medical microbiology will perhaps be the most profound legacy of the genomics revolution. The development of nucleotide sequence-based typing schemes (multilocus sequence typing (MLST)) represents a similar synthesis, for this technique both matches the practical requirements for a highly portable standard for strain characterisation whilst also being firmly grounded in the population biology principles of multilocus enzyme electrophoresis (MLEE). Contrary to recent claims that population biology analyses of public health-oriented MLST data 'obscures its utility in applied microbiology' [Maiden MC. Multilocus sequence typing of bacteria. Annu Rev Microbiol 2006;60:561-88.], we argue that such an emphasis is essential for full interpretation of the data. Here we note a pertinent case in point; how a consideration of the rates of genetic recombination can help to explain why MLST data tend to correlate with virulence properties in some species (Neisseria meningitidis) but not in others (Staphylococcus aureus). We also discuss how the argument applies to the identification of recently emerged methicillin-resistant S. aureus (MRSA) clones using MLST. We conclude with a speculative rationale for promoting the 'clonal complexes' of S. aureus to species status.

Bacterial Typing Techniques↗

Similar compositional biases are caused by very different mutational effects.

Compositional replication strand bias, commonly referred to as GC skew, is present in many genomes of prokaryotes, eukaryotes, and viruses. Although cytosine deamination in ssDNA (resulting in C-->T changes on the leading strand) is often invoked as its major cause, the precise contributions of this and other substitution types are currently unknown. It is also unclear if the underlying mutational asymmetries are the same among taxa, are stable over time, or how closely the observed biases are to mutational equilibrium. We analyzed nearly neutral sites of seven taxa each with between three and six complete bacterial genomes, and inferred the substitution spectra of fourfold degenerate positions in nonhighly expressed genes. Using a bootstrap procedure, we extracted compositional biases associated with replication and identified the significant asymmetries. Although all taxa showed an overrepresentation of G relative to C on the leading strand (and imbalances between A and T), widely variable substitution asymmetries are noted. Surprisingly, all substitution types show significant asymmetry in at least one taxon, but none were universally biased in all taxa. Notably, in the two most biased genomes, A-->G, rather than C-->T, shapes the compositional bias. Given the variability in these biases, we propose that the process is multifactorial. Finally, we also find that most genomes are not at compositional equilibrium, and suggest that mutational-based heterotachy is deeply imprinted in the history of biological macromolecules. This shows that similar compositional biases associated with the same essential well-conserved process, replication, do not reflect similar mutational processes in different genomes, and that caution is required in inferring the roles of specific mutational biases on the basis of contemporary patterns of sequence composition.

Bacteria↗

Protein evolution: causes of trends in amino-acid gain and loss.

Understanding how proteins evolve is important for determining the molecular basis of adaptation, for inferring phylogenies and for engineering novel proteins. It has been suggested that some amino acids were incorporated into the genetic code more recently than others and, after comparing pairs of closely related genomes, Jordan et al. report that 'recent' amino acids are becoming more common. They argue that this process has been going on since the genetic code first evolved to encompass all 20 amino acids. Here we provide evidence that the patterns observed conform with standard, nearly neutral theoretical expectations and require no new explanation. This reinforces the need for caution in the interpretation of results derived from closely related taxa.

Amino Acids↗

The phylogeny of Staphylococcus aureus - which genes make the best intra-species markers?

The ability to make informed decisions on the suitability of alternative marker loci is central for population and epidemiological investigations. This issue was addressed using Staphylococcus aureus as a model population by generating nucleotide sequence data from 33 gene fragments in a representative sample of 30 strains. Supplementing the data with pre-existing multilocus sequence typing data, an intra-species tree based on approximately 17.8 kb of sequence was reconstructed and the goodness of fit of each individual gene tree was computed. No strong association was noted between gene function per se and phylogenetic reliability, but it is suggested that candidate loci should possess at least the average degree of nucleotide diversity for all genes in the genome. In the case of S. aureus this threshold is >1 % mean pairwise diversity.

DNA, Bacterial↗

Nonrandom distribution of Burkholderia pseudomallei clones in relation to geographical location and virulence.

Burkholderia pseudomallei is a soil-dwelling saprophyte and the causative agent of melioidosis, a life-threatening human infection. Most cases are reported from northeast Thailand and northern Australia. Using multilocus sequence typing (MLST), we have compared (i) soil and invasive isolates from northeast Thailand and (ii) invasive isolates from Thailand and Australia. A total of 266 Thai B. pseudomallei isolates were characterized (83 soil and 183 invasive). These corresponded to 123 sequence types (STs), the most abundant being ST70 (n=21), ST167 (n=15), ST54 (n=12), and ST58 (n=11). Two clusters of related STs (clonal complexes) were identified; the larger clonal complex (CC48) did not conform to a simple pattern of radial expansion from an assumed ancestor, while a second (CC70) corresponded to a simple radial expansion from ST70. Despite the large number of STs, overall nucleotide diversity was low. Of the Thai isolates, those isolated from patients with melioidosis were overrepresented in the 10 largest clones (P<0.0001). There was a significant difference in the classification index between environmental and disease isolates (P<0.001), confirming that genotypes were not distributed randomly between the two samples. MLST profiles for 158 isolates from Australia (mainly disease associated) contained a number of STs (96) similar to that seen with the Thai invasive isolates, but no ST was found in both populations. There were also differences in diversity and allele frequency distribution between the two populations. This analysis reveals strong genetic differentiation on the basis of geographical isolation and a significant differentiation on the basis of virulence potential.

Australia↗

Nematode symbiont for Photorhabdus asymbiotica.

Photorhabdus asymbiotica is an emerging bacterial pathogen that causes locally invasive soft tissue and disseminated bacteremic infections in the United States and Australia. Although the source of infection was previously unknown, we report that the bacterium is found in a symbiotic association with an insect-pathogenic soil nematode of the genus Heterorhabditis.

Animals↗

Comparisons of dN/dS are time dependent for closely related bacterial genomes.

The ratio of non-synonymous (dN) to synonymous (dS) changes between taxa is frequently computed to assay the strength and direction of selection. Here we note that for comparisons between closely related strains and/or species a second parameter needs to be considered, namely the time since divergence of the two sequences under scrutiny. We demonstrate that a simple time lag model provides a general, parsimonious explanation of the extensive variation in the dN/dS ratio seen when comparing closely related bacterial genomes. We explore this model through simulation and comparative genomics, and suggest a role for hitch-hiking in the accumulation of non-synonymous mutations. We also note taxon-specific differences in the change of dN/dS over time, which may indicate variation in selection, or in population genetics parameters such as population size or the rate of recombination. The effect of comparing intra-species polymorphism and inter-species substitution, and the problems associated with these concepts for asexual prokaryotes, are also discussed. We conclude that, because of the critical effect of time since divergence, inter-taxa comparisons are only possible by comparing trajectories of dN/dS over time and it is not valid to compare taxa on the basis of single time points.

Bacteria↗

Opinion: Re-evaluating prokaryotic species.

There is no widely accepted concept of species for prokaryotes, and assignment of isolates to species is based on measures of phenotypic or genome similarity. The current methods for defining prokaryotic species are inadequate and incapable of keeping pace with the levels of diversity that are being uncovered in nature. Prokaryotic taxonomy is being influenced by advances in microbial population genetics, ecology and genomics, and by the ease with which sequence data can be obtained. Here, we review the classical approaches to prokaryotic species definition and discuss the current and future impact of multilocus nucleotide-sequence-based approaches to prokaryotic systematics. We also consider the potential, and difficulties, of assigning species status to biologically or ecologically meaningful sequence clusters.

Bacteria↗

Databases and software for the comparison of prokaryotic genomes.

The explosion in the number of complete genomes over the past decade has spawned a new and exciting discipline, that of comparative genomics. To exploit the full potential of this approach requires the development of novel algorithms, databases and software which are sophisticated enough to draw meaningful comparisons between complete genome sequences and are widely accessible to the scientific community at large. This article reviews progress towards the development of computational tools and databases for organizing and extracting biological meaning from the comparison of large collections of genomes.

Computational Biology↗

Evolutionary genetics of the accessory gene regulator (agr) locus in Staphylococcus aureus.

The accessory gene regulator (agr) locus influences the expression of many virulence genes in the human pathogen Staphylococcus aureus. Four allelic groups of agr, which generally inhibit the regulatory activity of each other, have been identified within the species. Interference in virulence gene expression caused by different agr groups has been suggested to be a mechanism for isolating bacterial populations and a fundamental basis for subdividing the species. To test the hypothesis that the species is phylogenetically structured according to agr groups, we mapped agr groups onto a clone phylogeny inferred from partial sequences of 14 genes from 27 genetically diverse strains. Shimodaira-Hasegawa and parametric bootstrap tests rejected the hypotheses that the species is subdivided into three or five monophyletic agr groups but failed to reject the hypothesis that the species is subdivided into two groups that each consist of multiple clonal complexes and multiple agr groups. Additional evidence for agr recombination is found from clustered polymorphisms in complete agr sequences. However, agr recombination has not occurred frequently or randomly through time, because the topology and branch lengths of the clone phylogeny are reflected within each agr group. To account for these observations, we propose a new evolutionary model that involves a genetically polymorphic ancestral population of S. aureus that horizontally transferred agr groups between two subspecies groups near the time that these subspecies groups diverged.

Bacterial Proteins↗

Displaying the relatedness among isolates of bacterial species -- the eBURST approach.

Determining the most appropriate way to represent the relationships between bacterial isolates is complicated by the differing rates of recombination within species. In many cases, a bifurcating tree can be positively misleading. The recently described program eBURST can be used with multilocus data to define groups or clonal complexes of related isolates derived from a common ancestor, the patterns of descent linking them together, and the ancestral genotype. eBURST has recently been extensively updated to include additional tools for exploring the relationships between isolates. We discuss the advantages of this approach and describe its use to explore patterns of descent within clonal complexes identified using multilocus sequence typing.

Algorithms↗

Complete genomes of two clinical Staphylococcus aureus strains: evidence for the rapid evolution of virulence and drug resistance.

Staphylococcus aureus is an important nosocomial and community-acquired pathogen. Its genetic plasticity has facilitated the evolution of many virulent and drug-resistant strains, presenting a major and constantly changing clinical challenge. We sequenced the approximately 2.8-Mbp genomes of two disease-causing S. aureus strains isolated from distinct clinical settings: a recent hospital-acquired representative of the epidemic methicillin-resistant S. aureus EMRSA-16 clone (MRSA252), a clinically important and globally prevalent lineage; and a representative of an invasive community-acquired methicillin-susceptible S. aureus clone (MSSA476). A comparative-genomics approach was used to explore the mechanisms of evolution of clinically important S. aureus genomes and to identify regions affecting virulence and drug resistance. The genome sequences of MRSA252 and MSSA476 have a well conserved core region but differ markedly in their accessory genetic elements. MRSA252 is the most genetically diverse S. aureus strain sequenced to date: approximately 6% of the genome is novel compared with other published genomes, and it contains several unique genetic elements. MSSA476 is methicillin-susceptible, but it contains a novel Staphylococcal chromosomal cassette (SCC) mec-like element (designated SCC(476)), which is integrated at the same site on the chromosome as SCCmec elements in MRSA strains but encodes a putative fusidic acid resistance protein. The crucial role that accessory elements play in the rapid evolution of S. aureus is clearly illustrated by comparing the MSSA476 genome with that of an extremely closely related MRSA community-acquired strain; the differential distribution of large mobile elements carrying virulence and drug-resistance determinants may be responsible for the clinically important phenotypic differences in these strains.

Drug Resistance, Bacterial↗

Analyses of clonality and the evolution of bacterial pathogens.

The existence of bacterial clones was evident in early phenotypic studies that recognised high levels of similarity in geographically and temporally separated isolates. Multilocus sequence typing (MLST) has become the most common method for genetically characterizing clones of several bacterial pathogens, allowing the tracking of hypervirulent/antibiotic-resistant lineages. MLST has also been used to examine the way that bacterial populations, and in particular, bacterial clones evolve. Visualisation of MLST datasets has required the development of novel tools, such as 'eBURST', a key program in constructing evolutionary models that detail how methicillin resistant Staphylococcus aureus (MRSA) and other clones emerge and spread.

Algorithms↗

Multilocus sequence typing--what is resolved?

Nucleotide sequence-based methods for bacterial typing (multilocus sequence typing; MLST) allow rapid and global comparisons between results from different laboratories. Combining this advantage with the reduced cost of high throughput sequencing, increasing automation and the amenability of sequence data for evolutionary analysis, it seems inevitable that sequence-based typing will eventually predominate over gel-based methods such as pulsed-field gel electrophoresis (PFGE) for most bacterial species. The increasing availability of multiple genome sequences for single pathogenic species, and the recent development of many new MLST schemes, means that a re-examination of the utility of multilocus sequencing, and in particular the choice of gene loci, is now appropriate.

Bacteria↗

eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data.

The introduction of multilocus sequence typing (MLST) for the precise characterization of isolates of bacterial pathogens has had a marked impact on both routine epidemiological surveillance and microbial population biology. In both fields, a key prerequisite for exploiting this resource is the ability to discern the relatedness and patterns of evolutionary descent among isolates with similar genotypes. Traditional clustering techniques, such as dendrograms, provide a very poor representation of recent evolutionary events, as they attempt to reconstruct relationships in the absence of a realistic model of the way in which bacterial clones emerge and diversify to form clonal complexes. An increasingly popular approach, called BURST, has been used as an alternative, but present implementations are unable to cope with very large data sets and offer crude graphical outputs. Here we present a new implementation of this algorithm, eBURST, which divides an MLST data set of any size into groups of related isolates and clonal complexes, predicts the founding (ancestral) genotype of each clonal complex, and computes the bootstrap support for the assignment. The most parsimonious patterns of descent of all isolates in each clonal complex from the predicted founder(s) are then displayed. The advantages of eBURST for exploring patterns of evolutionary descent are demonstrated with a number of examples, including the simple Spain(23F)-1 clonal complex of Streptococcus pneumoniae, "population snapshots" of the entire S. pneumoniae and Staphylococcus aureus MLST databases, and the more complicated clonal complexes observed for Campylobacter jejuni and Neisseria meningitidis.

Algorithms↗