Next Article in Journal
Association of Long Non-Coding RNA Polymorphisms with Gastric Cancer and Atrophic Gastritis
Next Article in Special Issue
Carriage and Gene Content Variability of the pESI-Like Plasmid Associated with Salmonella Infantis Recently Established in United States Poultry Production
Previous Article in Journal
Seascape Genomics of the Sugar Kelp Saccharina latissima along the North Eastern Atlantic Latitudinal Gradient
Previous Article in Special Issue
Genotypic and Phenotypic Characterization of Incompatibility Group FIB Positive Salmonella enterica Serovar Typhimurium Isolates from Food Animal Sources
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 12
10.3390/genes11121504
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genomic Investigation into the Virulome, Pathogenicity, Stress Response Factors, Clonal Lineages, and Phylogenetic Relationship of Escherichia coli Strains Isolated from Meat Sources in Ghana

by
Frederick Adzitey
1,*,
Jonathan Asante
2,
Hezekiel M. Kumalo
3,
Rene B. Khan
3,
Anou M. Somboro
3 and
Daniel G. Amoako
2,*
1
Department of Veterinary Science, University for Development Studies, P. O. Box TL 1882, Tamale 00233, Ghana
2
Antimicrobial Research Unit, Infection Genomics and Applied Bioinformatics Division, University of KwaZulu-Natal, Durban 4000, South Africa
3
School of Laboratory Medicine and Medical Science, University of KwaZulu-Natal, Durban 4000, South Africa
*
Authors to whom correspondence should be addressed.
Genes 2020, 11(12), 1504; https://doi.org/10.3390/genes11121504
Submission received: 5 November 2020 / Revised: 29 November 2020 / Accepted: 3 December 2020 / Published: 14 December 2020
(This article belongs to the Special Issue Genetics and Genomics of Zoonotic Foodborne Pathogens)
Browse Figures
Versions Notes

Abstract

:
Escherichia coli are among the most common foodborne pathogens associated with infections reported from meat sources. This study investigated the virulome, pathogenicity, stress response factors, clonal lineages, and the phylogenomic relationship of E. coli isolated from different meat sources in Ghana using whole-genome sequencing. Isolates were screened from five meat sources (beef, chevon, guinea fowl, local chicken, and mutton) and five areas (Aboabo, Central market, Nyorni, Victory cinema, and Tishegu) based in the Tamale Metropolis, Ghana. Following microbial identification, the E. coli strains were subjected to whole-genome sequencing. Comparative visualisation analyses showed different DNA synteny of the strains. The isolates consisted of diverse sequence types (STs) with the most common being ST155 (n = 3/14). Based Upon Related Sequence Types (eBURST) analyses of the study sequence types identified four similar clones, five single-locus variants, and two satellite clones (more distantly) with global curated E. coli STs. All the isolates possessed at least one restriction-modification (R-M) and CRISPR defence system. Further analysis revealed conserved stress response mechanisms (detoxification, osmotic, oxidative, and periplasmic stress) in the strains. Estimation of pathogenicity predicted a higher average probability score (Pscore ≈ 0.937), supporting their pathogenic potential to humans. Diverse virulence genes that were clonal-specific were identified. Phylogenomic tree analyses coupled with metadata insights depicted the high genetic diversity of the E. coli isolates with no correlation with their meat sources and areas. The findings of this bioinformatic analyses further our understanding of E. coli in meat sources and are broadly relevant to the design of contamination control strategies in meat retail settings in Ghana.

1. Introduction

Escherichia coli (E. coli) are Gram-negative, rod shape, facultative anaerobic bacteria of the Enterobacteriaceae family [1,2]. They are naturally found in the gastrointestinal tract of animals and their droppings [2,3], and can cross-contaminate meat, especially under faulty processing conditions [4]. E. coli are classified into two major groups, namely, intestinal or diarrheagenic E. coli (DEC) and extraintestinal E. coli. DEC is divided into six pathotypes, namely, Enteropathogenic E. coli, Enterohaemorrhagic E. coli, Shiga toxin producing E. coli, Enteroinvasive E. coli, Enteroaggregative E. coli, Enterotoxigenic E. coli, Verotoxigenic E. coli, and diffusely adherent E. coli [5]. The extraintestinal E. coli group includes subtypes, such as neonatal meningitis associated E. coli, uropathogenic E. coli, and sepsis-associated E. coli [6,7]. Their various virulence genes differentiate these pathotypes. Their pathogenesis is initiated by the adherence to the host’s epithelial cell, and, then colonization with some pathotypes invaded the cells, thereby establishing infections [8].
E. coli survive by overcoming stressors released by their host or substrate [9]. In meats, E. coli must overcome stressors such as heat, especially during cold storage under refrigerated conditions, starvation, and acid when meats are treated with ingredients/spices for preservation and cooking. Stressors result in the release of regulator proteins and genes by E. coli to confer resistance and to promote survival [10,11]. This resistance and tolerance can lead to enhanced virulence.
E. coli of meat origin harbouring virulence genes have been reported [12,13,14]. Virulence genes can increase the pathogenicity of E. coli by means such as inhibition/evasion of the immune system of the host and obtaining its nutrition from the host, thereby, depriving the host of some nutrients [15]. Caruso et al. [14] found that 28% of extraintestinal E. coli of meat origin harboured virulence genes, which included kpsMII, iutA, and PapA genes. Lyhs et al. [13] characterised 207 E. coli isolates from poultry meat products and reported that a virulent gene was present in each isolate. In addition, all E. coli strains from raw meat and shell fish contained at least one virulence gene of the 16 genes detected [12].
Meat is any edible part of a slaughtered animal. It can be obtained from cattle, goat, sheep, poultry, and pigs, among other animals. Meat is an important source of protein and other nutrients for humans worldwide [16,17]. In Ghana, meat is an essential component of the diet of most people [18] and its association with pathogens represent a risk to public health. The most frequent mode of transmission for E. coli to man is through the consumption of food and water [19]. Meat has been implicated in several foodborne outbreak infections. Outbreaks of E. coli linked to ground beef and poultry have also been reported globally [20,21,22].
Molecular investigations have been carried using a number of typing techniques including enterobacterial repetitive intergenic consensus (ERIC), random amplified polymorphic deoxyribonucleic acid (RAPD), repetitive extragenic palindromic (REP), pulse-field gel electrophoresis, and multi-locus sequence typing [23,24,25,26,27]. Recent advances have promoted the use of whole-genome sequencing and metagenomics. Nonetheless, whole-genome sequencing (WGS) of pathogens has become the reference standard due to its accessibility and affordability, and has revolutionised the field of outbreak investigations [28]. WGS enables the genome of a bacterium to be sequenced and helps in the in silico prediction of various traditional typing methods, which explore short DNA strands with a quick turnaround time comparatively [29]. Reports on the genomic characterisation of E. coli from meat sources in Ghana are limited. The study used whole-genome sequencing to identify virulome, pathogenicity, stress response factors, clonal lineages, and phylogenomic relationship of E. coli isolated from various meat types in Ghana. This will provide useful information on the risk factors related to meat consumption and adaptive characteristics of E. coli that favour their clonal spread and survival of various niches.

2. Materials and Methods

2.1. Sample Collection and Identification of E. coli

2.1.1. Sample Collection

A total of two hundred and twenty-five (225) meat samples comprising of beef (n = 45), chevon (n = 45), chicken (n = 45), guinea fowl (n = 45), and mutton (n = 45) were sampled in Tamale (Ghana) between April and December 2016. Sterile cotton swabs were used to swab an area of 10 cm2 of each meat sample. The swabs were transported under 4 °C and analysed immediately upon reaching the laboratory.

2.1.2. Isolation of E. coli

The procedure in the Food and Drug Administration-Bacteriological Analytical Manual [2] with slight modification was used. The swabs were dipped into 10 mL of Buffered Peptone Water and incubated at 37 °C for 24 h. Afterward, 0.1 mL of the aliquots were streaked on Levine’s Eosin-methylene Blue Agar and incubated at 37 °C for 24 h. Presumptive E. coli colonies appeared as dark-centered and flat with or without metallic sheen. Presumptive E. coli colonies were purified on Trypticase Soy Agar and incubated at 37 °C for 24 h. They were identified and confirmed using Gram staining, growth on MacConkey Agar, growth in Brilliant Green Bile Broth, and the E. coli latex agglutination test. All media and reagents used were purchased from Oxoid Limited, Basingstoke, UK. A total of 189 E. coli isolates were recovered from the various meat sources (Table S1).

2.2. Selection of Isolates for Genome Sequencing and Assembly

A combination of cluster and simple random sampling were used to ensure that 14 isolates from the different meat sources and areas were chosen for WGS. A Genomic DNA (gDNA) of the E. coli isolates was extracted and purified using the QIAampR DNA Mini Kit (QIAGEN, Hilden, Germany). Following extraction, DNA quantification and the quality check was performed on a Qubit® 2.0 fluorometer (Life Technologies, Carlsbad, CA, USA) and Nanodrop 8000 (Thermo Scientific, Waltham, MA, USA), respectively [30]. The Nextera XT DNA Sample Preparation Kit was used to prepare the paired-end library and WGS was performed using the Illumina MiSeq machine (Illumina, San Diego, CA, USA). The obtained raw reads were de-novo assembled with the Shovill pipeline version 0.9.0 that uses SPAdes version 3.11.0. The assembled contigs were deposited in GenBank under project number PRJNA484345.

2.3. Genome Annotation and Visualisation

The Pathogen watch platform was used to confirm the identity of the E. coli isolates generated from the WGS data [31]. The genomes of the strains were visualised using the GView Server [32]. The National Center for Biotechnology Information (NCBI) Prokaryotic Genome Annotation Pipeline (PGAP; version 4.3) [33] and Rapid Annotation using Subsystem Technology (RAST) Server (version 2.0) [34] were used for annotation of the isolates.

2.4. WGS-Based Molecular Typing of Isolates

Multilocus sequence typing (MLST) v2.0 (https://cge.cbs.dtu.dk/services/MLST/) was used to predict the sequence types (STs) from the assembled genomes. An eBURST [35] analysis was done in the MLST database (https://pubmlst.org/escherichia/) to find out whether the obtained STs were satellite clones, double-locus variant (DLV), or single-locus variant (SLV) of the known STs using the Achtman scheme. The reference E. coli online platform, CHTyper v1.0 (https://cge.cbs.dtu.dk/services/CHTyper/) was used to infer the FumC and FimH types.

2.5. Detection of the Stress Response Mechanisms, CRISPR Array, and Restriction-Modification System (R-M System)

The CRISPRCasFinder available at: https://crisprcas.i2 bc.paris-saclay.fr/CrisprCasFinder/Index [36] was used to infer common CRISPR and Cas loci in the assembled genomes. The Restriction Modification Finder at: https://cge.cbs.dtu.dk/services/Restriction-ModificationFinder/ and Pathosystems Resource Integration Center (PATRIC) platform (https://www.patricbrc.org/) were used to predict the restriction modification system in the isolates.

2.6. Assessment of Pathogenic Potential

The automated mode of the PathogenFinder web service was used to estimate the pathogenic potential in the E. coli genomes [37]. All E. coli isolates were subjected to the pathogenicity prediction using fasta formatted genome data.

2.7. Virulome Analysis

Virulence factors corresponding to main bacterial virulence determinants (toxin, adherence, motility, iron uptake, and secretion system) associated with E. coli were searched with the Virulence Finder 2.0 at: https://cge.cbs.dtu.dk/services/VirulenceFinder/ [38].

2.8. Comparative Phylogenomic Analysis and Metadata Insights

The CSI Phylogeny-1.4 (https://cge.cbs.dtu.dk/services/CSIPhylogeny-1.2) was used to infer the phylogenetic relationship between E. coli genomes [39]. The genome of K. quasi-pneumoniae strain P27-02 (accession number: NXHG00000000.1) served as the outgroup to root the tree for easy configuration of the phylogenetic distance between the isolates. The pipeline was run with default parameters. The phylogeny was viewed with annotations for metadata (virulome, sequence types, meat source, and area) using Phandango [40] to provide vital insights of the phylogenetic tree.

3. Results

3.1. E. coli Strains and Data Source for Comparative Genome Analysis

A total of 14 isolates obtained from five meat sources [beef (n = 3), chevon (n = 3), guinea fowl (n = 3), local chicken (n = 3), and mutton (n = 2)] and five areas [Aboabo (n = 2), Central market (n = 3), Nyorni (n = 3), Victory cinema (n = 2), and Tishegu (n = 4)] based in the Tamale Metropolis, Ghana (Table 1) were selected for further bioinformatic and comparative analyses.

3.2. Genomic Visualisation and Annotation

The global Pathogen watch platform confirmed the phenotypic identity of the E. coli genomes. Comparative visualisation analysis via the GView server showed different DNA synteny of the strains with a coverage range of 90–97% from the reference strain (ATCC25922) in all the E. coli genomes (Figure 1).

3.3. In Silico Typing of the E. coli Isolates

Eleven different sequence types were identified in the isolates (Table 1). The predominant sequence types were the ST155 (n = 3) and ST540 (n = 2). Of note, most of the sequence types were singletons (n = 9) (Table 1). There was high diversity in the clonal lineages of strains with respect to their source and area of origin (Table 1). Based Upon Related Sequence Types (eBURST) analyses of the study sequence types identified four similar clones, five single-locus variants (SLV), and two satellite clones (more distantly) with global curated E. coli STs from animals (food), humans, and environmental sources (Figure 2 and Table S3). The E. coli strains possessed diverse serotypes with no association with the meat sources, area, and sequence types (Table 1). Of note, not all the O-antigen and H-antigen loci in the strains were predicted. The FumC and FimH types were identified in the E. coli isolates, which correlated with the sequence types to a large extent except for two isolates (TG5 and CM4) that differed in their STs and FumC/FimH types (Table 1). The FumC7/FimH54 (n = 3) and FumC4/FimH32 (n = 2) were most frequently identified. Several singletons FumC and FimH types were also identified.

3.4. Genomic Prediction of Defence Systems (CRISPR-Cas Elements and R-M System) and Stress Response Factors

The genomes of the strains contained CRISPR arrays (nrange = 3–9) with a Cas element belonging to the same subcategory (Type1 E) except for two isolates, known as CB1 (Beef from Central Market) and TG5 (Guinea fowl from Tishegu) that did not encode any Cas element (Table 1). All the isolates possessed at least one restriction modification (R-M) system within the genome. The Type II R-M system was the most common in the E. coli isolates with M.EcoKII being the most dominant (n = 12), which is followed by M.EcoJA03 PDcm (n = 7). The isolate SG6 obtained from Guinea fowl (Victory Cinema) contained the highest number of R-M systems (n = 3, Type I, II, and IV) (Table 1). Conserved stress response factors including the detoxification system, osmotic, oxidative, and periplasmic stress were found in all the isolates (Figure 3). The oxidative stress response was the most predominant mechanism (n = 7) and putative stress response genes/proteins (n = 38) in the isolates.

3.5. Pathogenic Potential and Putative Virulence Factor Predictions

The genomes of the E. coli predicted a 0.937 mean probability (Pscore) of being a human host-pathogen and was found to match pathogenic families ranging between 525 and 1001 (Table 1). The whole-virulome analysis predicted differences in the virulence factors possessed by the isolates. Nineteen putative virulence genes of the five main classes (adhesins, toxins, immune evasion, iron uptake, and microcin) of E. coli (Table 1 and Figure 4). Of note, adhesins (n = 5 genes) and toxins (n = 5 genes) were the most predominant virulence factor class. The most common virulence factors were gad (n = 14), iss (n = 8), and lpfA (n = 8), with eilA, iha, cdtB, mchB, mchC, mchF, and mcmA occurring in single isolates (Table 1 and Figure 4). The isolates CC6 and SG6 harbored the highest number of genes while all the isolates contained at least one virulence gene. As shown in Figure 4, diverse virulence genes that were mostly clonal-specific were identified. A specimen source-virulome association comparison was made (Figure 4) and there was little evidence to suggest that strains from a specific source had more virulence genes than others.

3.6. Phylogenomic Analysis

The tree analyses via the whole genome sequencing single nucleotide polymorphism (WGS SNP) tree differentiated the strains into three clusters (A, B, and C) with each cluster containing sub-clusters (Figure 4). Cluster A had the largest number of isolates (n = 9). The E. coli isolates showed significant phylogenetic diversity (Figure 4) with the phylogenomics showing greater resolution than the MLST technique. For instance, the sub-cluster C2 contained isolates (CM4, SLC2, and TLC13) from the same clone (ST155), but phylogenetically found on different sub-branches. The tree analyses coupled with metadata showed no patterns among the isolates with respect to their sources and areas (Figure 4).

4. Discussion

E. coli are among the most common foodborne pathogens associated with infections reported from meat sources. However, information on the virulome, pathogenic potential, and adaptive characteristics associated with its persistence using high throughput technologies in Africa are lacking. Therefore, this study investigated the virulome, pathogenicity, stress response factors, clonal lineages, and the phylogenomic relationship of 14 E. coli isolated from different meat types in Ghana using WGS. Of the 14 isolates, eight (57.1%) were multi-drug resistance (defined as resistance to three or more different classes of antibiotics) and possessed multiple antibiotic resistome [41].
The GView Server was used for the comparative circular visualisation of the multiple genomes to explore the inter-relationship between the genomic data points [32]. The circular plots generated for genomic data is tailored to gain the best biological insights into the data [42]. The inspection of the generated circular maps revealed variations in the genomic synteny (90–97%) as well as distinct trends of majorly gapped regions from the reference (Figure 1), indicating a possible closer and diverse association between the study strains.
The diversity of sequence types identified in the E. coli isolates is consistent with findings reported from China, USA, Korea, and Ghana [43,44,45,46]. Whiles four sequence types, namely ST155, ST69, ST29, and ST1727, found in the study have already been reported in human, animal, and environmental samples globally [43,47,48]. Two sequence types known as ST6646(CB1) and ST7483(NB12) were not associated with any globally curated ancestry lineages (satellite-variant) (Figure 2). Interestingly, globally reported ST10 (n = 3/5) was the most common ancestry (single-locus variant) for the five remaining sequence types (ST44, ST469, ST540, ST1141, and ST7473) (Figure 2). The results of the FumC-FimH types corroborated with the MLST scheme and affirmed by the high genetic diversity between the isolates could be reliable to indicate the degree of clonality between E. coli strains.
The CRISPRCas defence system that controls the adaptive immunity mechanisms responsible for protecting bacteria against viral invasion was predicted in silico (Table 1) [49,50]. The CRISPRCas has been advanced as a critical tool for the survival of bacteria in the microbial environment. It is, thus, likely this system contributed to the survival of these strains in different settings [49,51]. Additionally, the Type I-E CRISPR/Cas cluster detected in the strains has been reported to be involved in transcription regulation of endogenous genes with a possibility to rewire the regulatory interactions with the selection of targets done through native adaptation [52,53]. Further analysis showed other defence systems such as the restriction-modification (R-M) types I, II, and IV offering protection against exogenous DNA and, thus, increasing their survival [54]. The detection of the 2 R-M systems in most strains was not unusual considering that about 12.5% of bacteria commonly harbour more than one R-M system [54]. Moreover, the presence of the Type II system (Table 1) in all the isolates affirmed its abundance in E. coli [54,55]. However, there is little information about its regulation and further studies should be done to understand its regulation.
Bacteria deal with continuous stress and, thus, develop mechanisms to tolerate different stressors and adapt to hostile micro-environments [9]. Several conserved cellular mechanisms aiding in the survival of these strains in extreme biological niches were identified (Figure 3 and Table S3). The E. coli isolates showed genomic features to control the pH level and osmolarity to adapt in the different environment. For instance, to withstand drastic exchanges in osmolarity, they contained molecules such as Aquaporin Z and glycerol uptake protein facilitator (glpF) found to play a vital role in water and glycerol movement, respectively [56,57] (Table S3). They also harboured the biosynthesis system (Osmoprotectant) and choline-glycine betaine uptake, which manage the feedback mechanism by regulating these molecules to control high salinity in their micro-environment. Moreover, they serve as sources of nitrogen and carbon for the isolates to grow [58]. In the case of oxidative stress, the strains contained enzymes (peroxidase enzyme, reductase, and super-oxidase dismutase) that can neutralise harmful products such as the reactive oxygen species (ROS) before damaging vital cellular components, such as proteins, DNA, and membrane lipids (Figure 3 and Table S3) [59]. The periplasmic tolerant response mechanisms reported to aid Enterobacteria in surviving acidic environments were found in the strains (Table S3) [60,61]. More so, the periplasmic membrane proteases (RasP/YluC) involved in the splitting of anti-sigma factor regulons (RpoE/RseA) to prevent protein misfolding caused in cell stress factors were found in the isolates [62]. The E. coli strains were armed against toxicants such as the harboured detoxification proteins (YbaT, CysA, DedA, and TehA/TehB) to fight the harmful effects on proteins caused by heavy metals (tellurite and selenite) in the microenvironment (Figure 3 and Table S3) [63,64]. However, biology bacterial stress response mechanisms are still part of on-going research. Therefore, more studies will be needed to ascertain the exact mechanisms controlling the phenotypes of the stress response to aid in the discovery of potential targets for development of new anti-bacterial treatments.
The pathogenic potential (Pscore) with the probability ranging from 0 to 1 is employed in the automatic estimation of the chance of an isolate being a human pathogen [37,65]. This theoretical prediction of the pathogenicity using machine-learning algorithms to distinguish among commensal or pathogenic strains estimated a higher mean probability (Pscore ≈ 0.937), implying the strains from non-human origin to be pathogenic to humans (Table 1). Nevertheless, given the complexity between host-pathogen interactions, other experiments (including cell culture assays) are required to find out the practical relation of this theoretical measure (pathogenic potential) as this machine learning algorithm can be overestimated and it should be interpreted with caution [66,67].
The diversity of the virulome shown in the isolates is concerning as they are involved in virulence and play a significant role in their survival and pathogenesis. The virulome patterns were clonal-specific. For instance, the iron acquisition genes (ireA, iroN) and microcin H47 system (mcmA, mchB, mchC and mchF genes) were only possessed by the ST469 (Figure 4). The iron-acquisition system enhances the successful persistence and pathogenesis of the strain in the host [68,69]. The presence of these iron-uptake genes is very worrying, as these genes have been reported in various pathogenic E. coli [70,71]. The microcin H47 system acts as a bactericidal antibiotic secreted under conditions of nutrient depletion and exerts potent antibacterial action against closely related species, offering a survival advantage for the ST469 strain [72,73]. The adhesin-encoding (f17 A, f17 G) genes known for causing intestinal and extraintestinal disease in animals and humans were found in ST469 and ST1727 [74]. More so, a 2014 study by Morgan et al. [75] reported the production of F17 fimbriae by pathogenic E. coli implicated in septicaemia and diarrhoea outbreaks in lambs and calves. Moreover, the clonal-specific putative toxins (senB, sat, cdtB, and celb) were also observed in the isolates (Figure 4). The toxins offer the strains the ability to mediate cell membrane damage, inducing the formation of cytokines, and killing or reducing neutrophils [76]. These toxins have been identified in some strains (Table 1 and Figure 4), and are associated with causing illness in animals and humans [77,78]. The similarity of the virulence determinants between clones from the diverse settings imply that they were mainly acquired vertically and not via horizontal gene transfer. Moreover, no association could be established between the specimen sources, areas, and the virulome.
The phylogenetic relationships of E. coli strains were analysed by WGS SNPs tree. The greater resolution power of WGS sequencing over MLST was shown by the clustering of similar STs on the different branches and in sub-clusters (Figure 4, sub-cluster C2), giving credence to the need to shift to genomic investigation for a better understanding of the evolution of pathogens. Furthermore, tree analyses coupled with metadata insights depicted the diversity among the E. coli isolates with no correlation with their meat sources and sampling areas. The use of visual analytics increases confidence during molecular investigation of pathogens [65,79,80].
The WGS data could not infer the different E. coli pathotypes owing to a multiplicity of factors, such as genes, under review following concerns about specificity and unknown essential virulence determinants for some pathotypes, which has been reported in literature [81,82]. Hence, not all isolates comply with the standard pathotyping scheme. Additionally, the upsurge of hybrid strains of E. coli pathotypes owing to the mobility of most of the genes that encode virulence makes it easy for gene sharing and causing problems of specificity [81,82]. Finally, the particular limitation of WGS pathotyping is its limited capacity to accommodate new strains that do not comply with known categories since the current nomenclature is unwieldy and inflexible. The current challenges should be resolved to harness the full capacity of WGS to aid our understanding of the pathobiology and evolution of this highly versatile and adaptable E. coli species from different settings.

5. Conclusions

The findings of this present study revealed the diverse clonal lineages in the E. coli isolates with a rich repertoire of defence systems, conserved battery of stress tolerant factors, highly predicted pathogenic potential, and clonal-specific putative virulence factors in various meat types within the Tamale metropolis of Ghana, which is worrying for meat contamination management. Phylogenomic analyses combined with a metadata confirmed this high genetic diversity with no correlation with their meat sources and areas. Further molecular surveillance studies involving a large data set are recommended to help food safety and public health practitioners in the design of contamination control strategies in meat retail settings in Ghana.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/12/1504/s1. Table S1. Prevalence of E. coli in meat samples sold at the Tamale Metropolis. Table S2. A table showing the eBURST (Based Upon Related Sequence Types) analyses of the study sequence types with global curated STs in Escherichia PubMLST database. Table S3. In silico identification and characterisation of conserved stress response mechanisms in the E. coli strains.

Author Contributions

Conceptualization, F.A. and D.G.A. Methodology, F.A., J.A., H.M.K., R.B.K., A.M.S., and D.G.A. Validation, D.G.A. and J.A. Formal analysis, D.G.A. and J.A. Investigation, F.A. and D.G.A. Resources, F.A., H.M.K., R.B.K., A.M.S., and D.G.A. Data curation, F.A., J.A., H.M.K., R.B.K., A.M.S., and D.G.A. Writing—original draft preparation, F.A. and D.G.A. Writing—review and editing, J.A., H.M.K., R.B.K., A.M.S., and D.G.A. Supervision, D.G.A. Project administration, F.A. and D.G.A. Funding acquisition, F.A., H.M.K., R.B.K., A.M.S., and D.G.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors are grateful to the University for Development Studies for making available a microbiology laboratory for this work. The sequencing was supported by the Commonwealth Science Conference (CSC) Follow-on Travel Grants (CSC\R1\170022) and Nanyang Technological University Research Initiative. We also acknowledged Moon Tay Yue Feng and Jorgen Schlundt of Nanyang Technological University Food Technology Centre (NAFTEC) for their assistance with the sequencing.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding this research.

References

  1. Octavia, S.; Lan, R. The Family Enterobacteriaceae. In The Prokaryotes; Springer: Berlin/Heidelberg, Germany, 2014; pp. 225–286. ISBN 9783642389221. [Google Scholar]
  2. Feng, P.; Weagant, S.D.; Grant, M.A.; Burkhardt, W. Enumeration of Escherichia coli and the Coliform Bacteria. In Bacteriological Analytical Manual; Food & Drug Administration (FDA): Silver Spring, MD, USA, 2018; pp. 21–34. [Google Scholar]
  3. Croxen, M.A.; Law, R.J.; Scholz, R.; Keeney, K.M.; Wlodarska, M.; Finlay, B.B. Recent Advances in Understanding Enteric Pathogenic Escherichia coli. Clin. Microbiol. Rev. 2013, 26, 822–880. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Luna-Guevara, J.J.; Arenas-Hernandez, M.M.P.; De La Peña, C.M.; Silva, J.L.; Luna-Guevara, M.L. The Role of Pathogenic E. coli in Fresh Vegetables: Behavior, Contamination Factors, and Preventive Measures. Int. J. Microbiol. 2019, 2019, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Nataro, J.P.; Kaper, J.B. Diarrheagenic Escherichia coli. Clin. Microbiol. Rev. 1998, 11, 142–201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Pitout, J.D.D. Extraintestinal Pathogenic Escherichia coli: A Combination of Virulence with Antibiotic Resistance. Front. Microbiol. 2012, 3, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Leimbach, A.; Hacker, J.; Dobrindt, U.E. coli as an All-Rounder: The Thin Line between Commensalism and Pathogenicity. In Current Topics in Microbiology and Immunology; Springer: Berlin/Heidelberg, Germany, 2013; pp. 3–32. ISBN 9783642365591. [Google Scholar]
  8. Bisi-Johnson, M.A.; Obi, C.L.; Vasaikar, S.D.; Baba, K.A.; Hattori, T. Molecular basis of virulence in clinical isolates of Escherichia coli and Salmonella species from a tertiary hospital in the Eastern Cape, South Africa. Gut Pathog. 2011, 3, 9. [Google Scholar] [CrossRef] [Green Version]
  9. Fang, F.C.; Frawley, E.R.; Tapscott, T.; Vázquez-Torres, A. Bacterial Stress Responses during Host Infection. Cell Host Microbe 2016, 20, 133–143. [Google Scholar] [CrossRef] [Green Version]
  10. Marles-Wright, J.; Lewis, R.J. Stress responses of bacteria. Curr. Opin. Struct. Biol. 2007, 17, 755–760. [Google Scholar] [CrossRef]
  11. Chung, H.J.; Bang, W.; Drake, M.A. Stress Response of Escherichia coli. Compr. Rev. Food Sci. Food Saf. 2006, 5, 52–64. [Google Scholar] [CrossRef]
  12. Van, T.T.H.; Chin, J.; Chapman, T.; Tran, L.T.; Coloe, P.J. Safety of raw meat and shellfish in Vietnam: An analysis of Escherichia coli isolations for antibiotic resistance and virulence genes. Int. J. Food Microbiol. 2008, 124, 217–223. [Google Scholar] [CrossRef]
  13. Lyhs, U.; Ikonen, I.; Pohjanvirta, T.; Raninen, K.; Perko-Mäkelä, P.; Pelkonen, S. Extraintestinal pathogenic Escherichia coli in poultry meat products on the Finnish retail market. Acta Vet. Scand. 2012, 54, 64. [Google Scholar] [CrossRef] [Green Version]
  14. Caruso, G.; Giammanco, A.; Cardamone, C.; Oliveri, G.; Mascarella, C.; Capra, G.; Fasciana, T. Extra-Intestinal Fluoroquinolone-Resistant Escherichia coli Strains Isolated from Meat. BioMed. Res. Int. 2018, 2018, 1–7. [Google Scholar] [CrossRef] [PubMed]
  15. Cross, A.S. What is a virulence factor? Crit. Care 2008, 12, 197. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Marangoni, F.; Corsello, G.; Cricelli, C.; Ferrara, N.; Ghiselli, A.; Lucchin, L.; Poli, A. Role of poultry meat in a balanced diet aimed at maintaining health and wellbeing: An Italian consensus document. Food Nutr. Res. 2015, 59, 27606. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Bohrer, B.M. Review: Nutrient density and nutritional value of meat products and non-meat foods high in protein. Trends Food Sci. Technol. 2017, 65, 103–112. [Google Scholar] [CrossRef]
  18. Tay, M.Y.F.; Adzitey, F.; Sultan, S.A.; Tati, J.M.; Seow, K.L.G.; Schlundt, J. Whole-Genome Sequencing of Nontyphoidal Salmonella enterica Isolates Obtained from Various Meat Types in Ghana. Microbiol. Resour. Announc. 2019, 8, e00033-19. [Google Scholar] [CrossRef] [Green Version]
  19. Rwego, I.B.; Gillespie, T.R.; Isabirye-Basuta, G.; Goldberg, T.L. High Rates of Escherichia coli Transmission between Livestock and Humans in Rural Uganda. J. Clin. Microbiol. 2008, 46, 3187–3191. [Google Scholar] [CrossRef] [Green Version]
  20. Currie, A.; Honish, L.; Cutler, J.; Locas, A.; Lavoie, M.-C.; Gaulin, C.; Galanis, E.; Tschetter, L.; Chui, L.; Taylor, M.; et al. Outbreak of Escherichia coli O157:H7 Infections Linked to Mechanically Tenderized Beef and the Largest Beef Recall in Canada, 2012. J. Food Prot. 2019, 82, 1532–1538. [Google Scholar] [CrossRef]
  21. Ostroff, S.M.; Griffin, P.M.; Tauxe, R.V.; Shipman, L.D.; Greene, K.D.; Wells, J.G.; Lewis, J.H.; Blake, P.A.; Kobayashi, J.M. A statewide outbreak of Escherichia coli 0157: h7 infections in washington state. Am. J. Epidemiol. 1990, 132, 239–247. [Google Scholar] [CrossRef]
  22. Brzuszkiewicz, E.; Thürmer, A.; Schuldes, J.; Leimbach, A.; Liesegang, H.; Meyer, F.-D.; Boelter, J.; Petersen, H.; Gottschalk, G.; Daniel, R. Genome sequence analyses of two isolates from the recent Escherichia coli outbreak in Germany reveal the emergence of a new pathotype: Entero-Aggregative-Haemorrhagic Escherichia coli (EAHEC). Arch. Microbiol. 2011, 193, 883–891. [Google Scholar] [CrossRef] [Green Version]
  23. Alonso, R.; Martín, A.; Peláez, T.; Marín, M.; Rodríguez-Creixéms, M.; Bouza, E. An improved protocol for pulsed-field gel electrophoresis typing of Clostridium difficile. J. Med. Microbiol. 2005, 54, 155–157. [Google Scholar] [CrossRef]
  24. Hollmén, T.E.; Debroy, C.; Flint, P.L.; Safine, D.E.; Schamber, J.L.; Riddle, A.E.; Trust, K.A. Molecular typing of Escherichia coli strains associated with threatened sea ducks and near-shore marine habitats of south-west Alaska. Environ. Microbiol. Rep. 2011, 3, 262–269. [Google Scholar] [CrossRef] [PubMed]
  25. Adzitey, F.; Ali, G.R.R.; Huda, N.; Ahmad, R. Genotyping of Salmonella strains isolated from ducks, their rearing and processing environments in Penang, Malaysia, using RAPD. 3 Biotech 2013, 3, 521–527. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Molechan, C.; Amoako, D.G.; Abia, A.L.K.; Somboro, A.M.; Bester, L.A.; Essack, S.Y. Molecular epidemiology of antibiotic-resistant Enterococcus spp. from the farm-to-fork continuum in intensive poultry production in KwaZulu-Natal, South Africa. Sci. Total Environ. 2019, 692, 868–878. [Google Scholar] [CrossRef] [PubMed]
  27. Pillay, S.; Amoako, D.G.; Abia, A.L.K.; Somboro, A.M.; Shobo, C.O.; Perrett, K.; Bester, L.A.; Essack, S.Y. Characterisation of Campylobacter spp. Isolated from Poultry in KwaZulu-Natal, South Africa. Antibiotics 2020, 9, 42. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Nutman, A.; Marchaim, D. How to: Molecular investigation of a hospital outbreak. Clin. Microbiol. Infect. 2019, 25, 688–695. [Google Scholar] [CrossRef] [PubMed]
  29. Lynch, T.; Petkau, A.; Knox, N.; Graham, M.; Van Domselaar, G. A Primer on Infectious Disease Bacterial Genomics. Clin. Microbiol. Rev. 2016, 29, 881–913. [Google Scholar] [CrossRef] [Green Version]
  30. Osei Sekyere, J.; Amoako, D.G. Genomic and phenotypic characterisation of fluoroquinolone resistance mechanisms in Enterobacteriaceae in Durban, South Africa. PLoS ONE 2017, 12, e0178888. [Google Scholar] [CrossRef] [Green Version]
  31. Neher, R.A.; Bedford, T. Real-Time Analysis and Visualization of Pathogen Sequence Data. J. Clin. Microbiol. 2018, 56, 1–15. [Google Scholar] [CrossRef] [Green Version]
  32. Petkau, A.; Stuart-Edwards, M.; Stothard, P.; Van Domselaar, G. Interactive microbial genome visualization with GView. Bioinformatics 2010, 26, 3125–3126. [Google Scholar] [CrossRef]
  33. Tatusova, T.; DiCuccio, M.; Badretdin, A.; Chetvernin, V.; Nawrocki, P.; Zaslavsky, L.; Lomsadze, A.; Pruitt, K.D.; Borodovsky, M.; Ostell, J. NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res. 2016, 44, 6614–6624. [Google Scholar] [CrossRef]
  34. Overbeek, R.; Olson, R.; Pusch, G.D.; Olsen, G.J.; Davis, J.J.; Disz, T.; Edwards, R.A.; Gerdes, S.; Parrello, B.; Shukla, M.; et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res. 2014, 42, D206–D214. [Google Scholar] [CrossRef] [PubMed]
  35. Feil, E.J.; Li, B.C.; Aanensen, D.M.; Hanage, W.P.; Spratt, B.G. eBURST: Inferring Patterns of Evolutionary Descent among Clusters of Related Bacterial Genotypes from Multilocus Sequence Typing Data. J. Bacteriol. 2004, 186, 1518–1530. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Couvin, D.; Bernheim, A.; Toffano-Nioche, C.; Touchon, M.; Michalik, J.; Néron, B.; Rocha, E.P.C.; Vergnaud, G.; Gautheret, D.; Pourcel, C. CRISPRCasFinder, an update of CRISRFinder, includes a portable version, enhanced performance and integrates search for Cas proteins. Nucleic Acids Res. 2018, 46, W246–W251. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Cosentino, S.; Voldby Larsen, M.; Møller Aarestrup, F.; Lund, O. PathogenFinder—Distinguishing Friend from Foe Using Bacterial Whole Genome Sequence Data. PLoS ONE 2013, 8, e77302. [Google Scholar] [CrossRef]
  38. Kleinheinz, K.A.; Joensen, K.G.; Larsen, M.V. Applying the ResFinder and VirulenceFinder web-services for easy identification of acquired antibiotic resistance and E. coli virulence genes in bacteriophage and prophage nucleotide sequences. Bacteriophage 2014, 4, e27943. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Ahrenfeldt, J.; Skaarup, C.; Hasman, H.; Pedersen, A.G.; Aarestrup, F.M.; Lund, O. Bacterial whole genome-based phylogeny: Construction of a new benchmarking dataset and assessment of some existing methods. BMC Genom. 2017, 18, 19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Hadfield, J.; Croucher, N.J.; Goater, R.J.; AbuDahab, K.; Aanensen, D.M.; Harris, S.R. Phandango: An interactive viewer for bacterial population genomics. Bioinformatics 2018, 34, 292–293. [Google Scholar] [CrossRef]
  41. Adzitey, F.; Assoah-Peprah, P.; Teye, G.A. Whole-genome sequencing of Escherichia coli isolated from contaminated meat samples collected from the Northern Region of Ghana reveals the presence of multiple antimicrobial resistance genes. J. Glob. Antimicrob. Resist. 2019, 18, 179–182. [Google Scholar] [CrossRef]
  42. O’Donoghue, S.I.; Baldi, B.F.; Clark, S.J.; Darling, A.E.; Hogan, J.M.; Kaur, S.; Maier-Hein, L.; McCarthy, D.J.; Moore, W.J.; Stenau, E.; et al. Visualization of Biomedical Data. Annu. Rev. Biomed. Data Sci. 2018, 1, 275–304. [Google Scholar] [CrossRef] [Green Version]
  43. Castellanos, L.R.; Donado-Godoy, P.; León, M.; Clavijo, V.; Arevalo, A.; Bernal, J.F.; Timmerman, A.J.; Mevius, D.J.; Wagenaar, J.A.; Hordijk, J. High Heterogeneity of Escherichia coli Sequence Types Harbouring ESBL/AmpC Genes on IncI1 Plasmids in the Colombian Poultry Chain. PLoS ONE 2017, 12, e0170777. [Google Scholar] [CrossRef]
  44. Yamaji, R.; Friedman, C.R.; Rubin, J.; Suh, J.; Thys, E.; McDermott, P.; Hung-Fan, M.; Riley, L.W. A Population-Based Surveillance Study of Shared Genotypes of Escherichia coli Isolates from Retail Meat and Suspected Cases of Urinary Tract Infections. mSphere 2018, 3, e00179-18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Tadesse, D.A.; Li, C.; Mukherjee, S.; Hsu, C.-H.; Bodeis Jones, S.; Gaines, S.A.; Kabera, C.; Loneragan, G.H.; Torrence, M.; Harhay, D.M.; et al. Whole-Genome Sequence Analysis of CTX-M Containing Escherichia coli Isolates from Retail Meats and Cattle in the United States. Microb. Drug Resist. 2018, 24, 939–948. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Rasmussen, M.M.; Opintan, J.A.; Frimodt-Møller, N.; Styrishave, B. Beta-Lactamase Producing Escherichia coli Isolates in Imported and Locally Produced Chicken Meat from Ghana. PLoS ONE 2015, 10, e0139706. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Doumith, M.; Day, M.; Ciesielczuk, H.; Hope, R.; Underwood, A.; Reynolds, R.; Wain, J.; Livermore, D.M.; Woodford, N. Rapid Identification of Major Escherichia coli Sequence Types Causing Urinary Tract and Bloodstream Infections. J. Clin. Microbiol. 2015, 53, 160–166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Ogura, Y.; Gotoh, Y.; Itoh, T.; Sato, M.P.; Seto, K.; Yoshino, S.; Isobe, J.; Etoh, Y.; Kurogi, M.; Kimata, K.; et al. Population structure of Escherichia coli O26: H11 with recent and repeated stx2 acquisition in multiple lineages. Microb. Genom. 2017, 3. [Google Scholar] [CrossRef] [PubMed]
  49. Bakr Shabbir, M.A.; Hao, H.; Shabbir, M.Z.; Hussain, H.I.; Iqbal, Z.; Ahmed, S.; Sattar, A.; Iqbal, M.; Li, J.; Yuan, Z. Survival and evolution of CRISPR-Cas system in prokaryotes and its applications. Front. Immunol. 2016, 7, 1–14. [Google Scholar] [CrossRef] [Green Version]
  50. Barrangou, R.; Marraffini, L.A. CRISPR-Cas Systems: Prokaryotes Upgrade to Adaptive Immunity. Mol. Cell 2014, 54, 234–244. [Google Scholar] [CrossRef] [Green Version]
  51. Klompe, S.E.; Sternberg, S.H. Harnessing “A Billion Years of Experimentation”: The Ongoing Exploration and Exploitation of CRISPR–Cas Immune Systems. Cris. J. 2018, 1, 141–158. [Google Scholar] [CrossRef] [Green Version]
  52. Bozic, B.; Repac, J.; Djordjevic, M. Endogenous Gene Regulation as a Predicted Main Function of Type I-E CRISPR/Cas System in E. coli. Molecules 2019, 24, 784. [Google Scholar] [CrossRef] [Green Version]
  53. Chang, Y.; Su, T.; Qi, Q.; Liang, Q. Easy regulation of metabolic flux in Escherichia coli using an endogenous type I-E CRISPR-Cas system. Microb. Cell Factories 2016, 15, 195. [Google Scholar] [CrossRef] [Green Version]
  54. Vasu, K.; Nagaraja, V. Diverse Functions of Restriction-Modification Systems in Addition to Cellular Defense. Microbiol. Mol. Biol. Rev. 2013, 77, 53–72. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Pleška, M.; Qian, L.; Okura, R.; Bergmiller, T.; Wakamoto, Y.; Kussell, E.; Guet, C.C. Bacterial Autoimmunity Due to a Restriction-Modification System. Curr. Biol. 2016, 26, 404–409. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Lu, D.; Grayson, P.; Schulten, K. Glycerol Conductance and Physical Asymmetry of the Escherichia coli Glycerol Facilitator GlpF. Biophys. J. 2003, 85, 2977–2987. [Google Scholar] [CrossRef] [Green Version]
  57. Borgnia, M.J.; Kozono, D.; Calamita, G.; Maloney, P.C.; Agre, P. Functional reconstitution and characterization of AqpZ, the E. coli water channel protein. J. Mol. Biol. 1999, 291, 1169–1179. [Google Scholar] [CrossRef]
  58. Figueroa-Soto, C.G.; Valenzuela-Soto, E.M. Glycine betaine rather than acting only as an osmolyte also plays a role as regulator in cellular metabolism. Biochimie 2018, 147, 89–97. [Google Scholar] [CrossRef]
  59. Ezraty, B.; Gennaris, A.; Barras, F.; Collet, J.-F. Oxidative stress, protein damage and repair in bacteria. Nat. Rev. Microbiol. 2017, 15, 385–396. [Google Scholar] [CrossRef]
  60. Bearson, S.; Bearson, B.; Foster, J.W. Acid stress responses in enterobacteria. FEMS Microbiol. Lett. 2006, 147, 173–180. [Google Scholar] [CrossRef]
  61. Guan, N.; Liu, L. Microbial response to acid stress: Mechanisms and applications. Appl. Microbiol. Biotechnol. 2020, 104, 51–65. [Google Scholar] [CrossRef] [Green Version]
  62. Kroos, L.; Akiyama, Y. Biochemical and structural insights into intramembrane metalloprotease mechanisms. Biochim. Biophys. Acta Biomembr. 2013, 1828, 2873–2885. [Google Scholar] [CrossRef] [Green Version]
  63. Turner, R.J.; Weiner, J.H.; Taylor, D.E. Selenium metabolism in Escherichia coli. BioMetals 1998, 11, 223–227. [Google Scholar] [CrossRef]
  64. Choudhury, H.G.; Beis, K. Tellurite-Resistance Protein TehA from Escherichia coli. In Encyclopedia of Metalloproteins; Springer: New York, NY, USA, 2013; pp. 2157–2160. [Google Scholar]
  65. Amoako, D.G.; Somboro, A.M.; Abia, A.L.K.; Allam, M.; Ismail, A.; Bester, L.A.; Essack, S.Y. Genome Mining and Comparative Pathogenomic Analysis of An Endemic Methicillin-Resistant Staphylococcus aureus (MRSA) Clone, ST612-CC8-t1257-SCCmec_IVd(2B), Isolated in South Africa. Pathogens 2019, 8, 166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Berg, G.; Martinez, J.L. Friends or foes: Can we make a distinction between beneficial and harmful strains of the Stenotrophomonas maltophilia complex? Front. Microbiol. 2015, 6, 1–6. [Google Scholar] [CrossRef] [PubMed]
  67. Deneke, C.; Rentzsch, R.; Renard, B.Y. PaPrBaG: A machine learning approach for the detection of novel pathogens from NGS data. Sci. Rep. 2017, 7, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Ratledge, C. Iron Metabolism and Infection. Food Nutr. Bull. 2007, 28, S515–S523. [Google Scholar] [CrossRef]
  69. Cassat, J.E.; Skaar, E.P. Iron in Infection and Immunity. Cell Host Microbe 2013, 13, 509–519. [Google Scholar] [CrossRef] [Green Version]
  70. Li, Y.; Dai, J.; Zhuge, X.; Wang, H.; Hu, L.; Ren, J.; Chen, L.; Li, D.; Tang, F. Iron-regulated gene ireA in avian pathogenic Escherichia coli participates in adhesion and stress-resistance. BMC Vet. Res. 2016, 12, 167. [Google Scholar] [CrossRef] [Green Version]
  71. Okeke, I.N.; Scaletsky, I.C.A.; Soars, E.H.; Macfarlane, L.R.; Torres, A.G. Molecular Epidemiology of the Iron Utilization Genes of Enteroaggregative Escherichia coli. J. Clin. Microbiol. 2004, 42, 36–44. [Google Scholar] [CrossRef] [Green Version]
  72. Duquesne, S.; Destoumieux-Garzón, D.; Peduzzi, J.; Rebuffat, S. Microcins, gene-encoded antibacterial peptides from Enterobacteria. Nat. Prod. Rep. 2007, 24, 708. [Google Scholar] [CrossRef]
  73. Rebuffat, S. Microcins in action: Amazing defence strategies of Enterobacteria. Biochem. Soc. Trans. 2012, 40, 1456–1462. [Google Scholar] [CrossRef]
  74. Cid, D.; Sanz, R.; Marín, I.; De Greve, H.; Ruiz-Santa-Quiteria, J.A.; Amils, R.; De La Fuente, R. Characterization of Nonenterotoxigenic Escherichia coli Strains Producing F17 Fimbriae Isolated from Diarrheic Lambs and Goat Kids. J. Clin. Microbiol. 1999, 37, 1370–1375. [Google Scholar] [CrossRef] [Green Version]
  75. Bihannic, M.; Ghanbarpour, R.; Auvray, F.; Cavalié, L.; Châtre, P.; Boury, M.; Brugère, H.; Madec, J.-Y.; Oswald, E. Identification and detection of three new F17 fimbrial variants in Escherichia coli strains isolated from cattle. Vet. Res. 2014, 45, 76. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Do Vale, A.; Cabanes, D.; Sousa, S. Bacterial Toxins as Pathogen Weapons Against Phagocytes. Front. Microbiol. 2016, 7, 42. [Google Scholar] [CrossRef]
  77. Krüger, A.; Lucchesi, P.M.A.; Sanso, A.M.; Etcheverría, A.I.; Bustamante, A.V.; Burgán, J.; Fernández, L.; Fernandez, D.; Leotta, G.; Friedrich, A.W.; et al. Genetic characterization of Shiga toxin-producing Escherichia coli O26:H11 strains isolated from animal, food, and clinical samples. Front. Cell. Infect. Microbiol. 2015, 5, 74. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Faïs, T.; Delmas, J.; Serres, A.; Bonnet, R.; Dalmasso, G. Impact of CDT Toxin on Human Diseases. Toxins 2016, 8, 220. [Google Scholar] [CrossRef] [PubMed]
  79. Amoako, D.G.; Somboro, A.M.; Abia, A.L.K.; Allam, M.; Ismail, A.; Bester, L.; Essack, S.Y. Genomic analysis of methicillin-resistant Staphylococcus aureus isolated from poultry and occupational farm workers in Umgungundlovu District, South Africa. Sci. Total Environ. 2019, 670, 704–716. [Google Scholar] [CrossRef]
  80. Ramsamy, Y.; Mlisana, K.P.; Amoako, D.G.; Allam, M.; Ismail, A.; Singh, R.; Abia, A.L.K.; Essack, S.Y. Pathogenomic Analysis of a Novel Extensively Drug-Resistant Citrobacter freundii Isolate Carrying a blaNDM-1 Carbapenemase in South Africa. Pathogens 2020, 9, 89. [Google Scholar] [CrossRef] [Green Version]
  81. Robins-Browne, R.M.; Holt, K.E.; Ingle, D.J.; Hocking, D.M.; Yang, J.; Tauschek, M. Are Escherichia coli Pathotypes Still Relevant in the Era of Whole-Genome Sequencing? Front. Cell. Infect. Microbiol. 2016, 6, 141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Tagini, F.; Greub, G. Bacterial genome sequencing in clinical microbiology: A pathogen-oriented review. Eur. J. Clin. Microbiol. Infect. Dis. 2017, 36, 2007–2020. [Google Scholar] [CrossRef]
Genes 11 01504 g001 550
Figure 1. Comparative visualisation of the isolates (n = 14) with the E. coli reference strain (ATCC25922, Accession number: CP009074). The map was constructed using the GView online server (https://server.gview.ca/). The concentric circles represent comparisons between ATCC25922 and, starting with the inner circle, genome assemblies from E. coli genomes (strain ID: AB1, AC1, CB1, CC6, CM4, NB12, NC3, NM11, SG6, SLC2, TG1, TG5, TLC1, and TLC13). Colour codes are given for each strain with a synteny identity, ranging from 90–97%.
Figure 1. Comparative visualisation of the isolates (n = 14) with the E. coli reference strain (ATCC25922, Accession number: CP009074). The map was constructed using the GView online server (https://server.gview.ca/). The concentric circles represent comparisons between ATCC25922 and, starting with the inner circle, genome assemblies from E. coli genomes (strain ID: AB1, AC1, CB1, CC6, CM4, NB12, NC3, NM11, SG6, SLC2, TG1, TG5, TLC1, and TLC13). Colour codes are given for each strain with a synteny identity, ranging from 90–97%.
Genes 11 01504 g001
Genes 11 01504 g002 550
Figure 2. A graphical depiction showing the eBURST (Based Upon Related Sequence Types) analyses of the study sequence types with global curated STs in Escherichia PubMLST database using the Achtman scheme. eBURST identified four similar clones (have been reported), five single-locus variants (SLV), and two satellite lineages (more distantly) with global curated E. coli STs from animal (food), human, and environmental sources. The colour codes depicted the type of variant (blue), study isolates (STs) (orange), and closet global ancestry (green).
Figure 2. A graphical depiction showing the eBURST (Based Upon Related Sequence Types) analyses of the study sequence types with global curated STs in Escherichia PubMLST database using the Achtman scheme. eBURST identified four similar clones (have been reported), five single-locus variants (SLV), and two satellite lineages (more distantly) with global curated E. coli STs from animal (food), human, and environmental sources. The colour codes depicted the type of variant (blue), study isolates (STs) (orange), and closet global ancestry (green).
Genes 11 01504 g002
Genes 11 01504 g003 550
Figure 3. A colour gradient generated phylogenetic tree of distribution of stress response factors across the E. coli isolates (n = 14). The green colour represents the presence of the gene, and the yellow colour represents the absence of the gene. Conserved stress response factors (osmotic stress, oxidative stress, detoxification, and periplasmic stress) were found in all the strains irrespective of their sources, area, and sequence types.
Figure 3. A colour gradient generated phylogenetic tree of distribution of stress response factors across the E. coli isolates (n = 14). The green colour represents the presence of the gene, and the yellow colour represents the absence of the gene. Conserved stress response factors (osmotic stress, oxidative stress, detoxification, and periplasmic stress) were found in all the strains irrespective of their sources, area, and sequence types.
Genes 11 01504 g003
Genes 11 01504 g004 550
Figure 4. A phylogenomic branch coupled with metadata of isolate information (source, area, sequence types (STs), and virulence factors) showing the relationship between the study isolates (n = 14) using Phandango. The colour codes differentiated the different sources, areas, and STs. The virulence factors are represented by Roman numerals I: immune evasion, II: adhesins, III: toxins, IV: iron uptake, V: microcins.
Figure 4. A phylogenomic branch coupled with metadata of isolate information (source, area, sequence types (STs), and virulence factors) showing the relationship between the study isolates (n = 14) using Phandango. The colour codes differentiated the different sources, areas, and STs. The virulence factors are represented by Roman numerals I: immune evasion, II: adhesins, III: toxins, IV: iron uptake, V: microcins.
Genes 11 01504 g004
Table
Table 1. Summary of the population, the source of specimens, sample type, and genotypic characteristics of the isolates.
Table 1. Summary of the population, the source of specimens, sample type, and genotypic characteristics of the isolates.
1 DDemographic InformationIn Silico TypingMolecular Defence SystemsPathogenic Determinants
SourceAreaMLSTSerotype 1CHTyperCRISPRs (Cas)
(Cluster)		Restriction-Modification (R-M) SystemPathogenicity Score
(Pathogenic Families)		Virulence Factors
FumCFimH
AB1BeefAboabo540O11fumC7fimH547 (1)
CAS-TypeIE		Type II (M.EcoGVI, M.Eco4255 Dcm)
Type IV (EcoAPECGmrSD)		0.923 (525)gad, iss
CB1BeefCentral
Market		6646H39fumC203fimH5623 (0)Type I (S.Sso30807 I)
Type II (M.EcoKII, M.EcoGI, M.Eco3609 Dcm)		0.938 (698)air, cdtB, eilA, gad, iss
NB12BeefNyorni7483O8fumC23fimH383 (1)
CAS-TypeIE		Type I (S.EcoNIH1 II, EcoO127 I, M.EcoNIH1 II)
Type II (M.EcoKII, M.Eco3740 Dcm)		0.942 (710)astA, gad, lpfA
AC1ChevonAboabo44O162fumC11fimH548 (1)
CAS-TypeIE		Type I (M.Eco6409 I, S.Eco3325 I, M.Eco2747 II)
Type II (M.EcoGVII, M.EcoKII, M.EcoJA03 PDcm)		0.938 (698)gad
CC6ChevonCentral
Market		469O8fumC65fimH686 (1)
CAS-TypeIE		Type I (M.Sso30807 I),
Type II (M.EcoKII, M.Eco3740 Dcm)		0.940 (678)celb, f17 A, f17 G, gad, iss, lpfA mchB, mchC, mchF, mcmA
NC3ChevonNyorni1727O88fumC19fimH313 (1)
CAS-TypeIE		Type II (M.EcoKII, M.EcoGVII, M.EcoJA03 PDcm)0.940 (686)f17 A, f17 G, gad, iss, lpfA
SG6Guinea
fowl		Victory
Cinema		69O15fumC35fimH276 (2)
CAS-TypeIE		Type I (EcoKI, M.EcoJA69 PI, S.EcoJA69 PI),
Type II (M.EcoGI, M.EcoKII, M.Eco3317 Dcm)
Type IV (EcoKMrr)		0.930 (1001)air, eilA, gad, iha, ireA, iroN, iss, IpfA, sat, senB
TG1Guinea
fowl		Tishegu540O9fumC7fimH548 (1)
CAS-TypeIE		Type II (M.EcoGVII, M.EcoJA03 PDcm)0.934 (564)gad, iss
TG5Guinea
fowl		Tishegu7473O61fumC7fimH5410 (0)Type II (M.EcoKII, M.EcoJA03 PDcm)0.929 (604)astA, gad
SLC2Local ChickenVictory
Cinema		155H9fumC4fimH324 (1)
CAS-TypeIE		Type II (M.EcoKII, M.EcoGVII, M.EcoJA03 PDcm)0.943 (651)gad, lpfA
TLC1Local ChickenTishegu297H9fumC65fimH388 (1)
CAS-TypeIE		Type I (M.Eco3609 I)
Type II (M.EcoKII, M.EcoGI, M.Eco3609 Dcm)		0.943 (691)gad, lpfA
TLC13Local ChickenTishegu155O132fumC4fimH324 (1)
CAS-TypeIE		Type II (M.EcoKII, M.EcoGVII, M.EcoJA03 PDcm)0.943 (648)gad, lpfA
CM4MuttonCentral
Market		155H40fumC4fimH3666 (1)
CAS-TypeIE		Type I (M.Eco2747 II, EcoR124 II),
Type II (M.EcoE1140 Dcm, M.EcoKII)		0.940 (655)gad, iss, lpfA
NM11MuttonNyorni1141O113fumC11fimH259 (1)
CAS-TypeIE		Type I (M.EcoJA65 PI),
Type II (M.EcoJA03 PDcm, M.EcoKII)		0.936 (623)gad, iss
1 Not all the O-antigen and H-antigen loci in the strains were predicted.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Adzitey, F.; Asante, J.; Kumalo, H.M.; Khan, R.B.; Somboro, A.M.; Amoako, D.G. Genomic Investigation into the Virulome, Pathogenicity, Stress Response Factors, Clonal Lineages, and Phylogenetic Relationship of Escherichia coli Strains Isolated from Meat Sources in Ghana. Genes 2020, 11, 1504. https://doi.org/10.3390/genes11121504

AMA Style

Adzitey F, Asante J, Kumalo HM, Khan RB, Somboro AM, Amoako DG. Genomic Investigation into the Virulome, Pathogenicity, Stress Response Factors, Clonal Lineages, and Phylogenetic Relationship of Escherichia coli Strains Isolated from Meat Sources in Ghana. Genes. 2020; 11(12):1504. https://doi.org/10.3390/genes11121504

Chicago/Turabian Style

Adzitey, Frederick, Jonathan Asante, Hezekiel M. Kumalo, Rene B. Khan, Anou M. Somboro, and Daniel G. Amoako. 2020. "Genomic Investigation into the Virulome, Pathogenicity, Stress Response Factors, Clonal Lineages, and Phylogenetic Relationship of Escherichia coli Strains Isolated from Meat Sources in Ghana" Genes 11, no. 12: 1504. https://doi.org/10.3390/genes11121504

APA Style

Adzitey, F., Asante, J., Kumalo, H. M., Khan, R. B., Somboro, A. M., & Amoako, D. G. (2020). Genomic Investigation into the Virulome, Pathogenicity, Stress Response Factors, Clonal Lineages, and Phylogenetic Relationship of Escherichia coli Strains Isolated from Meat Sources in Ghana. Genes, 11(12), 1504. https://doi.org/10.3390/genes11121504

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop