Molecular Typing Of Diarrheagenic Escherichia Coli Strains Isolated From Children Under Five Years Using ERIC-PCR in Koula-Moutou, Gabon

Abstract

Introduction

Diarrheagenic Escherichia coli (DEC) are enteric pathogens responsible for diarrhea in children. The aim of this study was to assess the genetic diversity of DEC in children under five years of age in Koula-Moutou, a topic that has not been previously investigated in Gabon.

Methods

A total of 41 DEC strains, previously isolated and characterized, were included in this study. Genetic diversity was assessed using the Enterobacterial Repetitive Intergenic Consensus-Polymerase Chain Reaction (ERIC-PCR) method.

Results

The ERIC-PCR method enabled typing of 34 out of the 41 DEC isolates. These 34 strains were differentiated into 27 distinct ERIC-PCR genetic profiles, with fingerprints consisting of 1 to 6 bands ranging in size from 250 to 1700 bp. The profiles were grouped into 10 clusters with 90% similarity, each cluster containing between 1 and 7 strains. The DEC isolates from children exhibited diverse genetic profiles, indicating a high degree of genotypic polymorphism.

Discussion

27 different and repeated genetic profiles revealed, suggesting a great diversity of circulating clones, which present different resistance profiles, and which may make the treatment ineffective or inappropriate, reflecting the genetic evolution of DEC in this region of Gabon.

Conclusion

ERIC-PCR is a simple and rapid method for assessing the genetic diversity of DEC. This study revealed significant genetic variation among DEC strains and demonstrated that ERIC-PCR is a valuable tool for epidemiological surveillance and health-related research.

Keywords: Diarrheagenic E. coli, ERIC-PCR, Genotyping, Children, Koula-Moutou, Gabon.

1. INTRODUCTION

Diarrhea is one of the most commonly occurring waterborne diseases in developing countries [1]. It mainly affects children and causes morbidity and mortality in infants [2]. This disease is caused by a wide range of viruses, parasites, and bacteria [3]. Among the bacterial pathogens involved, diarrheagenic Escherichia coli (DEC) is an important causative agent of endemic and epidemic diarrhea worldwide [2, 4]. There are many strains of DEC, categorized into six main intestinal pathotypes according to the generated clinical signs and expressed pathogenicity factors. The various pathotypes described are: Enterotoxigenic Escherichia coli (ETEC), Enteropathogenic Escherichia coli (EPEC), Enteroaggregative Escherichia coli (EAggEC), Enterohemorrhagic Escherichia coli (EHEC), Enteroinvasive Escherichia coli (EIEC), and diffuse adherent Escherichia coli (DAEC).

The prevalence of different enteric pathogens, including DEC, varies according to geographical locations [5]. In addition, several studies have suggested that DECs are responsible for about 30-40% of acute diarrheal episodes in children. These studies have shown that DECs were the most predominant pathotypes, playing an important role in diarrheal morbidity [6-9]. In Africa, and particularly in Nigeria, prevalence rates of 18.4% and 35% were recorded in Abuja and Ile-Ife, respectively [10, 11]. In Gabon, a 2021 study highlighted the phenotyping and molecular typing in the city of Koula-Moutou with DEC detection of 63.5% and 68.5%, respectively [12].

The identification by an accurate and rapid diagnosis can help prevent alarming clinical developments and epidemics [13, 14]. Many methods for strain-typing and origin-tracking of infections exist and are essential to help understand the epidemiology of infectious diseases [15]. Among these, the Enterobacterial Repetitive Intergenic Consensus-Polymerase Chain Reaction (ERIC-PCR) offers a clear and reliable discriminatory power, as demonstrated by numerous studies [15-17]. Among all methods that are used to study the genetic diversity of bacteria in community and hospital settings, ERIC-PCR stands out because of its proven strong discriminatory power, reproducibility, accuracy, specificity, speed, and ease of use, making it a useful tool for routine epidemiological investigations and surveillance in developing countries. Consequently, several studies have used ERIC-PCR as a genotyping method to identify the molecular diversity and genetic relatedness of bacterial strains, particularly E. coli [18]. Similar observations related to the effectiveness of ERIC-PCR were reported, and circulating βLSE genes among E. coli isolates were screened in another study [14]. Likewise, in a study conducted in India, 14 clusters of multi-resistant E. coli strains and βLSE producers were reported by Jena et al. [19]. Equally, Ateba and Mbewe's work on monitoring E. coli O157:H7 infections isolated from several origins revealed remarkable similarities (72% to 91%) between obtained ERIC profiles [20]. Thus, ERIC-PCR gives more sensitive typing results than conventional methods and can be used to discriminate among clones or different members of bacterial species [14]. Although most strains of E. coli are non-pathogenic, those that are considered pathogenic, like DEC, can be the causes of some clinical complications.

This work, which forms part of the study of the genetic relatedness of DEC isolated from a rural city in Gabon, is a pioneering study. Using ERIC-PCR, this study aimed to investigate the genetic diversity of DEC isolates from children under 5 years in Koula-Moutou, a middle-eastern city of Gabon.

2. MATERIALS AND METHODS

2.1. Sample Source and DEC Isolation

The biological material used was supplied and obtained from the collection of strains available at the Laboratory of Molecular and Cellular Biology (LABMC) of the University of Sciences and Technology of Masuku (USTM), Franceville, Gabon. This material was made up of 41 diarrheagenic Escherichia coli (DEC) strains, isolated from childhood diarrheal feces at the Paul Moukambi Regional Hospital Centre in Koula-Moutou, Gabon. The characterization of these DECs has been carried out and described in previous work [21].

2.2. DNA Extraction and Quantification

The extraction of the total DNA was carried out using the heating technique. Young bacterial colonies of 18 h were collected, washed, and homogenized in 300 µL of ultra-pure water. Then, the samples were incubated for 20 min at -20°C and transferred to a dry water bath at 100°C for 10 min. Thereafter, the samples were cooled down for 10 min at -20°C. A volume of 40 µL of protein K at 100 µg/mL was added, and the preparations were incubated at 56°C for 1 h. The resulting lysates were centrifuged at 15,000 rpm for 10 min; the supernatants were transferred to new tubes and kept at -20°C until use.

These supernatants were then quantified using the Invitrogen Qubit™ HS Assay (ThermoFisher, France) with the Qubit Fluorometer. Then, a volume of 1400 µL of the mixture, made of 1393 µL of Buffer Qubit dsDNA HS and 7 µL of Qubit reagents, was freshly prepared and homogenized using a vortex. Previously, 198 and 190 µL of the working solution were collected and allocated into the sample quantification tubes and the standard tubes, respectively. Two microliters (2 µL) of each sample were added to 198 µL of working solution, and 10 µL of standard (S1/S2) to 190 µL of working solution in the appropriate tubes. Then, each sample was shaken for about 2–3 s, and the tubes were incubated for 2 min at room temperature. The sample concentrations were determined on the basis of the calibration curve (S1/S2) and were valid for the volume of 2 µL.

2.3. ERIC-PCR Amplification

Genetic diversity was investigated using the ERIC-PCR method on the 41 DEC strains. The sequences of the different primers used were (5’ATGTAAGCTCCTGGGGA TTCAC3’) and (5’AAGTAAGTGACTGGGGTGAGCG3’) for ERIC-F and -R, respectively.

From total DNA diluted to 1/3, the amplification program was carried out as described in previous studies [22]. The PCR reactions were carried out using a thermocycler (Bio-RAD, T100™, USA) in a final reaction volume of 20 µL containing 2 µL of DNA, 10 µL of master mix at 2X (Golden 360 Taq Amp, ThermoFisher), 2 µL of sense and anti-sense primers (10 pmol), 1.2 µL MgCl2 to 1.5 nM, and 4.8 µL ultra-pure water. The PCR reaction conditions and amplification programs were as follows: initial denaturation at 95°C for 5 min, followed by 30 denaturation cycles at 94°C for 1 min; hybridization at 52°C for 1 min, elongation at 72°C for 2 min, and final elongation at 72°C for 5 min.

2.4. Migration and visualization on Agarose Gel

Amplicons were separated by electrophoresis on a 1.5% (w/v) agarose gel prepared in 0.5X Tris borate EDTA (TBE) solution, containing ethidium bromide (BET). Migration was performed at 200 V and 300 mA for 1 h using the Bio-RAD apparatus (Power Pac™ Basic, USA). In each agarose gel electrophoresis of the amplified products, DNA size markers of 100 bp and 1 kb (Invitrogen, DNA Ladder, ThermoFisher Scientific) were incorporated.

2.5. Data Analysis of Electrophoretic Profiles and Creation of Clusters

The gels were read by the Quantum software (ST4-1100/26MX, France). The manual zero-one method was used to analyze the models using the online data analysis service (http://insilico.ehu.es/dice_upgma/), which uses the Dice index and UPGMA match analysis. It allows comparison of the strains according to their similarities and groups of the most similar individuals, which is to say, with the smallest genetic distance. A dendrogram was then produced according to the obtained clusters.

3. RESULTS

The main objective of this study was to determine the distribution of ERIC sequences within strains of DEC isolated from childhood diarrhea. The analysis of the genetic diversity of the 41 DEC strains by ERIC-PCR revealed a total of 82.9% (34/41) positivity (Fig. 1). Out of the 41 tested isolates, 34 could be successfully typed.

3.1. ERIC-PCR Amplification and Profiling Identification

Fig. 1 presents the analysis of results, revealing DNA fingerprints or profiles of amplified DNA bands ranging from 1 to 7 bands with sizes between 250 and 1700 bp. Common bands were those with sizes of 1450, 1300, 100, 500, and 400 bp, while less frequent bands included those weighing ≥1500, 900, 600, 300, and 250 bp. Seven of the 41 DEC strains did not yield any band, so they could not be genotyped (Fig. 1).

Profile identification showed 27 genotypes with 21 non-repeated and 6 repeated profiles (Table 1). These 27 profiles represented 34 of the 41 DEC strains.

The distribution of profiles according to the pathovars and environment indicate the existence of different genetic profiles (Fig. 2).

Table 1.
Identification of profiles, bands and clusters in association with the DEC strains from Koula-Moutou, Gabon.
Stains Diversity of bands Minimum bands Maximum bands Dominant bands (%) Number of profiles Unrepeated profiles Repeated profiles Number of clusters
41 13 250 1700 1450 (73.2) 27 21 6 10
1000 (39.0)
Fig (1).

Electrophoretic gel presenting ERIC-PCR profiles of DEC from Koula-Moutou, Gabon. PM: 100 bp DNA ladder; T-: negative control (UltraPure™ DNase/RNase-Free Distilled Water) T+: positive control (E. coli35218); E2-E20: some amplified DEC isolates.

Fig (2).

UPMA dendrogram image obtained from cluster analysis showing the relationship of 41 DEC strains isolated from Koula-Moutou, Gabon, using the ERIC-PCR technique.

3.2. Profile and Cluster Analysis

Dendrogram analysis showed that ERIC-PCR differentiated between 34 DEC isolates and categorized them into 10 clusters (C1–C10) (Fig. 2). Results also showed a high degree of similarity among the isolates, with similarity indices ranging from 90% for most isolates to 100% for identical isolates. Depending on their specificity, clusters 1, 2, and 10 had the same type of pathovars (EHEC and ETEC) but in different proportions. Clusters 3, 4, and 7 consisted of the same pathovar type (EPEC); clusters 6 and 8 contained the pathovar EHEC, while cluster 5 contained three different pathovars (EAEC, EHEC, and ETEC).

Of the 34 positive ERIC-PCR strains and the distribution of pathovars by cluster (C), 7 were contained in C5, 6 in C1 and C10, 5 in C2, 4 in C8, 2 in C7 and 1 in C3, C4, C6 and C9, respectively (Fig. 2).

4. DISCUSSION

Several molecular typing studies have been conducted on strains of Escherichia coli, but fewer studies have been done on molecular typing of DEC strains by ERIC-PCR [13, 17, 22, 23]. Using the ERIC-PCR method, the purpose of this study was to determine the genetic diversity of DEC strains isolated in children in Koula-Moutou, Gabon. This genetic prospection of DEC was done on a collection of 41 previously screened and characterized DEC strains [12].

The typing of DEC strains in this study revealed 27 different and repeated genetic profiles, suggesting a wide diversity of circulating clones, which were similar to those reported in the study by Roy et al. [24]. These 27 different genetic profiles can have different resistance profiles and make the treatment ineffective or unsuitable. As a result, patient care can be complicated. In addition, these different clones can be spread through pathways such as low sanitation, unprotected water sources, and non-compliance with good healthcare practices.

Although the band sizes greatly varied from 250 to 1700 bp between strains, the 1450 bp band was present in 30 of the 34 (88.24%) strains. Previous work on the typing of E. coli isolated from avian strains recorded electrophoretic band sizes ranging from 232 to 2690 bp [13]. The difference in sizes would be the result of deletions or acquisitions of sequences at the ERIC sequence level [25].

The high level of DNA profiles of DECs (82.9%) could be explained by the discriminatory power of two ERIC primers (ERIC1R, ERIC2, and ERIC1R+ERIC2). Indeed, it has been previously shown that the ERIC1R and ERIC2 used in our study have a greater ability to discriminate within the E. coli populations [26, 27]. These results are similar to those on 63 enteroaggregative Escherichia coli isolates from Saudi Arabia and revealed 63 unique profiles after typing by ERIC-PCR [23]. Ranjbar et al. [17] identified 115 strains of E. coli from 120 strains isolated from animal feces using ERIC-PCR. Strong discriminatory powers have also been observed with uropathogens [13, 22]. All these results show the existence of considerable variations between different DEC strains, particularly within the same pathovar type, correlated by the presence of one or more elements in the inter-gene region. As a result, the EPEC strains identified had different ERIC profiles, which were grouped in different clusters and originated from different areas of the city of Koula-Moutou. All these data seem to converge on the hypothesis that these strains originally had a common ancestor and then evolved over time by retaining some genetic material. Furthermore, some genetic fragments appear to have been preserved since the divergence of DEC [25]. This conservation of ERIC sequences would undoubtedly be linked to the acquisition of a function certainly related to the stability of the mRNA or a functional role in the bacterial cell, as suggested by other authors [25, 28]. Three DEC pathovars (EHEC, ETEC, EAEC) had similar repeated ERIC profiles, which were grouped into clusters sharing the type of pathovar and the type of patient (hospitalized or community), suggesting the existence of an endemic clone predominant in the city of Koula-Moutou. These results also reflect the genetic evolution of DEC in this region of Gabon. Moreover, the discriminatory power, reproducibility, precision, specificity, speed, and ease of use of ERIC-PCR make it a useful tool for routine epidemiological surveys and surveillance in this country as well as elsewhere.

STUDY LIMITATIONS

Despite the fact that this study showed relevant results on the relatedness of DEC in community and hospital settings, the sample size (41 strains tested) could be a limitation to this study. In addition, although molecular typing methods are commonly used to investigate epidemiological relationships among isolates and sources of infection, it would be advisable, before being used for those purposes, PCR methods for molecular typing require careful in-house validation of typeability, reproducibility, repeatability, stability, discriminatory power and epidemiologic concordance, which were not considered in this study. Also, it would have been interesting to use the MLST, pulsed field gel electrophoresis (PFGE) method or whole-genome sequencing (WGS), which are the gold-standard techniques and seem to be more discriminatory than the ERICPCR used. The uniqueness of these profiles could demonstrate that horizontal transmission and clonal dissemination can occur from clones of different or the same origin and over a greater distance [28].

CONCLUSION

This study showed a genetic difference and relatedness of DEC strains, and the ERIC-PCR method can be a tool for epidemiological health surveillance in Gabon. The genetic diversity revealed through techniques like ERIC-PCR further emphasizes the complexity of challenges in managing DEC infections, as similar clones circulate across different environments. This is the first work reporting on the genetic diversity of DECs in Gabon and revealing different genetic profiles after analyzing the feces of children in the city of Koula-Moutou.

AUTHORS’ CONTRIBUTIONS

The authors confirm their contribution to the paper as follows: J.F.Y., R.M.M., and F.M.: Coordinated the study and drafted the manuscript; R.B.M. and O.Z.M.: Carried out the molecular analyses; J.F.Y., R.M., and A.S:. Revised the article after it was drafted. All authors read and approved the final manuscript.

LIST OF ABBREVIATIONS

DAEC = Diffuse adhesion Escherichia coli
DEC = Diarrheagenic Echerichia coli
EAggEC = Enteroaggregative Escherichia coli
EHEC = Enterohemorrhagic Escherichia coli
EIEC = Enteroinvasive Escherichia coli
EPEC = Enteropathogenic Escherichia coli
ERIC-PCR = Enterobacterial Repetitive Intergenic Consensus-Polymerase Chain Reaction
ETEC = Enterotoxinogenic Escherichia coli
LABMC = Laboratory of Molecular and Cellular Biology
USTM = University of Sciences and Technology of Masuku

ETHICS APPROVAL AND CONSENT TO PARTICIPATE

No applicable

HUMAN AND ANIMAL RIGHTS

No applicable

CONSENT FOR PUBLICATION

No applicable.

FUNDING

None.

CONFLICT OF INTEREST

The authors declare no conflict of interest financial or otherwise.

ACKNOWLEDGEMENTS

The authors are grateful to the bacteriology immunology team at LABMC, Gabon for providing the bacterial materials necessary for this study and their availability.

REFERENCES

1
Graf S, Pichard C, Thibault R. Conduite à tenir lors d’une diarrhée sous nutrition entérale. Nutr Clin Metab 2013; 27(1): 51-3.
2
Somda NS, Bonkoungou OJI, Zongo C, et al. Prevalence of escherichia coli virulence genes in patients with Diarrhoea in Ouagadougou, Burkina Faso. Afr J Clin Exp Microbiol 2017; 18(4): 179-85.
3
Bonkoungou IJO, Haukka K, Österblad M, et al. Bacterial and viral etiology of childhood diarrhea in Ouagadougou, Burkina Faso. BMC Pediatr 2013; 13(1): 36.
4
Jafari A. Escherichia coli : Une brève revue des pathotypes diarrhéiques et leur rôle dans les maladies diarrhéiques en Iran Résumé. Iran J Microbiol 2012; 4(3): 102-17.
5
Patzi-Vargas S, Zaidi MB, Perez-Martinez I, et al. Diarrheagenic Escherichia coli carrying supplementary virulence genes are an important cause of moderate to severe diarrhoeal disease in Mexico. PLoS Negl Trop Dis 2015; 9(3): 0003510.
6
Thakur N, Jain S, Changotra H, Shrivastava R, Kumar Y, Grover N. Molecular characterization of diarrheagenic Escherichia coli pathotypes: Association of virulent genes, serogroups, and antibiotic resistance among moderate-to-severe diarrhea patients. J Clin Lab Anal 2018; 32(5): e22388.
7
Natarajan M, Kumar D, Mandal J, Biswal N, Stephen S. A study of virulence and antimicrobial resistance pattern in diarrhoeagenic Escherichia coli isolated from diarrhoeal stool specimens from children and adults in a tertiary hospital, Puducherry, India. J Health Popul Nutr 2018; 37(1): 17.
8
Canizalez-Roman A, Flores-Villaseñor HM, Gonzalez-Nuñez E, Velazquez-Roman J, Vidal JE, Muro-Amador S. Surveillance of diarrheagenic Escherichia coli strains isolated from diarrhea cases from children, adults and elderly at northwest of Mexico. Front Microbiol 2016; 7: 1924.
9
Coelho CVL, Gomes TAT, Vieira MAM, et al. Brazilian Indigenous children as carriers of diarrheagenic Escherichia coli pathotypes. Adv Infect Dis 2018; 8(4): 200-16.
10
Onanuga A, Igbeneghu O, Lamikanra A. A study of the prevalence of diarrhoeagenic Escherichia coli in children from Gwagwalada, Federal Capital Territory, Nigeria. Pan Afr Med J 2014; 17: 146.
11
Odetoyin BW, Hofmann J, Aboderin AO, Okeke IN. Diarrhoeagenic Escherichia coli in mother-child Pairs in Ile-Ife, South Western Nigeria. BMC Infect Dis 2015; 16(1): 28.
12
Mabika Mabika R, Mounioko F, Souza A, Yala JF. Phenotypic and molecular characterization of quinolone resistance of enteropathogens isolated from diarrhea in young children in Koula-Moutou/Gabon. Egypt J Basic Appl Sci 2021; 8(1): 345-53.
13
Movahedi M, Zarei O, Hazhirkamal M, Karami P, Shokoohizadeh L, Taheri M. Molecular typing of Escherichia coli strains isolated from urinary tract infection by ERIC-PCR. Gene Rep 2021; 23: 101058.
14
Tanko N, Tolulope OA, Olajumoke BR, Ong EBB, Yahaya M, Olalekan OB. Genetic relatedness in extended-spectrum beta-lactamase-producing Escherichia coli from clinical isolates using enterobacterial repetitive intergenic consensus polymerase chain reaction. Int J Health Sci 2021; 15(5): 18-27.
15
Simar SR, Hanson BM, Arias CA. Techniques in bacterial strain typing: Past, present, and future. Curr Opin Infect Dis 2021; 34(4): 339-45.
16
Behzadi P, Ranjbar R, Alavian SM. Nucleic Acid-based approaches for detection of viral hepatitis. Jundishapur J Microbiol 2014; 8(1): 17449.
17
Ranjbar R, Tabatabaee A, Behzadi P, Kheiri R. Enterobacterial repetitive intergenic consensus polymerase chain reaction (ERIC-PCR) genotyping of Escherichia coli strains isolated from different animal stool specimens. Iran J Pathol 2017; 12(1): 25-34.
18
Mokhtari S, Tahamtan Y, Kargar M, Tadayon K, Moazamian E. Genotyping Escherichia coli O157:H7 strains isolated from cattle and small ruminants in Fars province (Iran) using Enterobacterial Repetitive Intergenic Consensus PCR Analysis. Rev Élev Méd Vét Pays Trop 2024; 77(37397): 1-6.
19
Jena J, Debata NK, Sahoo RK, Gaur M, Subudhi E. Genetic diversity study of various β-lactamase-producing multidrug-resistant Escherichia coli isolates from a tertiary care hospital using ERIC-PCR. Indian J Med Res 2017; 146(Suppl. 1): S23-9.
20
Ateba C, Mbewe M. Genotypic characterization of Escherichia coli O157:H7 isolates from different sources in the North-West Province, South Africa, using enterobacterial repetitive intergenic consensus PCR analysis. Int J Mol Sci 2014; 15(6): 9735-47.
21
Mabika RM, Liabagui SLO, Moundounga HK, Mounioko F, Souza A, Yala JF. Molecular Prevalence and Epidemiological Characteristics of Diarrheagenic E. coli in Children under 5 Years Old in the City of Koula-Moutou, East-Central Gabon. Open J Med Microbiol 2021; 11(3): 157-75.
22
Afkhami Ardakani M, Ranjbar R. Molecular typing of uropathogenic E. coli strains by the ERIC-PCR method. Electron Physician 2016; 8(4): 2291-5.
23
Elhadi N, Aljindan R, Alsamman K, Alomar A, Aljeldah M. Antibiotic resistance and molecular characterization of enteroaggregative Escherichia coli isolated from patients with diarrhea in the Eastern Province of Saudi Arabia. Heliyon 2020; 6(4): 03721.
24
Roy S, Bhattacharjee B, Mazumder PB, Bhattacharjee M, Dhar D, Bhattacharjee A. Molecular characterization of enteropathogenic Escherichia coli isolated from patients with gastroenteritis in a tertiary referral hospital of northeast India. Indian J Med Microbiol 2024; 47: 100535.
25
Wilson LA, Sharp PM. Enterobacterial repetitive intergenic consensus (ERIC) sequences in Escherichia coli: Evolution and implications for ERIC-PCR. Mol Biol Evol 2006; 23(6): 1156-68.
26
Alsultan A, Elhadi N. Evaluation of ERIC-PCR method for determining genetic diversity among Escherichia coli isolated from human and retail imported frozen shrimp and beef. Int J Food Contam 2022; 9(1): 12.
27
Zong Minko O, Yala JF, Mabika Mabika R, Franck M, Ondjiangui LF, Armesto Paula J-M. Molecular typing of Campylobacter jejuni diarrhoeal isolates in a hospital in Makokou, Gabon by enterobacterial repetitive intergenic consensus PCR. Access Microbiol 2025; 7(3): 000947.v3.
28
Secundo de Souza AI, Freitas Neto OC, Batista DFA, et al. ERIC-PCR genotyping of field isolates of Salmonella enterica subsp. enterica serovar Gallinarum biovars Gallinarum and Pullorum. Avian Pathol 2015; 44(6): 475-9.