# Occurrence of Extended Spectrum Cephalosporin-, Carbapenem- and Colistin-Resistant Gram-Negative Bacteria in Fresh Vegetables, an Increasing Human Health Concern in Algeria

Widad Chelaghma, Lotfi Loucif, Esma Bendjama, Zineb Cherak, Mourad Bendahou, Jean-Marc Rolain

► To cite this version:

Widad Chelaghma, Lotfi Loucif, Esma Bendjama, Zineb Cherak, Mourad Bendahou, et al. Occurrence of Extended Spectrum Cephalosporin-, Carbapenem- and Colistin-Resistant Gram-Negative Bacteria in Fresh Vegetables, an Increasing Human Health Concern in Algeria. Antibiotics, 2022, 11 (8), pp.988. 10.3390/antibiotics11080988. [https://doi.org/10.3390/antibiotics11080988](https://doi.org/10.3390/antibiotics11080988)

HAL Id: hal-03941608

https://amu.hal.science/hal-03941608

Submitted on 16 Jan 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Occurrence of Extended Spectrum Cephalosporin-, Carbapenem- and Colistin-Resistant Gram-Negative Bacteria in Fresh Vegetables, an Increasing Human Health Concern in Algeria

Widad Chelaghma<sup>1,†</sup>, Lotfi Loucif<sup>2,*</sup>, Esma Bendjama<sup>2,3</sup>, Zineb Cherak<sup>4,†</sup>, Mourad Bendahou<sup>1</sup>
and Jean-Marc Rolain<sup>5,6</sup>

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Chelaghma, W.; Loucif, L.; Bendjama, E.; Cherak, Z.; Bendahou, M.; Rolain, J.-M. Occurrence of Extended Spectrum Cephalosporin-, Carbapenem- and Colistin-Resistant Gram-Negative Bacteria in Fresh Vegetables, an Increasing Human Health Concern in Algeria.

*Antibiotics* **2022**, *11*, 988. https://doi.org/10.3390/antibiotics11080988

Academic Editor: Corinna Kehrenberg

Received: 31 May 2022

Accepted: 4 July 2022

Published: 22 July 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** The aim of this study was to screen for extended spectrum cephalosporin-, carbapenem- and colistin-resistant Gram-negative bacteria in fresh vegetables in Batna, Algeria. A total of 400 samples of fresh vegetables were collected from different retail stores. Samples were immediately subjected to selective isolation, then the representative colonies were identified using matrix-assisted laser desorption ionisation time-of-flight mass spectrometry (MALDI-TOF-MS). Phenotypic and genotypic analyses were carried out in terms of species identification and relative antibiotic resistance. Transferability of the carbapenemase and *mcr*-bearing plasmids was verified by conjugation. The clonal relationships of carbapenemase and *mcr*-positive *Escherichia coli* isolates were studied by multi-locus sequence typing (MLST). Sixty-seven isolates were characterised and were mostly isolated from green leafy vegetables, where the dominant species identified included *Citrobacter freundii*, *Klebsiella pneumoniae*, *Enterobacter cloacae*, *Stenotrophomonas maltophilia*, *E. coli* and *Citrobacter braakii*. PCR and sequencing results showed that *E. coli* was the bacterial species presenting the highest antibiotic resistance level in parallel to *bla*<sub>TEM</sub> (*n* = 16) and *bla*<sub>CTX-M-15</sub> (*n* = 11), which were the most detected genes. Moreover, five isolates carried carbapenemase genes, including the *bla*<sub>OXA-48</sub> and/or *bla*<sub>VIM-4</sub> genes. The *mcr-1* gene was detected in two *E. coli* isolates. MLST analysis revealed three different *E. coli* sequence types: ST101 (*n* = 1), ST216 (*n* = 1) and ST2298 (*n* = 1). Conjugation assays confirmed the transferability of the *bla*<sub>OXA-48</sub> and *mcr-1* genes. In this study we report, for the first time, the detection of the *bla*<sub>OXA-48</sub> gene in *E. coli* and *C. braakii* isolates and the *bla*<sub>VIM-4</sub> gene in vegetables. To the best of our knowledge, this is the first report on the detection of *mcr-1* genes from vegetables in Algeria.

**Keywords:** fresh produce; ESBL; colistin resistance; carbapenem resistance; Batna

1. Introduction

Vegetables are considered to be a good source of essential nutrients for human health, particularly vitamins, dietary fibre, minerals and phytonutrients. In the last decade, the

<footnote n="1">¹ Laboratoire de Microbiologie Appliquée à l'Agroalimentaire, au Biomédical et à l'Environnement, Département de Biologie, Faculté des Sciences de la Nature et de la Vie et des Sciences de la Terre et de l'Univers, Université Abou Bekr Belkaid, Tlemcen 13000, Algeria; widadc2014@hotmail.com (W.C.); bendahou63@yahoo.fr (M.B.)</footnote>

<footnote n="2">² Laboratoire de Biotechnologie des Molécules Bioactives et de la Physiopathologie Cellulaire (LBMBPC), Faculté des Sciences de la Nature et de la Vie, Université Batna 2, Batna 05000, Algeria; bendjamaesma@hotmail.fr</footnote>

<footnote n="3">³ Département de Technologie Alimentaire, Institut des Sciences Vétérinaires et des Sciences Agronomiques, Université Batna 1, Batna 05000, Algeria</footnote>

<footnote n="4">⁴ Faculté des Sciences de la Nature et de la Vie, Université Batna 2, Batna 05000, Algeria;
cherak-zineb@hotmail.com</footnote>

<footnote n="5">5 Institut de Recherche pour le Développement (IRD), Microbes, Evolution, Phylogénie et Infection (MEPHI), Faculté de Médecine et de Pharmacie, Aix Marseille Université, 13005 Marseille, France; jean-marc.rolain@univ-amu.fr</footnote>

<footnote n="6">⁶ IHU Méditerranée Infection, Assistance Publique des Hôpitaux de Marseille, 13000 Marseille, France</footnote>

<footnote>* Correspondence: lotfiloucif@hotmail.fr; Tel.: +213-(0)-540-92-5400</footnote>
consumption of fresh produce, including vegetables, has increased, as part of healthier diets due to high nutrient density, correlated with low energy density [1,2]. In 2003, the World Health Organization (WHO) and the Food and Agriculture Organization (FAO) launched an initiative to support the consumption of fruit and vegetables to promote health around the world, with a recommended minimum intake of 400 g per day [2]. Because vegetables are often eaten raw, their consumption may result in the ingestion of bacteria able to pose a serious threat to consumer health [1]. Recently, fresh vegetables have been considered as potential vehicles of foodborne pathogens with different drug resistance levels and various, or most often unknown, sources of contamination [2,3]. Thus, vegetables represent a possible route of human exposure to antibiotic-resistant bacteria [1,4]. The antibiotic resistance threat has focused attention on all uses of antibiotics, including plant agriculture, where antibiotics play a primordial role as an excellent chemical tool for the control of various bacterial diseases [5]. In addition, the soil may well represent an important source of antibiotic-resistant bacteria for the plant due to several reasons such as the use of natural manure on crop fields [6]. Vegetables may also be exposed to putative sources of contamination, including: irrigation water as well as poor or inappropriate personal hygienic practices linked to post-harvest periods such as transport, market environment, methods of storage, processing and preparation [2,7].

One of the major concerns regarding antibiotic resistance worldwide is the dissemi-
nation of Gram-negative bacteria (GNB) which display resistance to β-lactam antibiotics,
particularly through the production of β-lactamases, including extended spectrum β-
lactamases (ESBLs), cephalosporinases (AmpC) and carbapenemases [1,8,9]. Carbapene-
mase production is one of the main mechanisms of resistance to carbapenems, which are a
powerful group of antimicrobial agents widely regarded by clinicians as the last therapeutic
line, particularly for the treatment of multidrug-resistant (MDR) bacterial infections [8].
The most frequent carbapenemase enzymes are KPC, NDM, VIM and OXA types, espe-
cially the OXA-48 variant. All of these carbapenemases have been reported in clinical and
environmental isolates of Gram-negative bacteria [9,10]. With the rising use of carbapenem
over the past two decades, a rapid worldwide increase has been found in the populations
of carbapenem-resistant GNB, making colistin the last therapeutic line for the treatment
of infections caused by such organisms. Nevertheless, the efficacy of this drug has been
challenged by the recent appearance and dissemination of mobilised colistin-resistant (mcr)
genes [8,11].

In this regard, the emergence of carbapenem- and colistin-resistant GNB is of great con-
cern to public health. These bacteria have been reported from different sources worldwide,
including humans, aquatic environments, animals and food products [12–14], but little is
still known about the carriage of carbapenem and colistin resistance genes in Gram-negative
bacteria on fresh vegetables [15]. The few reports focused on the detection of such levels of
drug resistance in GNB isolates on vegetables indicated that the majority of these bacteria
consist of environmental species, whereas faecal *Enterobacterales* species were occasionally
reported [4,16]. The ingestion of antibiotic-resistant bacteria is a potential public health
threat since they are able to colonise the gut or pass through the intestines and exchange
resistance genes with intestinal bacteria, facilitating their widespread dissemination in the
environment [17].

In Algeria, two studies have been published on antimicrobial resistance in veg-
etable samples in Bejaia city. The published studies showed the presence of OXA-48 and
ESBL (CTX-M-15, SHV-101, SHV-28 and OXA-1)-producing *K. pneumoniae* isolates [18,19].
Thus, the aim of the present study was to screen for extended spectrum cephalosporin-
carbapenem- and colistin-resistant GNB isolates on fresh vegetables in the city of Batna in
Algeria as well as to investigate their resistance mechanisms.
## 2. Results

### 2.1. Isolation and Identification of Presumptive Antibiotic-Resistant GNB

A total of sixty-seven GNB isolates were obtained; twenty isolates were originally from coriander (29.85%), fifteen were from carrot (22.39%) and parsley (22.39%), respectively, seven were from zucchini (10.45%), seven were from lettuce (10.45%), two were from turnip (2.98%) and one was from cucumber (1.49%). In this context, the Pearson chi-square test revealed a statistically significant effect of vegetable type (source) on the rate of positive isolated strains (positivity) ($\chi^2 = 118,004; p < 0.0001$).

In addition, season months have also a significant effect on the isolation rate, where 21 strains were isolated from 29 samples collected in summer months, while, in spring and fall months, 29 and 17 isolates were obtained from 216 and 112 vegetable samples, respectively. However, no isolates were obtained from samples (n = 43) collected in winter months.

The obtained isolates identified with a high proportion were *Citrobacter freundii* (20.90%; n = 14), *Klebsiella pneumoniae* (11.94%; n = 8), *Enterobacter cloacae* (11.94%; n = 8), *Stenotrophomonas maltophilia* (10.45%; n = 7), *Escherichia coli* (8.96%; n = 6) and *Citrobacter braakii* (8.96%; n = 6). The other identified species were *Shewanella putrefaciens* (2.99%; n = 2), *Aeromonas hydrophila* (2.99%; n = 2), *Serratia marcescens* (2.99%; n = 2), *Aeromonas caviae* (1.49%; n = 1), *Aeromonas eucrenophila* (1.49%; n = 1), *Acinetobacter pittii* (1.49%; n = 1), *Pseudomonas aeruginosa* (1.49%; n = 1), *Pseudomonas rhodesiae* (1.49%; n = 1), *Pseudomonas alcaligenes* (1.49%; n = 1), *Raoultella ornithinolytica* (1.49%; n = 1), *Proteus mirabilis* (1.49%; n = 1), *Serratia rubidæa* (1.49%; n = 1), *Serratia odorifera* (1.49%; n = 1), *Providencia rettgeri* (1.49%; n = 1) and *Rhizobium radiobacter* (1.49%; n = 1). Species identification results of the obtained GNB isolates are presented in Table 1, Figures 1 and 2.

<table><caption>**Table 1.** Bacterial species abundance and antibiotic resistance percentages.</caption><thead><tr><th>Bacterial Categories</th><th>Species Identification</th><th>Species Abundance</th><th>Antibiotics</th><th>Antibiotic Resistance Percentage</th></tr></thead><tbody><tr><td rowspan="14">Fermenting isolates</td><td>Citrobacter freundii</td><td>20.90%</td><td>Amoxicillin</td><td>96.23%</td></tr><tr><td>Klebsiella pneumoniae</td><td>11.94%</td><td>Cefoxitin</td><td>60.38%</td></tr><tr><td>Enterobacter cloacae</td><td>11.94%</td><td>Cefotaxime</td><td>75.47%</td></tr><tr><td>Escherichia coli</td><td>8.96%</td><td>Ceftazidime</td><td>62.26%</td></tr><tr><td>Citrobacter braakii</td><td>8.96%</td><td>Cefepime</td><td>39.62%</td></tr><tr><td>Aeromonas hydrophila</td><td>2.99%</td><td>Aztreonam</td><td>50.94%</td></tr><tr><td>Serratia marcescens</td><td>2.99%</td><td>Amoxicillin/clavulanate</td><td>79.25%</td></tr><tr><td>Aeromonas caviae</td><td>1.49%</td><td>Ertapenem</td><td>20.75%</td></tr><tr><td>Aeromonas eucrenophila</td><td>1.49%</td><td>Imipenem</td><td>1.89%</td></tr><tr><td>Raoultella ornithinolytica</td><td>1.49%</td><td>Tobramycin</td><td>20.75%</td></tr><tr><td>Proteus mirabilis</td><td>1.49%</td><td>Gentamicin</td><td>16.98%</td></tr><tr><td>Serratia rubidæa</td><td>1.49%</td><td>Amikacin</td><td>3.77%</td></tr><tr><td>Serratia odorifera</td><td>1.49%</td><td>Ciprofloxacin</td><td>26.42%</td></tr><tr><td>Providencia rettgeri</td><td>1.49%</td><td>Colistin</td><td>37.73%</td></tr><tr><td rowspan="10">Non-fermenting isolates</td><td>Stenotrophomonasmaltophilia</td><td>10.45%</td><td>Ticarcillin</td><td>92.86%</td></tr><tr><td>Shewanella putrefaciens</td><td>2.99%</td><td>Ceftazidime</td><td>78.57%</td></tr><tr><td>Acinetobacter pittii</td><td>1.49%</td><td>Cefepime</td><td>50%</td></tr><tr><td>Pseudomonas aeruginosa</td><td>1.49%</td><td>Aztreonam</td><td>85.71%</td></tr><tr><td>Pseudomonas rhodesiae</td><td>1.49%</td><td>Ticarcillin/clavulanate</td><td>92.86%</td></tr><tr><td>Pseudomonas alcaligenes</td><td>1.49%</td><td>Imipenem</td><td>78.57%</td></tr><tr><td></td><td></td><td>Tobramycin</td><td>21.43%</td></tr><tr><td></td><td></td><td>Gentamicin</td><td>7.14%</td></tr><tr><td></td><td></td><td>Amikacin</td><td>7.14%</td></tr><tr><td></td><td></td><td>Ciprofloxacin</td><td>0%</td></tr><tr><td></td><td>Rhizobium radiobacter</td><td>1.49%</td><td>Colistin</td><td>0%</td></tr></tbody></table>
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Antibiogram results, β-lactams and colistin resistance mechanisms and origins of our fermenting isolates clustered using the MultiExperimentViewer (MEV) software version 4_6_2. AX, amoxicillin; FOX, cefoxitin; CTX, cefotaxime; CAZ, ceftazidime; FEP, cefepime; ATM, aztreonam; AMC, amoxicillin/clavulanate; ETP, ertapenem; IPM, imipenem; TOB, tobramycin; GN, gentamicin; AK, amikacin; CIP, ciprofloxacin; R, resistant; S, susceptible; I, intermediate; Medium: the selective antibiotic added to the isolation culture medium. K: city of Kechida; P: city of Parc à fourrage; B: city of Bouakal; CH: city of Chouhada; D: downtown; M: 1020 household residence; H: city of Hamla; 1: 4/3/2019; 2: 18/3/2019; 3: 8/4/2019; 4: 15/4/2019; 5: 29/4/2019; 6: 9/6/2019; 7: 3/11/2019; 8: 17/11/2019.</figcaption></figure>

## 2.2. Antimicrobial Susceptibility Testing

The obtained fermenting isolates exhibited resistance to amoxicillin (96.23%; 51/53), cefoxitin (60.38%; 32/53), cefotaxime (75.47%; 40/53), ceftazidime (62.26%; 33/53), cefepime (39.62%; 21/53), aztreonam (50.94%; 27/53) and amoxicillin-clavulanic acid (79.25%; 42/53). Carbapenem resistance was detected on the basis of reduced susceptibility to ertapenem (20.75%; 11/53) or imipenem (1.89%; 1/53). Moreover, some of the isolates were detected to be resistant to aminoglycoside- and fluoroquinolone-tested antibiotics, including tobramycin (20.75%; 11/53), gentamicin (16.98%; 9/53), amikacin (3.77%; 2/53) and ciprofloxacin (26.42%; 14/53). Twenty isolates were resistant to colistin (37.73%; 20/53), with six isolates that presented an intrinsic resistance phenotype.

The fourteen remaining isolates identified as non-fermenting GNB were resistant to ticarcillin and ticarcillin-clavulanic acid, respectively (92.86%; 13/14), followed by aztreonam (85.71%; 12/14), ceftazidime (78.57%; 11/14), imipenem (78.57%; 11/14) (seven out of them presented a wild resistance phenotype to imipenem) and cefepime (50%; 7/14). The obtained isolates also presented a resistance to tobramycin (21.43%; 3/14), gentamicin (7.14%; 1/14) and amikacin (7.14%; 1/14). Ciprofloxacin and colistin were found to be the most active antibiotics, with 100% efficiency. The antibiogram results of the obtained GNB isolates clustered using the MultiExperimentViewer (MEV) software version 4_6_2 are presented in Figures 1 and 2 [20]. Antibiotic resistance percentages of the obtained species are presented in Table 1.
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Antibiogram results, β-lactam resistance mechanisms and origins of our non-fermenting isolates clustered using the MultiExperimentViewer (MEV) software version 4_6_2. TIC, ticarcillin; CAZ, ceftazidime; FEP, cefepime; ATM, aztreonam; TIM, ticarcillin/clavulanate; IPM, imipenem; TOB, tobramycin; GN, gentamicin; AK, amikacin; CIP, ciprofloxacin; R, resistant; S, susceptible; I, intermediate; Medium: the selective antibiotic added to the isolation culture medium. K: city of Kechida; D: downtown; CH: city of Chouhada; M: 1020 household residence; P: city of Parc à fourrage; B: city of Bouakal; H: city of Hamla; 1: 4/3/2019; 2: 18/3/2019; 3: 8/4/2019; 4: 15/4/2019; 6: 9/6/2019; 7: 3/11/2019; 8: 17/11/2019.</figcaption></figure>

2.3. *Phenotypic and Molecular Characterisation of β-Lactamase and mcr Genes*

The double disk synergy test was positive for fifteen GNB isolates, while the modified carba NP test showed positive results for five carbapenem-resistant strains, suggesting the production of carbapenemase. PCR and sequencing results showed the presence of β-lactamase genes in 26 isolates (38.80%; *n* = 26). The most prevalent β-lactamase gene, *bla*<sub>TEM</sub>, was found in seven *C. freundii* isolates (10.45%; *n* = 7), followed by the combination of *bla*<sub>CTX-M-15</sub> with *bla*<sub>TEM</sub>, which was detected in five of the positive isolates (two *C. braakii* and three *E. coli*) (7.46%; *n* = 5). Three isolates (two *K. pneumoniae* and one *C. braakii*) that were identified as CTX-M-15 type ESBL tested positive (4.48%; *n* = 3). The combination of *bla*<sub>CTX-M-15</sub> with *bla*<sub>TEM</sub> and the *bla*<sub>SHV-168</sub> genes was identified in two *K. pneumoniae* isolates (2.99%; *n* = 2). In addition, two *K. pneumoniae* isolates carried the combination of *bla*<sub>TEM</sub> with *bla*<sub>SHV-168</sub> β-lactamase gene (2.99%; *n* = 2), while one *K. pneumoniae* isolate harboured the *bla*<sub>SHV-11</sub> gene (1.49%; *n* = 1). Furthermore, one *K. pneumoniae* isolate tested positive for the combination of *bla*<sub>CTX-M-15</sub>, *bla*<sub>CTX-M-14</sub> and *bla*<sub>SHV-11</sub> (1.49%; *n* = 1). Carbapenemase genes were detected in five isolates (7.46%; *n* = 5), three of them carried *bla*<sub>OXA-48</sub> gene (two *C. braakii* and one *E. coli*); one *S. putrefaciens* carried *bla*<sub>VIM-4</sub>, and the remaining *S. putrefaciens* isolate co-harboured *bla*<sub>OXA-48</sub> and *bla*<sub>VIM-4</sub> genes. However, none of the carbapenemase genes were detected in *S. maltophilia*, *P. rhodesiae* and *P. aeruginosa* isolates, confirming the modified carba NP test results.

Out of 400 vegetable samples, the *mcr-1* gene was detected in two isolates from two coriander samples (2.99%; *n* = 2). However, no other *mcr* gene was identified. In
addition, the *P. aeruginosa* strain with imipenem resistance had modifications in its *oprD* gene sequence, with deletion at position 614 leading to a premature stop codon. The genetic characterisations of the obtained isolates are summarised in Figures 1 and 2.

## 2.4. Conjugation Experiments

Two transconjugants (TCL17 and TCL38) were successfully obtained and were resistant to amoxicillin-clavulanic acid and ertapenem or colistin (4 µg/mL), respectively. The PCR results showed that the transconjugant TCL17 harboured the *bla*<sub>OXA-48</sub> gene, and TCL38 carried the *mcr*-1-encoding gene.

## 2.5. Multi-Locus Sequence Typing

Multi-locus sequence typing analysis revealed three *E. coli* sequence types, including ST2298 in OXA-48-producing isolate, ST216 and the epidemic clone ST101 in *mcr*-1-positive strains.

# 3. Discussion

This study has demonstrated the occurrence of extended spectrum cephalosporin-, carbapenem- and colistin-resistant GNB in fresh vegetables. In recent years, there has been a significant expansion of scientific research focused on food safety. Our results correlate with other studies indicating vegetables as a possible route for the dynamic diffusion of antibiotic-resistance genes in the community [2,15]. It has been speculated that the extensive application of antibiotics on crops might augment the incidence of antibiotic resistance in bacteria existing on plant surfaces [5]. In addition to their use for the control of bacterial diseases of plants, antibiotics can be introduced into agriculture via soil fertilisation through the use of manure from animal farming. Furthermore, the soil may well represent an important source of antibiotic-resistant bacteria and genes to the plant [3]. In spite of this, soil fertilisation with animal manure is a common agricultural practice worldwide. Antibiotic-resistant bacteria associated with manure and soil may enter the plant microbiome through colonising the roots, which are in direct contact with soil, or the above ground parts of plants, possibly via the motility of root endophytes or air particulates. In this context, manure from dairy cows, usually used as a soil fertiliser, may harbour different resistant bacteria and resistance encoding genes from the gut microbiota of cattle [6]. This phenomenon has become a major problem for public health, not only in underdeveloped countries but also in high-performing, socio-economically developed countries [3]. The statistical analysis result showed that the vegetable type (source) influences the presence and rate of the targeted resistant GNB, where in our study coriander and parsley samples were the most contaminated. The types of nutrients that are frequently present on the leaf surface of leafy green vegetables can be categorised into two categories: inorganic (ions) and organic molecules (organic acids and carbohydrates). Carbohydrates (dominant phyllosphere sugars are glucose, fructose and sucrose) are of especial interest due to their capacity to readily support growth of enteric bacteria such as *E. coli* [21]. In addition, season months also showed a significant effect on the isolation rate, which is consistent with previous published studies, where the authors showed that bacterial contamination was higher for the fruit and vegetable samples during the dry seasons than in the rainy ones [22].

Most of the resistant Gram-negative bacteria species detected in the current study were previously identified on various vegetables in other studies [16]. In this study, C. freundii (20.90%; n = 14), K. pneumoniae (11.94%; n = 8), E. cloacae (11.94%; n = 8), S. maltophilia (10.45%; n = 7), E. coli (8.95%; n = 6) and C. braakii (8.95%; n = 6) were the most frequently detected. The bacterial species detected in our study have already been isolated from different sources: either environment, usually from the aquatic environments and/or soil such as R. radiobacter, A. caviae, A. eucrenophila and A. pittii, or those described as commensal and clinically relevant, such as E. coli, C. braakii, C. freundii, E. cloacae and P. mirabilis. In spite of this, the majority of the identified environmental species are also able
to cause infections, for example: *S. maltophilia*, *A. caviae*, *A. eucrenophila*, *P. alcaligenes* and *S. putrefaciens* [23–25]. The *bla*~ESBL~ and other resistance genes on vegetables have been described as existing predominantly in opportunistic and saprophytic bacteria which are thought to represent a background reservoir of antibiotic resistance genes [15,26,27]. In our study, different β-lactamase associations were reported, including *bla*~CTX-M-15~ with *bla*~TEM~, *bla*~CTX-M-15~ in combination with the *bla*~TEM~ and *bla*~SHV-168~ genes, *bla*~TEM~ in combination with the *bla*~SHV-168~ gene and the combination of the *bla*~CTX-M-15~, *bla*~CTX-M-14~ and *bla*~SHV-11~ genes; this latter combination has never been reported in Algeria. In the present study, CTX-M-15 ESBL was detected in eleven isolates and has been identified as the most common type of ESBL in Gram-negative bacteria worldwide from different sources, including humans, animal, wildlife and wastewater samples [28,29]. However, strains with no resistance genes detected in this study may harbour other known or unknown mechanisms not investigated in this study, where whole genome sequencing could be beneficial for their characterisation.

For extended spectrum cephalosporin-resistant faecal *Enterobacterales*, the main contamination sources are linked to faecally contaminated surface water used for irrigation and presumably soil contaminated with animal faeces [1]. Relatively few publications have described the presence of extended spectrum β-lactamase in GNB from vegetables. In this context, different ESBL variants have been reported on vegetable samples, including SHV, TEM, FONA, RAHN, CTX-M and OXA enzymes, where the most detected was CTX-M type [15]. To the best of our knowledge, this study has also reported the first occurrence of *bla*~VIM-4~ in *S. putrefaciens* and *bla*~OXA-48~ in *E. coli*, *C. braakii* and *S. putrefaciens* from fresh vegetables. *S. putrefaciens* is considered to be an important agent of human diseases. It commonly causes otitis media and soft tissue infections, often after trauma or exposure to water sources, bacteraemia or hepatobiliary infections and it can, in rare cases, lead to pneumonia, gastrointestinal infections and hospital-acquired infections, particularly in immunocompromised patients. It is a part of the aquatic environment microbiota, and can also be found in other sources such as soil, poultry, dairy products and medical devices [25]. Thus, the presence of such bacterial species in water and soil as their primary habitat underscores the notion that the agricultural environment may be a source of contamination by antibiotic-resistant GNB. Although the origin of *bla*~VIM-4~ and *bla*~OXA-48~ genes in vegetables is unknown and contamination may occur through animal manure fertilisation. Furthermore, the presence of carbapenem-resistant bacteria on fresh vegetables not only raises questions about the role of agriculture but could be linked to human contamination during post-harvest stages [13]. Various studies have investigated the potential sources of produce contamination in the supply chain at the post-harvest periods linked to humans, where they have reported that poor or inappropriate personal hygienic practices during transport, methods of storage, processing and preparation by handlers including sellers and consumers and the market environment also contribute to vegetable contamination [7]. The potential implications of wastewater used for irrigation in the dissemination of antibiotic-resistant human pathogens should also not be neglected as suggested by recent studies [30,31].

The first study describing a carbapenemase-producing GNB from fresh vegetable samples was in 2015, from a *Klebsiella variicola* isolate harbouring the *bla*~OXA-181~ gene in Switzerland [32]; after this report, different carbapenemase genes have been reported from vegetable samples [4,15]. More interestingly, our study also documented the presence of fourteen isolates with acquired colistin resistance, two of which were positive for the *mcr-1* gene. The remaining twelve isolates may have harboured other colistin resistance mechanisms such as the chromosomal mutation where whole genome sequencing could bring more insights into these mechanisms [33]. The two *E. coli* strains carrying the *mcr-1* gene were isolated from coriander samples. Currently, several studies have revealed that the occurrence of *mcr-1* in GNB from vegetables was lower than that from food animals [34,35]. The detection in fresh vegetables of commensal *Enterobacterales* species such as *E. coli*, known for their recent high levels of antimicrobial resistance and their involvement in
several outbreaks of foodborne diseases, is of great public health concern [13], where recently, *mcr-1*-positive *E. coli* isolates were found in different vegetable types [4,16,35,36].

In vegetables, different sequence types were identified in *E. coli* isolates carrying carbapenemase and *mcr* genes, such as ST10, ST167, ST156 and ST744 [4,13,16,35–39]; however, in our study, MLST analysis identified, for the first time, the following STs: ST2298 in OXA-48-producing isolate; ST216 and ST101 in *mcr-1*-positive isolates, respectively (Figure 3). The international high-risk resistant lineage ST101 has already been reported in wastewater isolates with ciprofloxacin resistance from an Algerian hospital [40]. It was also found in *E. coli* isolates carrying NDM-1 from municipal wastewater in Saudi Arabia [41], in NDM-1- and NDM-3-positive *E. coli* from environmental waters in Bangladesh and in a CTX-M-1-producing isolate from a wastewater treatment plant in Tunisia [42,43]. Furthermore, *mcr-1*-positive *E. coli* isolates (ST101) have been described in fish guts collected from a rainbow trout aquaculture farm in Lebanon [44]. The second detected clone, ST216, has also been reported in hospital wastewater in the United Kingdom as a carrier of the *bla*<sub>KPC</sub> resistance gene [45]. In this manner, the two detected sequence type clones of *E. coli* isolates (ST101 and ST216) are mostly associated with aquatic environments, particularly wastewater. However, no epidemiological links between the obtained isolates have been signalled. We speculate that the MCR-1-producers obtained in our study might have originated from wastewater. In addition, several studies have reported foodborne human outbreaks linked to the consumption of fresh vegetables irrigated with wastewater, and have highlighted that the type of irrigation practice plays a vital role in the contamination of farm produce [46].

<figure><img src="image_7.png" /><figcaption>**Figure 3.** MLST concatenated gene-sequence-based phylogenetic tree of our carbapenemase or *mcr*-producing *E. coli* isolates with those positive for carbapenemase or *mcr* genes occurring worldwide from vegetables. The phylogenetic tree was built using Mega 7 software with the neighbour-joining method and evolutionary distances were computed using the Kimura two-parameter method. The number above the nodes is the bootstrap level from 1000 replicates. All sequence types are presented with their relative data (study period, country, source, strain number, drug resistance determinants with the described mobile element).</figcaption></figure>
# 4. Materials and Methods

## 4.1. *Sample Collection*

During the period between March and December 2019, 400 fresh vegetable samples including zucchini (n = 82), carrot (n = 75), turnip (n = 75), cucumber (n = 70), tomato (n = 45), onion (n = 20), lettuce (n = 17), parsley (n = 8) and coriander (n = 8) were purchased from 27 different retail stores in food markets, supermarkets and street-trading greengrocers from different locations in the city of Batna, including: city of Kechida (n = 140), Parc à fourrage (n = 38), Bouakal (n = 63), Chouhada (n = 25), downtown (n = 38), 800 household residence (n = 43), 1020 household residence (n = 49) and the city of Hamla (n = 4).

The collected samples were immediately packaged aseptically in sterile polyethylene zip bags and were then transported in a cold box (4 °C) for analysis within two to three hours. Vegetable samples were first pooled before the isolation procedure where each pooled sample contained two to five samples of the same vegetable category collected from the same food market, supermarket or street-trading greengrocer and at the same period. Parsley and coriander samples were analysed individually.

## 4.2. *Isolation and Identification of β-Lactams and Colistin-Resistant GNB*

Each vegetable sample was cut into small pieces using a sterile knife and then placed in a sterile plastic bag. The samples were diluted in sterilised saline solution (0.9% NaCl) with a ratio of 1:4 (w/v). The obtained mixtures were homogenised separately and then incubated for three to four hours at 37 °C. For the isolation of extended spectrum cephalosporin-, carbapenem- and colistin-resistant GNB, each homogenised sample was inoculated in five tubes containing 10 mL of Brain-Heart Infusion broth with 64 µg/mL vancomycin and supplemented with five different selective antibiotics: 2 µg/mL cefotaxime (1), 9 µg/mL ceftazidime (2), 2 µg/mL ertapenem (3), 9 µg/mL imipenem (4) and 3 µg/mL colistin (5), respectively, and then incubated at 37 °C for 24 h. Ten microlitres from each enrichment tube was streaked onto MacConkey agar with the corresponding selective antibiotic combination and incubated at 37 °C for 24 h. *E. coli* ATCC25922 strain and strains with resistance to the used selective antibiotics were used as negative and positive controls, respectively, for quality control purposes.

Representative colonies were identified using matrix-assisted laser desorption ionisation time-of-flight mass spectrometry (MALDI-TOF-MS) (Microflex LTII, Bruker Daltonics, Bremen, Germany), as previously described [47]. Non-inoculated matrix solution was used as negative control; in addition, *Escherichia coli* DH5alpha (ref 255343, Bruker Daltonics) was used as positive control.

## 4.3. *Antimicrobial Susceptibility Testing*

The obtained isolates were subjected to antimicrobial susceptibility testing using the disc diffusion method, according to the protocols of the Antibiogram Committee of the French Society for Microbiology (CASFM 2019), (https://www.sfm-microbiologie.org/wp-content/uploads/2019/02/CASFM2019_V1.0.pdf (accessed on March 2019) [48]. The antibiotic panel presented in Figures 1 and 2 included amoxicillin (20 µg), ticarcillin (75 µg), amoxicillin-clavulanic acid (20–10 µg), ticarcillin-clavulanic acid (75–10 µg), cefoxitin (30 µg), cefotaxime (30 µg), ceftazidime (30 µg), cefepime (30 µg), aztreonam (30 µg), ertapenem (10 µg), imipenem (10 µg), tobramycin (10 µg), gentamicin (10 µg), amikacin (30 µg) and ciprofloxacin (5 µg).

The minimum inhibitory concentration (MIC) for colistin was determined by the broth microdilution method in cation-adjusted Mueller–Hinton broth, in accordance with the 2019 guidelines from the European Committee on Antimicrobial Susceptibility Testing (https://www.eucast.org/) (accessed on December 2019) [49]. Antibiogram results were interpreted according to the CA-SFM (2019) as well as the Clinical and Laboratory Standards Institute (CLSI, 2017) breakpoints. *E. coli* ATCC 25922 strain was used for quality control.
4.4. *Phenotypic and Molecular Characterisation of β-Lactamases and mcr Genes*

The double disk synergy test, using ceftazidime and cefotaxime in proximity to amoxicillin-clavulanic acid, was performed for the screening of ESBL production. The carbapenemase production of the obtained strains was verified phenotypically by the modified carba NP test [50].

Genomic DNA extraction was carried out with the EZ1 biorobot (Qiagen, Hilden, Germany), using the EZ1 DNA tissue kit (Qiagen) and the bacterial protocol card. Thereafter, the presence of ESBL determinants (*bla*~TEM~, *bla*~SHV~, *bla*~CTX-M~), carbapenemase genes (*bla*~KPC~, *bla*~VIM~, *bla*~NDM~, *bla*~OXA-48~, *bla*~OXA-23~, *bla*~OXA-24~ and *bla*~OXA-58) and *mcr* genes (*mcr*-1, *mcr*-2, *mcr*-3, *mcr*-4, *mcr*-5 and *mcr*-8) in the selected isolates was determined by real-time PCR assays performed with a CFX96 Touch real-time PCR detection system (Bio-Rad, Marnes-la-Coquette, France) using TaqMan technology. Positive qPCR strains were subjected to standard PCR and sequencing. Negative control (PCR mix) and positive control templates were included in each qPCR and conventional PCR experimental run. PCR products as well as positive and negative controls were analysed on a 1.5% SYBR Safe stained agarose gel, then sequenced using the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA, USA) on an ABI 3500xl automated sequencer (Applied Biosystems, Foster City, CA, USA) along with the positive and negative controls provided by the manufacturer. Molecular characterisation of the *oprD* gene in *P. aeruginosa* isolates was performed using PCR amplification and sequencing and then compared to the *oprD* sequence of the reference strain *P. aeruginosa* PAO1 (GenBank accession no. NC_002516.2), using Clustal W software.

The obtained sequences were analysed using CodonCode Aligner software, and then
aligned with those from the National Center for Biotechnology Information (NCBI) and
ARG-ANNOT sequence database using the BLAST program [51]. The primers and probes
used in this study are listed in Table 2.
<table><caption>**Table 2.** Oligonucleotide primers and probes used for polymerase chain reaction.</caption><thead><tr><th>Gene Name</th><th>Type of PCR</th><th>Primers</th><th>Primer Sequence (5'-&gt;3')</th><th>References</th></tr></thead><tbody><tr><td rowspan="3">bla<sub>KPC</sub></td><td rowspan="3">Real-time PCR</td><td>KPC-F</td><td>GATACCACGTTCCGTCTGGA</td><td rowspan="3">[52]</td></tr><tr><td>KPC-R</td><td>GGTCGTGTTTCCCTTTAGCC</td></tr><tr><td>KPC-Probe</td><td>6-FAM-CGCGCGCCGTGACGGA AAGC-TAMRA</td></tr><tr><td rowspan="2">bla<sub>VIM</sub></td><td rowspan="2">Real-time PCR</td><td>VIM-F</td><td>CACAGYGGCMCTTCTCGCGGAGA</td><td rowspan="2">[52]</td></tr><tr><td>VIM-R</td><td>GCGTACGTYGCCACYCCAGCC</td></tr><tr><td rowspan="3">bla<sub>NDM</sub></td><td rowspan="3">Real-time PCR</td><td>VIM-Probe</td><td>6-FAM-AGTCTCCACGCACTTTCATGA CGACCGCGTCGGCG-TAMR</td><td rowspan="3">[53]</td></tr><tr><td>NDM-1 F</td><td>GCGCAACACAGCCTGACTTT</td></tr><tr><td>NDM-1 R</td><td>CAGCCACCAAAAGCGATGTC</td></tr><tr><td rowspan="3">bla<sub>OXA-48</sub></td><td rowspan="3">Real-time PCR</td><td>NDM-1-Probe</td><td>6-FAM-CAACCGCGCCCAACTTTGGC-TAMRA</td><td rowspan="3">[52]</td></tr><tr><td>OXA48-RT-F</td><td>TCTTAAACGGGCGAACCAAG</td></tr><tr><td>OXA48-RT-R</td><td>GCGTCTGTCCATCCCACTTA</td></tr><tr><td rowspan="3">bla<sub>OXA-23</sub></td><td rowspan="3">Real-time PCR</td><td>OXA48-RT-Probe</td><td>6-FAM-AGCTTGATCGCCCTCG ATTTGG-TAMRA</td><td rowspan="3">[54]</td></tr><tr><td>OXA-23-F</td><td>TGCTCTAAGCCGCGCAAATA</td></tr><tr><td>OXA-23-R</td><td>TGACCTTTTCTCGCCCTTCC</td></tr><tr><td rowspan="3">bla<sub>OXA-24</sub></td><td rowspan="3">Real-time PCR</td><td>OXA-23-probe</td><td>FAM- GCCCTGATCGGATTGGAGAACCA-TAMRA</td><td rowspan="3">[54]</td></tr><tr><td>OXA-24-F</td><td>CAAATGAGATTTTCAAATGGGATGG</td></tr><tr><td>OXA-24-R</td><td>TCCGTCTTGCAAGCTCTTGAT</td></tr><tr><td rowspan="3">bla<sub>OXA-58</sub></td><td rowspan="3">Real-time PCR</td><td>OXA-24-probe</td><td>FAM- GGTGAGGCAATGGCATTGTCAGCA-TAMRA</td><td rowspan="3">[55]</td></tr><tr><td>OXA-58-F</td><td>CGCAGAGGGAGAATCGTCT</td></tr><tr><td>OXA-58-R</td><td>TTGCCCATCTGCCTTTTCAA</td></tr><tr><td rowspan="3">mcr-1-2</td><td rowspan="3">Real-time PCR</td><td>OXA-58-probe</td><td>FAM-GGGGAATGGCTGTAGACCCGC-TAMRA</td><td rowspan="3">[55]</td></tr><tr><td>mcr-1-2-F</td><td>CTGTGCCGTGTATGTTCAGC</td></tr><tr><td>mcr-1-2-R</td><td>TTATCCATCACGCCTTTTGAG</td></tr><tr><td rowspan="3">mcr-3</td><td rowspan="3">Real-time PCR</td><td>Probe (mcr-1-2)</td><td>FAM-TATGATGTCGATACCGCCAAATACC-TAMRA</td><td rowspan="3">[55]</td></tr><tr><td>Probe (mcr-2)</td><td>VIC-TGACCGCTTGGGTGTGGGTA-TAMRA</td></tr><tr><td>mcr-3-F</td><td>TGAATCACTGGGAGCATTAGGGC</td></tr><tr><td rowspan="3">mcr-4</td><td rowspan="3">Real-time PCR</td><td>mcr-3-R</td><td>TGCTGCAAACACGCCATATCAAC</td><td rowspan="3">[55]</td></tr><tr><td>mcr-3-probe</td><td>FAM-TGCACCGGATGATCAGACCCGT-TAMRA</td></tr><tr><td>mcr-4-F</td><td>GCCAACCAATGCTCATACCCAAAA</td></tr><tr><td rowspan="3">mcr-5</td><td rowspan="3">Real-time PCR</td><td>mcr-4-R</td><td>CCGCCCCATTCGTGAAAACATAC</td><td rowspan="3">[55]</td></tr><tr><td>mcr-4-probe</td><td>FAM-GCCACGGCGGTGTCTCTACCC-TAMRA</td></tr><tr><td>mcr-5-F</td><td>TATCCGCAAGCTACCGACGC</td></tr><tr><td rowspan="2">mcr-5</td><td rowspan="2">Real-time PCR</td><td>mcr-5-R</td><td>ACGGGCAAGCACATGATCGGT</td><td rowspan="2">[55]</td></tr><tr><td>mcr-5-probe</td><td>FAM-TGCGACACCACCGATCTGGCCA-TAMRA</td></tr></tbody></table>
<table><caption>**Table 2. Cont.**</caption><thead><tr><th>Gene Name</th><th>Type of PCR</th><th>Primers</th><th>Primer Sequence (5'-&gt;3')</th><th>References</th></tr></thead><tbody><tr><td rowspan="2">mcr-8</td><td rowspan="2">Real-time PCR</td><td>mcr-8-F</td><td>TCCGGGATGCGTGACGTTGC</td><td rowspan="2">[56]</td></tr><tr><td>mcr-8-R</td><td>TGCTGCGCGAATGAAGACGA</td></tr><tr><td rowspan="6">bla<sub>CTX-M</sub></td><td rowspan="6">Real-time PCR group A</td><td>mcr-8-probe</td><td>FAMTCATGGAGAATCGCTGGGGGAAAGC-TAMRA</td><td rowspan="6">[57]</td></tr><tr><td>CTX-A-F</td><td>CGGGCRATGGCGCARAC</td></tr><tr><td>CTX-A-R</td><td>TGCRCGGGTSGTATTGCC</td></tr><tr><td>CTX-A-probe</td><td>Yakima Yellow-CCARCGGGCGCAGYTGGTGAC-BHQ1</td></tr><tr><td>CTX-B-F</td><td>ACCGAGCCSACGCTCAA</td></tr><tr><td>CTX-B-R</td><td>CCGCTGCCGGTTTTATC</td></tr><tr><td rowspan="3">bla<sub>SHV</sub></td><td rowspan="3">Real-time PCR</td><td>CTX-B-probe</td><td>Yakima Yellow-CCCGCGYGATACCACCACGC-BHQ1</td><td rowspan="3">[57]</td></tr><tr><td>SHV-F</td><td>TCCCATGATGAGCACCTTTAAA</td></tr><tr><td>SHV-R</td><td>TCCTGCTGGCGATAGTGGAT</td></tr><tr><td rowspan="2">bla<sub>TEM</sub></td><td rowspan="2">Real-time PCR</td><td>SHV-probe</td><td>Cy5-TGCCGGTGACGAACAGCTGGAG-BBQ-650</td><td rowspan="2">[58]</td></tr><tr><td>TEM-F</td><td>GCATCTTACGGATGGCATGA</td></tr><tr><td rowspan="2">OprD</td><td rowspan="2">Standard PCR</td><td>TEM-R</td><td>GTCCTCCGATCGTTGTCAGAA</td><td rowspan="2">[58]</td></tr><tr><td>TEM-probe</td><td>6-Fam CAGTG CTGCCATAACCA TGAGTGA-BHQ-1</td></tr><tr><td rowspan="2">bla<sub>VIM</sub></td><td rowspan="2">Standard PCR</td><td>oprD-F</td><td>GGAACCTCAACTATCGCCAAG</td><td rowspan="2">[59]</td></tr><tr><td>oprD-R</td><td>GTTGCCTGTCGGTCGATTAC</td></tr><tr><td rowspan="2">bla<sub>OXA-48</sub></td><td rowspan="2">Standard PCR</td><td>VIM-F</td><td>ATTGGTCTATTTGACCGCGTC</td><td rowspan="2">[59]</td></tr><tr><td>VIM-R</td><td>TGCTACTCAACGACTGCGCG</td></tr><tr><td rowspan="2">mcr-1</td><td rowspan="2">Standard PCR</td><td>OXA-48-F</td><td>TTGGTGGCATCGATTATCGG</td><td rowspan="2">[60]</td></tr><tr><td>OXA-48-R</td><td>GAGCACTTCTTTTGTGATGGC</td></tr><tr><td rowspan="4">bla<sub>CTX-M</sub></td><td rowspan="4">Standard PCR</td><td>mcr-1-F</td><td>GCAGCATACTTCTGTGTGGTAC</td><td rowspan="4">[61]</td></tr><tr><td>mcr-1-R</td><td>TATGCACGCGAAAGAAACTGGC</td></tr><tr><td>CTX-M-1-F</td><td>CCCATGGTTAAAAAATCACTGC</td></tr><tr><td>CTX-M-1-R</td><td>CAGCGCTTTTGCCGTCTAAG</td></tr><tr><td rowspan="2">bla<sub>SHV</sub></td><td rowspan="2">Standard PCR</td><td>CTX-M-9-F</td><td>GTTACAGCCCTTCGGCGATGATTC</td><td rowspan="2">[63]</td></tr><tr><td>CTX-M-9-R</td><td>GCGCATGGTGACAAAGAGAGTGCAA</td></tr><tr><td rowspan="2"></td><td rowspan="2"></td><td>SHV-F</td><td>ATTTGTCGCTTCTTTACTCGC</td><td rowspan="2">[64]</td></tr><tr><td>SHV-R</td><td>TTTATGGCGTTACCTTTGACC</td></tr></tbody></table>
## 4.5. Conjugation Experiments

The transferability of carbapenemase and *mcr* genes was tested by conjugation experiments using sodium-azide-resistant *Escherichia coli* J53 as the recipient strain. Briefly, single colonies of the donor (L17, L21, L22, L38 and L39) and recipient isolates were inoculated separately in BHI broth and grown overnight at 37 °C. Subsequently, a ratio of 1:10 volumes (donor: recipient) was mixed and incubated overnight at 37 °C without shaking. The mixtures were then plated on nutrient agar plates supplemented with 200 µg/mL of sodium azide and 2 µg/mL of ertapenem or 3 µg/mL of colistin [65]. For all conjugation experiments, the donor strains alone and recipient strain alone were used as controls to ensure the effectiveness of the selective plates used. Antimicrobial susceptibility testing, MCNP test and OXA-48 as well as *mcr-1* PCRs were subsequently performed on the obtained transconjugants.

## 4.6. Multi-Locus Sequence Typing Analysis

Multi-locus sequence typing of *E. coli* isolates carrying carbapenemase or *mcr* genes was carried out by PCR and sequencing targeting seven housekeeping genes (*adk*, *fumC*, *gyrB*, *icd*, *mdh*, *purA* and *recA*). The sequence types were determined through the *E. coli* MLST database website (http://mlst.warwick.ac.uk/mlst/dbs/Ecoli accessed on October 2020) [66].

MLST concatenated gene-sequence-based phylogenetic tree of our carbapenemase- or *mcr*-producing *E. coli* isolates with those that carried carbapenemase and/or *mcr* genes reported worldwide (published sequence types) from vegetables was constructed in order to illustrate their phylogenetic position and clusterisation [4,13,16,35–39]. The MLST sequences of *E. coli* strains used for the phylogenetic tree were retrieved from the *E. coli* MLST database website (http://mlst.warwick.ac.uk/mlst/dbs/Ecoli accessed on January 2021). The phylogenetic tree was built using Mega 7 software [67] and evolutionary distances were computed using the Kimura two-parameter method [68].

## 4.7. Statistical Analysis

The isolation rate of the targeted drug-resistant GNB (extended spectrum cephalosporin-, carbapenem- and colistin-resistant GNB) related to each type of vegetable (source) was analysed by performing the Pearson chi-square test using SPSS (version 26.0; SPSS, Inc., Chicago, IL, USA). The level of significance was set at a *p*-value < 0.05.

# 5. Conclusions

Our results show that fresh vegetables constitute possible reservoirs for carbapenemase- and MCR-1-producing GNB and highlighted that the vegetable type as well as the season influence the presence and rate of the targeted resistant GNB. In our study, coriander and parsley samples were the most contaminated vegetable types. The presence of these bacteria in fresh vegetables is alarming and constitutes a serious human health risk. Therefore, further investigations are required for monitoring such organisms in fresh vegetables to ensure food safety and consumer health. In addition, appropriate measures such as pre-treatment of animal excrement before being used as fertilisers and the quality of irrigation water need to be taken into consideration.

**Author Contributions:** W.C. performed the sampling procedure and the experiments, analysed the data and wrote the original draft. L.L. conceptualised and directed the study and corrected the manuscript. E.B. and Z.C. performed the experiments and analysed the data. M.B. directed the study, and corrected the manuscript. J.-M.R. corrected the manuscript. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was supported by the French Government under the "Investments for the Future" program managed by the National Agency for Research (ANR) (Méditerranée-Infection10-IAHU-03).
**Data Availability Statement:** The data presented in this study are available on request from the
corresponding author.

**Acknowledgments:** This work was supported by the DGRSDT of the Algerian Ministry of Higher Education and Scientific Research. The authors are grateful to Necer Abdeldjabar for the statistical analysis.

**Conflicts of Interest:** The authors declare no conflict of interest.

References

1. Van Hoek, A.H.; Veenman, C.; Van Overbeek, W.M.; Lynch, G.; De Roda Husman, A.M.; Blaak, H. Prevalence and characterization of ESBL- and AmpC-producing Enterobacteriaceae on retail vegetables. *Int. J. Food Microbiol.* **2015**, *204*, 1–8. [CrossRef] [PubMed]

2. Holzel, C.S.; Tetens, J.L.; Schwaiger, K. Unraveling the Role of Vegetables in Spreading Antimicrobial-Resistant Bacteria: A Need for Quantitative Risk Assessment. *Foodborne Pathog. Dis.* **2018**, *15*, 671–688. [CrossRef] [PubMed]

3. Iseppi, R.; De Nu, S.; Bondi, M.; Messi, P.; Sabia, C. Extended-Spectrum β-Lactamase, AmpC, and MBL-Producing Gram-Negative Bacteria on Fresh Vegetables and Ready-to-Eat Salads Sold in Local Markets. *Microb. Drug Resist.* **2018**, *24*, 1156–1164. [CrossRef]

4. Liu, B.T.; Song, F.J. Emergence of two *Escherichia coli* strains co-harboring mcr-1 and bla (NDM) in fresh vegetables from China. *Infect. Drug Resist.* **2019**, *12*, 2627–2635. [CrossRef]

5. Stockwell, V.O.; Duffy, B. Use of antibiotics in plant agriculture. *Rev. Sci. Tech.* **2012**, *31*, 199–210. [CrossRef]

6. Zalewska, M.; Blazejewska, A.; Czapko, A.; Popowska, M. Antibiotics and Antibiotic Resistance Genes in Animal Manure—Consequences of Its Application in Agriculture. *Front. Microbiol.* **2021**, *12*, 610656. [CrossRef]

7. Berger, C.N.; Sodha, S.V.; Shaw, R.K.; Griffin, P.M.; Pink, D.; Hand, P.; Frankel, G. Fresh fruit and vegetables as vehicles for the transmission of human pathogens. *Environ. Microbiol.* **2010**, *12*, 2385–2397. [CrossRef]

8. Van, D.D.; Doi, Y. The global epidemiology of carbapenemase-producing Enterobacteriaceae. *Virulence* **2017**, *8*, 460–469.

9. Touati, A.; Mairi, A. Epidemiology of carbapenemase-producing Enterobacterales in the Middle East: A systematic review. *Expert Rev. Anti. Infect. Ther.* **2020**, *18*, 241–250. [CrossRef]

10. Lutgring, J.D. Carbapenem-resistant Enterobacteriaceae: An emerging bacterial threat. *Semin. Diagn. Pathol.* **2019**, *36*, 182–186. [CrossRef] [PubMed]

11. Son, S.J.; Huang, R.; Squire, C.J.; Leung, I.K.H. MCR-1: A promising target for structure-based design of inhibitors to tackle polymyxin resistance. *Drug. Discov. Today* **2019**, *24*, 206–216. [CrossRef] [PubMed]

12. Caniaux, I.; Van, B.A.; Zambardi, G.; Poirrel, L.; Gros, M.F. MCR: Modern colistin resistance. *Eur. J. Clin. Microbiol. Infect. Dis.* **2017**, *36*, 415–420. [CrossRef] [PubMed]

13. Liu, B.T.; Zhang, X.Y.; Wan, S.W.; Hao, J.J.; Jiang, R.D.; Song, F.J. Characteristics of Carbapenem-Resistant Enterobacteriaceae in Ready-to-Eat Vegetables in China. *Front. Microbiol.* **2018**, *9*, 1147. [CrossRef] [PubMed]

14. Anyanwu, M.U.; Jaja, I.F.; Nwobi, O.C. Occurrence and Characteristics of Mobile Colistin Resistance (mcr) Gene-Containing Isolates from the Environment: A Review. *Int. J. Environ. Res. Public Health.* **2020**, *17*, 1028. [CrossRef]

15. Chelaghma, W.; Loucif, L.; Bendahou, M.; Rolain, J.M. Vegetables and Fruit as a Reservoir of β-Lactam and Colistin-Resistant Gram-Negative Bacteria: A Review. *Microorganisms* **2021**, *9*, 2534. [CrossRef]

16. Luo, J.; Yao, X.; Lv, L.; Doi, Y.; Huang, X.; Huang, S.; Liu, J.H. Emergence of mcr-1 in *Raoultella ornithinolytica* and *Escherichia coli* Isolates from Retail Vegetables in China. *Antimicrob. Agents Chemother.* **2017**, *61*, e01139-17. [CrossRef]

17. Blaak, H.; Van Hoek, A.H.; Veenman, C.; Docters Van Leeuwen, A.E.; Lynch, G.; Van Overbeek, W.M.; De Roda Husman, A.M. Extended spectrum β-lactamase- and constitutively AmpC-producing Enterobacteriaceae on fresh produce and in the agricultural environment. *Int. J. Food Microbiol.* **2014**, *169*, 8–16. [CrossRef]

18. Touati, A.; Mairi, A.; Baloul, Y.; Lalaoui, R.; Bakour, S.; Thighilt, L.; Gharout, A.; Rolain, J.M. First detection of *Klebsiella pneumoniae* producing OXA-48 in fresh vegetables from Béjaïa city, Algeria. *J. Glob. Antimicrob. Resist.* **2017**, *9*, 17–18. [CrossRef]

19. Mesbah, Z.F.; Granier, S.A.; Touati, A.; Millemann, Y. Occurrence of Third-Generation Cephalosporins-Resistant Klebsiella pneumoniae in Fresh Fruits and Vegetables Purchased at Markets in Algeria. *Microb. Drug. Resist.* **2020**, *26*, 353–359. [CrossRef]

20. Berrazeg, M.; Drissi, M.; Medjahed, L.; Rolain, J.M. Hierarchical clustering as a rapid tool for surveillance of emerging antibiotic-resistance phenotypes in Klebsiella pneumoniae strains. *J. Med. Microbiol.* **2013**, *62*, 864–874. [CrossRef] [PubMed]

21. Mogren, L.; Windstam, S.; Boqvist, S.; Vagsholm, I.; Söderqvist, K.; Rosberg, A.K.; Lindén, J.; Mulaosmanovic, E.; Karlsson, M.; Uhlig, E.; et al. The Hurdle Approach-A Holistic Concept for Controlling Food Safety Risks Associated with Pathogenic Bacterial Contamination of Leafy Green Vegetables. A Review. *Front. Microbiol.* **2018**, *9*, 1965. [CrossRef]

22. Erhirhie, E.O.; Omoirri, M.A.; Chikodiri, S.C.; Ujam, T.N.; Kesiena, E.E.; Oseyomon, J.O. Evaluation of microbial quality of vegetables and fruits in Nigeria: A review. *Int. J. Nutr. Sci.* **2020**, *5*, 99–108.

23. Fernandez-Bravo, A.; Figueras, M.J. An Update on the Genus Aeromonas: Taxonomy, Epidemiology, and Pathogenicity. *Microorganisms* **2020**, *8*, 129. [CrossRef] [PubMed]

24. Enoch, D.A.; Birkett, C.I.; Ludlam, H.A. Non-fermentative Gram-negative bacteria. *Int. J. Antimicrob. Agents* **2007**, *29* (Suppl. 3), S33–S41. [CrossRef]

25. Ullah, S.; Mehmood, H.; Pervin, N.; Zeb, H.; Kamal, K.R.; Liaqat, S. *Shewanella putrefaciens*: An Emerging Cause of Nosocomial Pneumonia. *J. Investig. Med. High Impact Case Rep.* **2018**, *6*, 2324709618775441. [CrossRef]
26. Zurfluh, K.; Nuesch-Inderbinen, M.; Morach, M.; Zihler, B.A.; Hachler, H.; Stephan, R. Extended-spectrum-β-lactamase-producing Enterobacteriaceae isolated from vegetables imported from the Dominican Republic, India, Thailand, and Vietnam. *Appl. Environ. Microbiol.* **2015**, *81*, 3115–3120. [CrossRef]

27. Marti, R.; Scott, A.; Tien, Y.C.; Murray, R.; Sabourin, L.; Zhang, Y.; Topp, E. Impact of manure fertilization on the abundance of antibiotic-resistant bacteria and frequency of detection of antibiotic resistance genes in soil and on vegetables at harvest. *Appl. Environ. Microbiol.* **2013**, *79*, 5701–5709. [CrossRef]

28. Lahlaoui, H.; Ben Haj, K.A.; Ben Moussa, M. Epidemiology of Enterobacteriaceae producing CTX-M type extended spectrum β-lactamase (ESBL). *Med. Mal. Infect.* **2014**, *44*, 400–404. [CrossRef]

29. Peirano, G.; Pitout, J.D.D. Extended-Spectrum β-Lactamase-Producing Enterobacteriaceae: Update on Molecular Epidemiology and Treatment Options. *Drugs* **2019**, *79*, 1529–1541. [CrossRef]

30. Antwi-Agyei, P.; Cairncross, S.; Peasey, A.; Price, V.; Bruce, J.; Baker, K.; Moe, C.; Ampofo, J.; Armah, G.; Ensink, J. A Farm to Fork Risk Assessment for the Use of Wastewater in Agriculture in Accra, Ghana. *PLoS ONE* **2015**, *10*, e0142346. [CrossRef] [PubMed]

31. Makkaew, P.; Miller, M.; Fallowfield, H.J.; Cromar, N.J. Microbial risk in wastewater irrigated lettuce: Comparing *Escherichia coli* contamination from an experimental site with a laboratory approach. *Water Sci. Technol.* **2016**, *74*, 749–755. [CrossRef] [PubMed]

32. Zurfluh, K.; Poirel, L.; Nordmann, P.; Klumpp, J.; Stephan, R. First detection of *Klebsiella variicola* producing OXA-181 carbapenemase in fresh vegetable imported from Asia to Switzerland. *Antimicrob. Resist. Infect. Control.* **2015**, *4*, 38. [CrossRef] [PubMed]

33. Moffatt, J.H.; Harper, M.; Boyce, J.D. Mechanisms of Polymyxin Resistance. *Adv. Exp. Med. Biol.* **2019**, *1145*, 55–71.

34. Kumar, H.; Chen, B.H.; Kuca, K.; Nepovimova, E.; Kaushal, A.; Nagraik, R.; Bhatia, S.K.; Dhanjal, D.S.; Kumar, V.; Kumar, A.; et al. Understanding of Colistin Usage in Food Animals and Available Detection Techniques: A Review. *Animals* **2020**, *10*, 1892. [CrossRef]

35. Oh, S.S.; Song, J.; Kim, J.; Shin, J. Increasing prevalence of multidrug-resistant mcr-1-positive *Escherichia coli* isolates from fresh vegetables and healthy food animals in South Korea. *Int. J. Infect. Dis.* **2020**, *92*, 53–55. [CrossRef]

36. Manageiro, V.; Jones-Dias, D.; Ferreira, E.; Caniça, M. Plasmid-Mediated Colistin Resistance (mcr-1) in *Escherichia coli* from Non-Imported Fresh Vegetables for Human Consumption in Portugal. *Microorganisms* **2020**, *8*, 429. [CrossRef]

37. Zurfuh, K.; Poirel, L.; Nordmann, P.; Nuesch-Inderbinen, M.; Hachler, H.; Stephan, R. Occurrence of the Plasmid-Borne mcr-1 Colistin Resistance Gene in Extended-Spectrum-β-Lactamase-Producing Enterobacteriaceae in River Water and Imported Vegetable Samples in Switzerland. *Antimicrob. Agents Chemother.* **2016**, *60*, 2594–2595. [CrossRef]

38. Wang, J.; Yao, X.; Luo, J.; Lv, L.; Zeng, Z.; Liu, J.H. Emergence of *Escherichia coli* co-producing NDM-1 and KPC-2 carbapenemases from a retail vegetable, China. *J. Antimicrob. Chemother.* **2018**, *73*, 252–254. [CrossRef] [PubMed]

39. Liu, B.T.; Li, X.; Zhang, Q.; Shan, H.; Zou, M.; Song, F.J. Colistin-Resistant mcr-Positive Enterobacteriaceae in Fresh Vegetables, an Increasing Infectious Threat in China. *Int. J. Antimicrob. Agents* **2019**, *54*, 89–94. [CrossRef]

40. Anssour, L.; Messai, Y.; Estepa, V.; Torres, C.; Bakour, R. Characteristics of ciprofloxacin-resistant Enterobacteriaceae isolates recovered from wastewater of an Algerian hospital. *J. Infect. Dev. Ctries.* **2016**, *10*, 728–734. [CrossRef] [PubMed]

41. Mantilla-Calderon, D.; Jumat, M.R.; Wang, T.; Ganesan, P.; Al-Jassim, N.; Hong, P.Y. Isolation and Characterization of NDM-Positive *Escherichia coli* from Municipal Wastewater in Jeddah, Saudi Arabia. *Antimicrob. Agents Chemother.* **2016**, *60*, 5223–5231. [CrossRef] [PubMed]

42. Ben Said, L.; Jouini, A.; Alonso, C.A.; Klibi, N.; Dziri, R.; Boudabous, A.; Ben Slama, K.; Torres, C. Characteristics of extended-spectrum β-lactamase (ESBL)- and pAmpC beta-lactamase-producing Enterobacteriaceae of water samples in Tunisia. *Sci. Total. Environ.* **2016**, *550*, 1103–1109. [CrossRef]

43. Toleman, M.A.; Bugert, J.J.; Nizam, S.A. Extensively drug-resistant New Delhi metallo-β-lactamase-encoding bacteria in the environment, Dhaka, Bangladesh, 2012. *Emerg. Infect. Dis.* **2015**, *21*, 1027–1030. [CrossRef]

44. Hassan, J.; Eddine, R.Z.; Mann, D.; Li, S.; Deng, X.; Saoud, I.P.; Kassem, I.I. The Mobile Colistin Resistance Gene, mcr-1.1, Is Carried on IncX4 Plasmids in Multidrug Resistant *E. coli* Isolated from Rainbow Trout Aquaculture. *Microorganisms* **2020**, *8*, 1636. [CrossRef]

45. Decraene, V.; Phan, H.T.T.; George, R.; Wyllie, D.H.; Akinremi, O.; Aiken, Z.; Cleary, P.; Dodgson, A.; Pankhurst, L.; Crook, D.W.; et al. A Large, Refractory Nosocomial Outbreak of *Klebsiella pneumoniae* Carbapenemase-Producing *Escherichia coli* Demonstrates Carbapenemase Gene Outbreaks Involving Sink Sites Require Novel Approaches to Infection Control. *Antimicrob. Agents Chemother.* **2018**, *62*, e01689-18. [CrossRef] [PubMed]

46. Adegoke, A.A.; Amoah, I.D.; Stenstrom, T.A.; Verbyla, M.E.; Mihelcic, J.R. Epidemiological Evidence and Health Risks Associated with Agricultural Reuse of Partially Treated and Untreated Wastewater: A Review. *Front. Public Health* **2018**, *6*, 337. [CrossRef] [PubMed]

47. Seng, P.; Drancourt, M.; Gouriet, F.; La Page, S.B.; Fournier, P.E.; Rolain, J.M.; Raoult, D. Ongoing revolution in bacteriology: Routine identification of bacteria by matrix-assisted laser desorption ionization time-of-flight mass spectrometry. *Clin. Infect. Dis.* **2009**, *49*, 543–551. [CrossRef]

48. Available online: https://www.sfm-microbiologie.org/wp-content/uploads/2019/02/CASFM2019_V1.0.pdf (accessed on March 2019).

49. Available online: https://www.eucast.org (accessed on December 2019).
50. Bakour, S.; Garcia, V.; Loucif, L.; Brunel, J.M.; Gharout-Sait, A.; Touati, A.; Rolain, J.M. Rapid identification of carbapenemase-producing Enterobacteriaceae, *Pseudomonas aeruginosa* and *Acinetobacter baumannii* using a modified Carba NP test. *New Microbes New Infect.* **2015**, *7*, 89–93. [CrossRef]

51. Gupta, S.K.; Padmanabhan, B.R.; Diene, S.M.; Lopez-Rojas, R.; Kempf, M.; Landraud, L.; Rolain, J.M. ARG-ANNOT, a new bioinformatic tool to discover antibiotic resistance genes in bacterial genomes. *Antimicrob. Agents Chemother.* **2014**, *58*, 212–220. [CrossRef]

52. Yousfi, H.; Hadjadj, L.; Dandachi, I.; Lalaoui, R.; Merah, A.; Amoura, K.; Dahi, A.; Dekhil, M.; Messalhi, N.; Diene, S.M.; et al. Colistin- and Carbapenem-Resistant *Klebsiella pneumoniae* Clinical_Isolates: Algeria. *Microb. Drug. Resist.* **2019**, *25*, 258–263. [CrossRef]

53. Diene, S.M.; Bruder, N.; Raoult, D.; Rolain, J.M. Real-time PCR assay allows detection of the New Delhi metallo-β-lactamase (NDM-1)-encoding gene in France. *Int. J. Antimicrob. Agents* **2011**, *37*, 544–546. [CrossRef] [PubMed]

54. Mesli, E.; Berrazeg, M.; Drissi, M.; Bekkhoucha, S.N.; Rolain, J.M. Prevalence of carbapenemase-encoding genes including New Delhi metallo-β-lactamase in *Acinetobacter* species, Algeria. *Int. J. Infect. Dis.* **2013**, *17*, e739–e743. [CrossRef] [PubMed]

55. Ly, T.D.A.; Hadjadj, L.; Hoang, V.T.; Louni, M.; Dao, T.L.; Badiaga, S.; Tissot-Dupont, H.; Raoult, D.; Rolain, J.M.; Gautret, P. Low prevalence of resistance genes in sheltered homeless population in Marseille, France, 2014–2018. *Infect. Drug. Resist.* **2019**, *12*, 1139–1151. [CrossRef] [PubMed]

56. Nabti, L.Z.; Sahli, F.; Ngaiganam, E.P.; Radji, N.; Mezaghcha, W.; Lupande-Mwenebitu, D.; Baron, S.A.; Rolain, J.M.; Diene, S.M. Development of real-time PCR assay allowed describing the first clinical *Klebsiella pneumoniae* isolate harboring plasmid-mediated colistin resistance mcr-8 gene in Algeria. *J. Glob. Antimicrob. Resist.* **2020**, *20*, 266–271. [CrossRef]

57. Roschanski, N.; Fischer, J.; Guerra, B.; Roesler, U. Development of a multiplex real-time PCR for the rapid detection of the predominant beta-lactamase genes CTX-M, SHV, TEM and CIT-type AmpCs in Enterobacteriaceae. *PLoS ONE* **2014**, *9*, e100956. [CrossRef]

58. Hong, J.S.; Yoon, E.J.; Lee, H.; Jeong, S.H.; Lee, K. Clonal Dissemination of *Pseudomonas aeruginosa* Sequence Type 235 Isolates Carrying blaIMP-6 and Emergence of blaGES-24 and blaIMP-10 on Novel Genomic Islands PAGI-15 and -16 in South Korea. *Antimicrob. Agents Chemother.* **2016**, *60*, 7216–7223. [CrossRef]

59. Jeon, B.C.; Jeong, S.H.; Bae, I.K.; Kwon, S.B.; Lee, K.; Young, D.; Lee, J.H.; Song, J.S.; Lee, S.H. Investigation of a nosocomial outbreak of imipenem-resistant *Acinetobacter baumannii* producing the OXA-23 beta-lactamase in korea. *J. Clin. Microbiol.* **2005**, *43*, 2241–2245. [CrossRef]

60. Mathlouthi, N.; Al-Bayssari, C.; El Salabi, A.; Bakour, S.; Ben Gwierif, S.; Zorgani, A.A.; Jridi, Y.; Ben Slama, K.; Rolain, J.M.; Chouchani, C. Carbapenemases and extended-spectrum β-lactamases producing Enterobacteriaceae isolated from Tunisian and Libyan hospitals. *J. Infect. Dev. Ctries.* **2016**, *10*, 718–727. [CrossRef]

61. Chabou, S.; Leangapichart, T.; Okdah, L.; Le Page, S.; Hadjadj, L.; Rolain, J.M. Real-time quantitative PCR assay with Taqman® probe for rapid detection of MCR-1 plasmid-mediated colistin resistance. *New Microbes New Infect.* **2016**, *13*, 71–74. [CrossRef]

62. Carattoli, A.; Garcia-Fernandez, A.; Varesi, P.; Fortini, D.; Gerardi, S.; Penni, A.; Mancini, C.; Giordano, A. Molecular epidemiology of *Escherichia coli* producing extended-spectrum beta-lactamases isolated in Rome, Italy. *J. Clin. Microbiol.* **2008**, *46*, 103–108. [CrossRef]

63. Yu, Y.; Ji, S.; Chen, Y.; Zhou, W.; Wei, Z.; Li, L.; Ma, Y. Resistance of strains producing extended-spectrum beta-lactamases and genotype distribution in China. *J. Infect. 2007*, *54*, 53–57. [CrossRef] [PubMed]

64. Yagi, T.; Kurokawa, H.; Shibata, N.; Shibayama, K.; Arakawa, Y. A preliminary survey of extended-spectrum beta-lactamases (ESBLs) in clinical isolates of *Klebsiella pneumoniae* and *Escherichia coli* in Japan. *FEMS Microbiol. Lett.* **2000**, *184*, 53–56. [PubMed]

65. Bendjama, E.; Loucif, L.; Chelaghma, W.; Attal, C.; Bellakh, F.Z.; Benaldjia, R.; Kahlat, I.; Meddour, A.; Rolain, J.M. First detection of an OXA-48-producing *Enterobacter cloacae* isolate from currency coins in Algeria. *J. Glob. Antimicrob. Resist.* **2020**, *23*, 162–166. [CrossRef]

66. Available online: http://mlst.warwick.ac.uk/mlst/dbs/Ecoli (accessed on October 2020).

67. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. *Mol. Biol. Evol.* **2016**, *33*, 1870–1874. [CrossRef] [PubMed]

68. Kimura, M.A. simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. *J. Mol. Evol.* **1980**, *16*, 111–120. [CrossRef] [PubMed]