The First Isolation and Molecular Characterization of Shiga Toxin-Producing Virulent Multi-Drug Resistant Atypical Enteropathogenic Escherichia coli O177 Serogroup From South African Cattle

Introduction

Enteropathogenic E. coli (EPEC) is a group of diarrhoeagenic E. coli that is reported to cause high morbidity and mortality in humans, especially in immune-compromised subjects, elderly individuals and young children. EPEC are characterized by the presence of intimin (eae) genes coupled with the absence of the stx genes (Martins et al., 2016; Alonso et al., 2017). The eae gene is responsible for attaching and effacing (A/E) lesions on the intestinal epithelial cell of the host (Kaper et al., 2004; Martins et al., 2016; Malik et al., 2017). Based on the presence or absence of the EPEC adherence factor (EAF) plasmid, EPEC is subdivided into two groups that include typical Enteropathogenic E. coli (tEPEC) and atypical Enteropathogenic E. coli (aEPEC) (Trabulsi et al., 2002; Alonso et al., 2017). The tEPEC strains possess EAF plasmid, which encodes a bundle forming pili (bfpA) while aEPEC strains lack the bfpA gene (Canizalez-Roman et al., 2013; Malik et al., 2017). It is on these bases that the virulence potentials of aEPEC is poorly understood and highly questioned. Despite the fact that tEPEC have been most often associated disease complications in humans, it is only recently that aEPEC been reported to cause diseases in both animals and humans (Malik et al., 2017). This may account for the reason why previous studies that have been documented worldwide and in the study area have focused on EHEC, especially E. coli O157 and non-O157 strains that received great attention due to its high pathogenicity (Ateba and Bezuidenhout, 2008; Ateba and Mbewe, 2011; Iwu et al., 2016; Jajarmi et al., 2017; Toro et al., 2018). The findings of most of these studies were in agreement with previous reports indicating that EHEC strains possessed stx, stx₂, hlyA, and eaeA gene determinants that enhance their potential to cause infections such as diarrhea, hemorrhagic colitis (HC) and hemolytic uremic syndrome (HUS) in humans (Farrokh et al., 2013). However, to the best of our knowledge there is currently no report on the occurrence of EPEC, and aEPEC in particular among South African food-producing animals. This therefore creates a knowledge gap on the virulence profiles of aEPEC strains in the area. Despite the fact that aEPEC strains were known to be less pathogenic when compared to EHEC counterparts, data from some recent studies have revealed the presence of virulence determinants in aEPEC strains thus making them to start receiving attention as pathogens of severe clinical significance in humans (Beutin et al., 2005; Bibbal et al., 2017) as well as in food-producing animals, especially sheep and wildlife (Otero et al., 2013; Álvarez-Suárez et al., 2016; Martins et al., 2016). In this study, we expand our previous investigations on aEPEC in the area by determining the occurrence of aEPEC O177 strains in cattle and also provide an overview of the putative virulence and antimicrobial resistance profiles of the isolates. This was motivated from the fact that cattle are the primary reservoir of both EHEC (stx⁺) and aEPEC (stx⁻) strains, thus providing opportunities for horizontal transfer of virulence and antibiotic resistant genes between species and/or strains living in the same ecological niches (Bibbal et al., 2017).

Materials and Methods

Samples Collection and Isolation of E. coli O177 Serogroup

A total of 780 cattle fecal samples were collected from eight commercial dairy and beef cattle farms in the North West Province, South Africa. Ethical clearance was obtained from the North West University Research Ethics Committee prior to the commencement of the study. Sampling was done between January 2017 and July 2017. The selection of the farms was based on the production system; either intensive, semi-intensive, and/or extensive farming system in the study area. The selection of all the three systems was based on the fact their management system is different. Furthermore, this was to avoid biasness in terms of selecting the farms. Ethical procedures such as restraining the animals using proper facilities and equipment were enforced during the collection of the samples. Fecal samples were collected directly from the rectum of individual animals using sterile arm-length gloves and in order to avoid duplication of sampling, the cattle were looked in their respective handling pens immediately after collection. Samples were placed in sterile sample collection bottles, labeled appropriately and immediately transported on ice to the Antimicrobial Resistance and Phage Biocontrol Laboratory (AREPHABREG), in the Department of Microbiology, North-West University for microbial analysis. Data on antibiotic type and treatment history were collected for the purpose of understanding antibiotic exposure histories of isolates from the study population.

Approximately, 2 g of fecal samples were dissolved in 10 mL of 10% (w/v) saline solution and homogenized. Aliquots of 5 μL were transferred into 10 mL buffered peptone water (Oxoid Ltd., Basingstok, Hampshire, UK). Ten-fold serial dilutions were prepared and aliquots of 100 μL from each dilution was spread-plated on Rainbow agarO157 plates (Kase et al., 2015) obtained from Biolog Inc., USA. The plates were incubated aerobically at 37°C for 24 h. Colonies with purple and pink colors were randomly picked and purified by streaking on sorbitol MacConkey agar (Merck, S.A) supplemented with 1 mg/L potassium tellurite (Merck, SA) and incubated aerobically at 37°C for 24 h. Pure isolates were preserved in 20 % (v/v) glycerol (Merck, SA) and stored at −80°C for future use.

Genomic DNA Extraction From Presumptive Isolates

Overnight culture for each sample was prepared and genomic DNA was extracted from the cultures using the Zymo Research Genomic DNA^TM-Tissue MiniPrep Kit (Biolab, South Africa) following the manufacturer's instructions. DNA concentration and purity was determined using the NanoDrop Lite 1,000 spectrophotometer (model: Thermo Fisher Scientific, USA) and the DNA was stored at −80°C until further analysis by PCR.

Identification of E. coli O177 Serogroup Using Multiplex PCR Assay

A singleplex PCR assay targeting the uidA E. coli O177-specific gene sequence was performed using a previous protocol (Anbazhagan et al., 2011). Escherichia coli O177 serogroup-specific primer pairs were designed using the Primer3 software (this study), based on rmlB and wzy gene sequences that encode for dTDP-glucose 4, 6-dehydratase, and O-antigen polymerase, respectively (Ye et al., 2012). The specificity of the PCR primers was tested using National Center for Biotechnology Information/Primer-Basic Alignment Search Tool (NCBI/Primer-BLAST) (https://www.ncbi.nlm.nih.gov/tools/primer-blast/). The oligonucleotides were synthesized and supplied by Inqaba Biotechnical Industry Ltd., Pretoria, South Africa. The newly developed multiplex PCR protocol was empirically validated for its specificity, sensitivity reproducibility, and robustness. PCR assays was performed to amplify rmlB and wzy gene fragments and the primer sequences, targeted genes, amplicon sizes as well as the PCR conditions are listed in Table 1. The PCR reactions constituted 12.5 μL of 2X DreamTaq Green Master Mix, 0.5 μM of each primer and 1 μL of template DNA. The volume of the reaction mixture was adjusted to 25 μL with RNase-nuclease free PCR water. A non-DNA template (nuclease-free water) reaction tube served as a negative control while a DNA sample from E. coli O177 (isolated during the preliminary study) was used as positive control. All the PCR reagents used were Fermentas USA products supplied by Inqaba Biotec Industry Ltd., Sunnyside, Pretoria, South Africa. Amplifications were performed using DNA thermal cycler (model- Bio-Rad C1000 Touch^TM Thermal Cycler, Singapore). All the PCR products were held at 4°C until gel electrophoresis.

TABLE 1

Table 1. Oligonucleotide primers used for amplification of the various targeted genes in E. coli O177 serogroup.

Detection of Virulence Genes in E. coli O177 Isolates

E. coli O177 isolates were subjected to a bfp gene PCR analysis in order to screen for characters of aEPEC strains. A highly virulent environmental E. coli O157:H7 isolate (bfp⁺) previously isolated from our research group (AREPHABREG) was used as positive control. The positive control E. coli O157:H7 possessed the stx₁, stx₂, eaeA, and hlyA genes. In addition, DNA extracted from an E. coli (ATCC 98222) that is a non-pathogenic strain was included in each PCR run as a negative control. Isolates that were positive for the bfp gene were further screened for the presence of an array of STEC virulence genes that included eaeA, hlyA, stx₁, stx₂, and stx₂ variants (stx_2a, stx_2b, stx_2c, stx_2d, stx_2e, stx_2f, stx_2g) using previously described protocols (Paton and Paton, 1998; He et al., 2012). Primer sequences, target genes, amplicon sizes as well as PCR cycling conditions for the different genes are listed in Table 1. All the PCR reactions were prepared as 25 μL standard volumes comprising 12.5 μL of 2X DreamTaq Green Master Mix, 0.5 μM of each primer, 1 μL of template DNA and RNase-nuclease free PCR water. Amplifications were performed using DNA thermal cycler (model- Bio-Rad C1000 Touch^TM Thermal Cycler, Singapore). All the PCR products were held at 4°C until they were separated by electrophoresis.

Antimicrobial Susceptibility Test

The antimicrobial susceptibility profile of the isolates was determined using the Kirby-Bauer disc diffusion technique (Bauer et al., 1966), making use of antibiotic discs obtained from Mast Diagnostics, UK. The antibiotics using in this study comprised Ampicillin (AMP), 10 μg; Chloramphenicol (C), 30 μg; Erythromycin (E), 15 μg; Kanamycin (K), 30 μg; Norfloxacin (NOR), 10 μg; Streptomycin (S), 10 μg and Tetracycline (TE), 30 μg. Some of these antimicrobial agents were selected due to the fact that they are widely used as prophylactics in both beef and dairy cattle farming in South Africa. Plates were incubated aerobically at 37°C for 24 h and antibiotic growth inhibition zone diameter data were compared with standard reference values in order to classify the isolates as sensitive, intermediate resistance or resistant to a particular antibiotic (Clinical Laboratory Standards Institute, 2016). Escherichia coli ATCC 25922 was used as a reference for quality control in antimicrobial susceptibility test.

Detection of Genetic Determinants for Antibiotic Resistance Genes by PCR

All confirmed E. coli O177 isolates that were resistant to three or more antibiotics were designated multi-drug resistant isolates and were screened for the presence of the tetA, tetW, aadA, strA, strB, ampC, cmlA, ermB, and kan antibiotic resistance determinants. Primer sequences, target genes, amplicon sizes as well as PCR cycling conditions for the different genes are listed in Table 2. All the PCR reactions were prepared as 25 μL standard volumes comprising 12.5 μL of 2X DreamTaq Green Master Mix, 0.5 μM of each primer, 1 μL of template DNA and RNase-nuclease free PCR water. Amplifications were performed using DNA thermal cycler (model- Bio-Rad C1000 Touch^TM Thermal Cycler, Singapore). All the PCR products were kept at 4°C and later resolved by electrophoresis.

TABLE 2

Table 2. Oligonucleotide primers used for amplification of the various targeted antibiotic resistance genes in E. coli O177 serogroup.

Agarose Gel Electrophoresis

All PCR amplicons were resolved by electrophoresis on a 2% (w/v) agarose gel containing 0.001 μg/mL ethidium bromide. A 100 bp DNA molecular weight DNA marker (Fermentas, USA) was included in each gel and used to confirm the sizes of the amplicons. A ChemiDoc Imaging System (Bio-Rad ChemiDoc^TM MP Imaging System, UK) was used to capture the images using Gene Snap software, version 6.0022.

Sequence Analysis

The amplified rmlB and wzy gene fragments were sequenced by Inqaba Biotec, Pretoria, South Africa, and sequences were subjected to a Blast Search Tool (https://blast.ncbi.nlm.nih.gov/Blast.cgi) in order to confirm the identities of the isolates (Altschul et al., 1997). Nucleotide sequence data were further deposited into GenBank using the BankIt NCBI web-based sequence submission tool (https://www.ncbi.nlm.nih.gov/genbank/) in order to obtain their accession numbers.

Statistical Analysis

Measured parameters were tested for normality using the NORMAL option in the Proc Univariate statement prior to analysis of variance. Data on the proportions of samples positive for E. coli O177 serogroup were analyzed for effect of farming system using the general linear models (GLM) procedures of SAS (2010) according to the following statistical model:

$\begin{array}{l}{Y}_{ij}=\mu +F_i+E_{ij}\end{array}$

where, Y_ij, response variable, μ, overall mean, F_i, farming system (intensive, semi-intensive and extensive) effect and E_ij, random error associated with observation ij, assumed to be normally and independently distributed. Data on number of isolates carrying virulent genes, number of antibiotic resistant isolates, number of isolates carrying antibiotic resistance genes, and number of multidrug resistant isolates were square root-transformed before statistical analysis using the GLM procedures (SAS, 2010). For all statistical tests, significance was declared at p < 0.05.

Results

Identification of E. coli O177 Serogroup Using Multiplex PCR Analysis

A total of 780 cattle fecal samples were analyzed for the presence of E. coli O177 using Rainbow agarO157 and a total of 1,272 non-repetitive presumptive isolates were obtained. Amplification of the uidA E. coli genus-specific gene sequence revealed that a large proportion (915; 71.93%) of the isolates were successfully identified as E. coli isolates. All the 915 E. coli isolates were further screened for the presence of E. coli O177 serogroup specific rmlB and wzy gene fragments using multiplex PCR analysis. Out of 915 isolates screened, 376 (41.09%) were identified as E. coli O177 serogroup and Figure 1 indicates a 2% (w/v) agarose gel of representative rmlB and wzy gene fragments amplified in the study. No significant (p > 0.05) differences were observed on the occurrence of E. coli O177 serogroup across intensive, semi-intensive, and extensive animal production systems. The occurrence of E. coli O177 serogroup was 33.3 ± 14.43, 33.3 ± 14.43, and 50.0 ± 14.18% in intensive semi-intensive extensive farming systems, respectively.

FIGURE 1

Figure 1. A 2% agarose gel image of the rmlB (645 bp) and wzy (457 bp) gene fragments amplified from E. coli O177 isolates. Lane M = 100 bp molecular marker; Lanes 1–15 = rmlB and wzy gene fragments, respectively, that were amplified from E. coli O177 isolates; Lane PC = rmlB and wzy gene fragments amplified from E. coli O177 serotype; Lane NC = negative control sample with DNA from E. coli O157:H7 environmental isolate.

Detection of Virulence Genotypes in E. coli O177 Isolates

A total of 376 atypical E. coli O177 isolates that were positively confirmed by PCR were screened for the presence of five virulence genes; bfpA, eaeA, hlyA, stx₁, and stx₂, using PCR analysis. None of the isolates possessed the bfpA gene fragment and thus were classified as aEPEC strains. On the contrary, all the other four virulence genes (eaeA, hlyA, stx₁, and stx₂) were successfully detected in this study. Generally, 350 isolates harbored these virulence genes. There were significant differences (p < 0.001) in terms of occurrence of various virulence genes in E. coli O177 isolates. The hlyA (12.74%) and stx₂ (11.20%) were the most commonly detected genes (p < 0.001), followed by stx₁ (9.07%), and eaeA (7.25%), Figure 2. Two stx₂ subtypes namely; the stx_2a and stx_2d, were detected among the isolates. A few of the isolates (2.60%) harbored stx_2a, while 0.63% possessed stx_2dgene. Some of the isolates carried a combination of the genes detected. The majority (9.74%) of the isolates carried a combination of hlyA/stx₂ while 7.40, 7.09, 6.80, and 5.55% possessed stx₁/stx₂, stx₁/stx₂/hlyA, hlyA/eaeA, and stx₂/hlyA/eaeA, respectively (p < 0.001). There was no difference (p > 0.05) between the occurrences of stx₁/eaeA (6.35%) and stx₁/stx₂/hlyA/eaeA (5.11%) in E. coli O177 isolates. Only one isolate possessed stx₁/stx₂/stx_2a/stx_2d/hlyA/eaeA. Despite the fact that none of the isolates possessed the stx₂ subtypes stx_2b, stx_2c, stx_2e, stx_2f, and stx_2g, the virulence gene profiles of these aEPEC strains especially the isolate with the stx₂ subtypes was a cause for concern.

FIGURE 2

Figure 2. Distribution of virulence genes in E. coli O177 serogroup isolated from cattle feces. The bars indicate the standard error (p < 0.05).

Antimicrobial Resistance Profiles

Antimicrobial susceptibility tests of the 376 isolates revealed varied antimicrobial resistance profiles against all the seven antibiotics tested. Most of the E. coli O177 isolates were resistant to erythromycin (63.84%) compared to other antibiotics (p < 0.001), Figure 3. No difference (p > 0.05) was observed in E. coli O177 resistance to ampicillin (21.54%), tetracycline (13.37%), streptomycin (17.01%), kanamycin (2.42%), chloramphenicol (1.97%), and norfloxacin (1.40%). The isolates showed resistance to at least two antibiotics tested. Most of the isolates (20.74%) exhibited high level of resistance to three or more antibiotics, Figure 4 (p < 0.001).

FIGURE 3

Figure 3. Antibiotic resistance pattern of E. coli O177 serogroup isolated from cattle feces. The bars indicate the standard error (p < 0.05).

FIGURE 4

Figure 4. AMP, ampicillin; C, chloramphenicol; E, erythromycin; K, kanamycin; NOR, norfloxacin; TE, tetracycline; S, streptomycin; MDR, multidrug-resistant; XDR, extensively drug-resistant. Multiple resistance patterns of E. coli O177 serogroup obtained from cattle feces. The bars indicate the standard error (p < 0.05).

Proportion of Isolates With Antimicrobial Resistance Genes

A total of 73 E. coli O177 isolates harbored either one or more different antimicrobial resistance genes investigated, however, only five of the ten genes were detected in the isolates. Resistance genes associated with streptomycin and erythromycin were the most frequently detected. Genes aadA (22.86%), strA (19.85%), strB (19.38%) ermB (19.68%), and tetA (18.23%) were detected. However, there was no significant differences in the occurrence of these genes in E. coli O177 isolates. No isolates possessed the tetW gene. Some isolates harbored a combination of three or more antibiotic resistance genes. Most (24%) of the isolates carried a combination of the aadA, strA, and strB genes while a smaller proportion 11% isolates possessed all the five antimicrobial resistance genes (aadA/strA/strB/ermB/tetA) (p < 0.05).

Sequence Identifier and Accession Numbers

Sequence data analysis of the rmlB and wzy gene fragments amplified from chromosomal DNA of isolates revealed a high percentage similarity (97%) sequence homology to E. coli serogroup O177. The nucleotide sequences for Seq1, Seq2, and Seq3 were assigned GenBank Accession numbers; MH389799, MH389800, and MH389801, respectively.

Discussion

In this study, the occurrence of aEPEC, E. coli O177 serogroup in cattle was investigated. Although the virulence profiles of aEPEC strains is poorly understood, recent reports suggest increased incidence of human infections caused by aEPEC strains worldwide (Ingle et al., 2016; Martins et al., 2016). Cattle are primary reservoir of E. coli species and contaminated meat as well as their associated food products have been implicated in foodborne infections in humans (Ateba and Mbewe, 2011). Studies have reported the occurrence of aEPEC strain in food-producing animals, foods of animal origin and farm environments (Horcajo et al., 2012; Otero et al., 2013; Ghanbarpour et al., 2017). Escherichia coli was successfully isolated from cattle feces using Rainbow agar O157. The identities of the isolates were confirmed as E. coli through amplification of uidA gene fragments. The rmlB and wzy O-antigen gene clusters specific for E. coli O177 were successfully amplified through Multiplex PCR analysis. The PCR assay was highly specific and reproducible for the detection of E. coli O177 serogroup from cattle feces. The occurrence of E. coli O177 in cattle was higher (41.09 %) in this study when compared to a previous report which detected this serogroup in ovine (Martins et al., 2016). Despite this, there were no significant (p > 0.05) differences on the occurrence of E. coli O177 serogroup across all the three farming systems.

Given the fact that aEPEC strains are known to lack medically important virulence factors, their ability to cause diseases particularly in humans is not yet well-understood (Ingle et al., 2016). Previous reports indicated that aEPEC strains do not harbor stx and bfpA genes. However, recent reports have revealed that aEPEC strains especially those from ruminants, harbor virulence gene profiles similar to those of the EHEC strains (Horcajo et al., 2012). In addition, some studies have reported that aEPEC strains may represent LEE-positive Shiga toxin-producing E. coli (STEC) that have lost the toxin-encoding prophage, or tEPEC that have lost the genes for BFP indicating the potential to have evolved from other pathogenic E. coli strains (Tennant et al., 2009). A motivation was the fact that previous studies conducted in the study area revealed a large proportion of virulent EHEC strains (Ateba and Bezuidenhout, 2008; Ateba and Mbewe, 2011, 2013, 2014). Despite the fact that we were unable at this stage to determine if aEPEC isolated from humans in South African cattle fall into either of these categories, an assessment of the bfpA and other primary STEC virulence genes (eaeA, hlyA, stx₁, stx₂, and stx₂ variants) was performed. Data of the shiga toxin subtypes among non-O157 STEC serogroups in cattle may provide an indication of the potential health risks these aEPEC isolates may cause humans. The bfpA gene was not detected in all the E. coli isolates analyzed and thus confirmed that the isolates belonged to aEPEC group, emerging diarrheagenic pathogens in both developing and developed countries (Ingle et al., 2016; González et al., 2017).

All the major virulence genes (stx₁, stx₂, hlyA, and eaeA) were detected in this study. Similar findings were reported by Momtaz et al. (2013a) and Momtaz et al. (2013b). Generally, the majority of E. coli O177 isolates possessed virulence genes. The hlyA was the most frequently detected virulence gene in this study. These results are similar to those reported in other studies (Ateba and Bezuidenhout, 2008; Ateba and Mbewe, 2011; Momtaz et al., 2013c), which described high prevalence of hlyA gene in E. coli O157 isolated from cattle, beef, vegetable and water intended for human consumption. The high prevalence of hlyA gene may be attributed to the fact that hlyA is a plasmid encoded gene and as a result, it may be transferred between strains of the same species (Ateba and Mbewe, 2011). The hlyA gene encodes α-hemolysin, a toxin that lyses mammalian erythrocytes (Toro et al., 2018).

Another observation was that the E. coli O177 isolates possessed stx genes. These genes are considered as primary virulence genes in shiga toxin-producing E. coli and their presence in aEPEC strains, especially E. coli O177 serogroup raises a serious public health concern. The stx₂ was the most prevalent gene detected in STEC/EHEC (Ateba and Mbewe, 2011; Rantsiou et al., 2012). Despite the fact that low frequency of stx₂ gene was observed in this study, the occurrence of this gene was significantly higher than stx₁ (p < 0.001). This observation was similar to the findings of the previous studies (Beutin et al., 2005; Jajarmi et al., 2017; Toro et al., 2018). Interestingly, stx₂ positive isolates also harbored stx_2a and stx_2d gene subtypes. However, the other stx₂ subtypes were not detected. Despite this, the presence of stx₂, stx_2a, and stx_2d genes in E. coli O177 serogroup is a serious concern. Strains harboring these genes have been implicated in hemolytic colitis and hemorrhagic uremic syndrome infections in humans (Farrokh et al., 2013). It is also reported that a strain harboring stx₂ gene is more virulent than a strain carrying either stx₁ or both stx₁ and stx₂ (Farrokh et al., 2013; Toro et al., 2018). Moreover, stx_2a and stx_2d subtypes are associated with severe HC and HUS infections in human (Farrokh et al., 2013; Cha et al., 2018). Therefore, the presence of stx₂gene and its subtypes in E. coli O177 may increase the spectrum of infections in humans and thus creating a public health concern.

The stx₁ and eaeA genes were also detected but at low levels compared to hlyA and stx₂ gene fragments. This was in contrast with the previous studies, which reported high occurrence of stx₁ and eaeA genes in E. coli isolated from poultry meat and humans (Momtaz et al., 2013b; Hemmatinezhad et al., 2015). Despite low detection of these genes, especially E. coli O177 serogroup, these findings cannot be overemphasized. The strain carrying stx₁ may cause diarrhea in immunocompromised individuals (Farrokh et al., 2013). In addition, the eaeA gene encodes for the outer membrane which mediates adherence between the STEC and/or EPEC and intestinal epithelial cells (Kaper et al., 2004; Farrokh et al., 2013). Although eaeA was detected in this study, this gene was not detected in some isolates. This was surprising because all the isolates were negative for bfpA gene and thus classified as aEPEC strains. These results differ significantly with other studies where the eaeA gene has been detected in aEPEC strains (Trabulsi et al., 2002; Beutin et al., 2005; Otero et al., 2013; Martins et al., 2016; Ghanbarpour et al., 2017; González et al., 2017). The possible explanation for this observation could be that the isolates may have lost the gene during sub-culturing process (Karch et al., 1992). Given that the strains carrying eaeA gene are considered potentially pathogenic, absence of this gene in some of the isolates in this study should not be underestimated. The eaeA-negative isolates may use other genes such as saa, aidA, agn43, ehaA, or iha to adhere to the epithelial cell of the host (Toro et al., 2018). However, the presence of these genes was not investigated in this study.

The co-expression of stx₁, stx₂, hlyA, and eaeA genes may increase the pathogenicity of E. coli strains (Jajarmi et al., 2017). It was observed that E. coli O177 isolates possessed combinations of virulence genes. There were twelve different combinations of virulence genes detected in this study. The most frequently detected combinations were stx₂/hlyA (9.74 %), stx₁/stx₂ (7.40%), stx₁/stx₂/hlyA (7.09%), and hlyA/eaeA (6.80%). Only 5.11% of the isolates harbored all the four genes (eaeA/hlyA/stx₁/stx₂). It was also remarkable that one isolate possessed all the virulence (stx₁/stx₂/hlyA/eaeA) genes and two stx₂ (stx_2a and stx_2d genes) subtypes. Simultaneous detection of these virulence genes in E. coli O177 serogroup in cattle pose a public health concern. Furthermore, these findings are similar to those reported in E. coli O26 in a previous study (Jajarmi et al., 2017; Ranjbar et al., 2017). The combinations of virulence genes in this study are higher than those reported by the previous studies (Ateba and Mbewe, 2011; Jajarmi et al., 2017; Toro et al., 2018). However, there was low occurrence of combination of stx₂ subtypes as compared to a previous study (Jajarmi et al., 2017).

Antimicrobial agents primarily play a vital role in the lives of both humans and animals worldwide (Hudson et al., 2017). South Africa, is one of the countries where the usage of antimicrobial agents in animal production is very high (Hudson et al., 2017). In addition, the overuse of antimicrobial agents in food-producing animals, often without either professional consultation or supervision, has resulted in the emergence of antimicrobial resistant bacteria and antibiotic resistant genes in human and animal pathogens (Qiao et al., 2017). Moreover, previous investigations in the study area revealed the presence of antimicrobial resistant gene determinants in E. coli O157 isolated from cattle, beef, pig, pork, vegetables, and water intended for human consumption (Ateba and Bezuidenhout, 2008). Despite this, farmers continue to use antimicrobial agents to maximize production and this presents severe public health challenges.

In this study, E. coli O177 isolates revealed phenotypic resistance against all the seven classes of antimicrobial agents tested. In contrary to the previous study which reported high prevalence of resistance against ampicillin (100%), gentamicine (100%), and tetracycline (96.87%) Ranjbar et al. (2018), this study revealed resistance against erythromycin (63.84%), ampicillin (21.54%), tetracycline (13.37%), streptomycin (17.01%), and kanamycine (2.42%). These disparities could be attributed to geographical location and the sample type. Furthermore, 20.74% of the isolates were resistant to at least three and antimicrobial agents tested. These results were similarly to those reported in the previous study (Ranjbar et al., 2017). Antimicrobial resistance genes encoding for three antibiotics (erythromycin, streptomycin, and tetracycline) were detected. It has been reported that antibiotic resistance is common among E. coli isolates obtained from animals and food of animal origin due to frequent use of antibiotics in animals (Ryu et al., 2012). Generally, all isolates harbored antibiotic resistance genes. It is also worth mentioning that the same isolates harbored STEC virulence genes and previous studies have shown that isolates with similar genetic determinants may pose severe complications in humans (Ateba and Bezuidenhout, 2008; Hemeg, 2018). The occurrence of antibiotic resistance genes detected in this study was higher as compared to the previous study (Ryu et al., 2012). This could be due to the relatively higher usage of antimicrobial agents in the intensive livestock farming system where the isolates were obtained.

Antibiotic resistance genes associated with some of the antibiotics tested in this study were detected. The aadA, strA, and strB were the most frequently detected genes. However, there were significant differences among all the genes. Most of the isolates possessed aadA (22.86%), strA (19.85%), and strB (19.38%) genes that code for resistance to streptomycin and similar observations have been reported in other studies (Srinivasan et al., 2007; Shahrani et al., 2014; Bibbal et al., 2017). Despite the fact that tet genes are commonly found in E. coli from cattle due to the high usage of tetracycline in dairy and feedlot cattle, the occurrence of tet resistance genes in E. coli O177 isolates was low. Only the tetA gene was detected individually or in combination with other genes. However, in a previous report (Olowe et al., 2013), the proportion of tetA (43.8%) was higher than tetB (32.0%) and a combination of both tetA and tetB (4.4%) among clinical E. coli isolates. In addition, Ranjbar et al. (2018) reported tetA (76.56%), tetB (20.31%), cat1 (18.75%) and cmlA (1.50%) in E. coli O157 and non-O157 serogroups from milk. Another study reported high detection of tetA and tetB in E. coli isolated from urinary infection patients, ruminants and donkey raw milk (Momtaz et al., 2012, 2013c). Similar to a previous report (Olowe et al., 2013), none of the isolates in this study harbored the tetW, amp, cmlA, and kan genes.

In conclusion, to the best of our knowledge this is the first study to report the occurrence of E. coli O177 strain in cattle, especially in South Africa. In this study, the search for virulence determinants in aEPEC clearly revealed that a majority of strains harbored DNA sequences that encode known virulence-associated determinants of other pathogenic E. coli such as O157 and non-O157 strains. In addition, E. coli O177 strains from this study displayed high levels of resistance to aminoglycoside antibiotic group thus have the potential to pose a public health concern to humans. Indeed, the fact that all aEPEC strains in this study expressed a variety of tEPEC virulence determinants supports the case for continued efforts to conduct a large scale study designed to determine the phylogenetic homology among aEPEC isolates from different sources. Furthermore, whole genome sequencing of E. coli O177 is required in order to understand the complete pathogenic determinants of this serogroup.

Data Availability Statement

The datasets generated for this study can be found in the NCBI.

Author Contributions

VM and CA contributed to project conception and contributed reagents, materials, and analysis tools. PM, VM, and CA designed the experiments. PM performed the experiments. All authors contributed to manuscript revision and read and approved the submitted version.

Funding

This research was supported financially by the National Research Foundation (Grant No. 112543) and the North-West University Postgraduate Bursary.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors wish to thank Mr. BJ Morapedi and Dr. A Kumar for their assistance with this project.

References

Alonso, C. A., Mora, A., Díaz, D., Blanco, M., González-Barrio, D., Ruiz-Fons, F., et al. (2017). Occurrence and characterization of stx and/or eae-apositive Escherichia coli isolated from wildlife, including a typical EPEC strain from a wild boar. Vet. Microbiol. 207, 69–73. doi: 10.1016/j.vetmic.2017.05.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Altschul, S. F., Madden, T. L., Schäffer, A. A., Zhang, J., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402. doi: 10.1093/nar/25.17.3389

PubMed Abstract | CrossRef Full Text | Google Scholar

Álvarez-Suárez, M. E., Otero, A., García-López, M. L., Dahbi, G., Blanco, M., Mora, A., et al. (2016). Genetic characterization of Shiga toxin-producing Escherichia coli (STEC) and atypical enteropathogenic Escherichia coli (EPEC) isolates from goat's milk and goat farm environment. Int. J. Food Microbiol. 236, 148–154. doi: 10.1016/j.ijfoodmicro.2016.07.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Anbazhagan, D., Mui, W. S., Mansor, M., Yan, G. O., Yusof, M. Y., and Sekaran, S. D. (2011). Development of conventional and real-time multiplex PCR assays for the detection of nosocomial pathogens. Braz. J. Microbiol. 42, 448–458. doi: 10.1590/S1517-83822011000200006

PubMed Abstract | CrossRef Full Text | Google Scholar

Ateba, C. N., and Bezuidenhout, C. C. (2008). Characterisation of Escherichia coli O157strains from humans, cattle, and pigs in the North-West Province, South Africa. Int. J. Food Microbiol. 128, 181–188. doi: 10.1016/j.ijfoodmicro.2008.08.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Ateba, C. N., and Mbewe, M. (2011). Detection of Escherichia coli O157: H7 virulence genes in isolates from beef, pork, water, human, and animal species in the northwest province, South Africa: public health implications. Res. Microbiol. 162, 240–248. doi: 10.1016/j.resmic.2010.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Ateba, C. N., and Mbewe, M. (2013). Determination of the genetic similarities of fingerprints from Escherichia coli O157:H7 isolated from different sources in the North West Province, South Africa using ISR, BOXAIR, and REP-PCR analysis. Microbiol. Res. 168, 438–446. doi: 10.1016/j.micres.2013.02.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Ateba, C. N., and Mbewe, M. (2014). Genotypic characterization of Escherichia coli O157:H7isolates from different sources in the north-west province, South Africa, using enterobacterial repetitive intergenic consensus PCR analysis. Int. J Mol. Sci. 15, 9735–9747. doi: 10.3390/ijms15069735

PubMed Abstract | CrossRef Full Text | Google Scholar

Bauer, A. W., Kirby, W. M., Sherris, J. C., and Turck, M. (1966). Antibiotic susceptibility testing by a standardized single disk method. Am. J. Clin. Pathol. 45, 493–496. doi: 10.1093/ajcp/45.4_ts.493

PubMed Abstract | CrossRef Full Text | Google Scholar

Bergeron, S., Boopathy, R., Nathaniel, R., Corbin, A., and LaFleur, G. (2015). Presence of antibiotic resistant bacteria and antibiotic resistance genes in raw source water and treated drinking water. Int. Biodeter. Biodegrad. 102, 370–374. doi: 10.1016/j.ibiod.2015.04.017

CrossRef Full Text | Google Scholar

Beutin, L., Kong, Q., Feng, L., Wang, Q., Krause, G., Leomil, L., et al. (2005). Development of PCR assays targeting the genes involved in synthesis and assembly of the new Escherichia coli O174 and O177 O antigens. J. Clin. Microbiol. 43, 5143–5149. doi: 10.1128/JCM.43.10.5143-5149.2005

CrossRef Full Text | Google Scholar

Bibbal, D., Um, M. M., Kérourédan, M., Dupouy, V., Toutain, P. L., Bousquet-Mélou, A., et al. (2017). Mixing of Shiga toxin-producing and enteropathogenic Escherichia coli in a wastewater treatment plant receiving city and slaughterhouse wastewater. Int. J. Hyg. Environ. Health 221, 355–363. doi: 10.1016/j.ijheh.2017.12.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Canizalez-Roman, A., Gonzalez-Nuñez, E., Vidal, J. E., Flores-Villaseñor, H., and León-Sicairos, N. (2013). Prevalence and antibiotic resistance profiles of diarrheagenic Escherichia coli strains isolated from food items in northwestern Mexico. Int. J. Food Microbiol. 164, 36–45. doi: 10.1016/j.ijfoodmicro.2013.03.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Cha, W., Fratamico, P. M., Ruth, L. E., Bowman, A. S., Nolting, J. M., Manning, S. D., et al. (2018). Prevalence and characteristics of Shiga toxin-producing Escherichia coli in finishing pigs: Implications on public health. Int. J. Food Microbiol. 264, 8–15. doi: 10.1016/j.ijfoodmicro.2017.10.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Clinical and Laboratory Standards Institute (2016). Performance Standards for Antimicrobial Susceptibility Testing, 26th Edn. CLSI supplement M100S. Wayne, PA: Clinical and Laboratory Standards Institute.

Google Scholar

Farrokh, C., Jordan, K., Auvray, F., Glass, K., Oppegaard, H., Raynaud, S., et al. (2013). Review of Shiga-toxin-producing Escherichia coli (STEC) and their significance in dairy production. Int. J. Food Microbiol. 162, 190–212. doi: 10.1016/j.ijfoodmicro.2012.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghanbarpour, R., Askari, N., Ghorbanpour, M., Tahamtan, Y., Mashayekhi, K., Afsharipour, N., et al. (2017). Genotypic analysis of virulence genes and antimicrobial profile of diarrhoeagenic Escherichia coli isolated from diseased lambs in Iran. Trop. Anim. Health Prod. 49, 591–597. doi: 10.1007/s11250-017-1234-7

CrossRef Full Text | Google Scholar

González, J., Cadona, J. S., Sanz, M., Bustamante, A. V., and Sanso, A. M. (2017). Molecular characterization of diarrheagenic Escherichia coli isolated from vegetables in Argentina. Int. J. Food Microbiol. 261, 57–61. doi: 10.1016/j.ijfoodmicro.2017.09.021

PubMed Abstract | CrossRef Full Text | Google Scholar

He, X., Quiñones, B., Mcmahon, S., and Mandrell, R. E. (2012). A single-step purification and molecular characterization of functional Shiga-toxin 2 variants from pathogenic Escherichia coli. Toxins 4, 487–504. doi: 10.3390/toxins4070487

PubMed Abstract | CrossRef Full Text | Google Scholar

Hemeg, H. A. (2018). Molecular characterization of antibiotic resistant Escherichia coli isolates recovered from food samples and outpatient Clinics, AKA. Saudi J. of Biol. Sci. 25, 928–931. doi: 10.1016/j.sjbs.2018.01.016

CrossRef Full Text | Google Scholar

Hemmatinezhad, B., Khamesipour, F., Mohammadi, M., Safarpoor Dehkordi, F., and Mashak, Z. (2015). Microbiological investigation of O-serogroups, virulence factors and antimicrobial resistance properties of shiga toxin-producing Escherichia coli isolated from ostrich, turkey and quail meats. J. Food Safety 35, 491–500. doi: 10.1111/jfs.12199

CrossRef Full Text | Google Scholar

Horcajo, P., Domínguez-Bernal, G., de La Fuente, R., Ruiz-Santa-Quiteria, J. A., Blanco, J. E., Blanco, M., et al. (2012). Comparison of ruminant and human attaching and effacing Escherichia coli (AEEC) strains. Vet. Microbiol. 155, 341–348. doi: 10.1016/j.vetmic.2011.08.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Hudson, J. A., Frewer, L. J., Jones, G., Brereton, P. A., Whittingham, M. J., and Stewart, G. (2017). The agri-food chain and antimicrobial resistance: a review. Trends Food Sci. Techn. 69, 131–147. doi: 10.1016/j.tifs.2017.09.007

CrossRef Full Text | Google Scholar

Ingle, D. J., Tauschek, M., Edwards, D. J., Hocking, D. M., Pickard, D. J., Azzopardi, K. I., et al. (2016). Evolution of atypical enteropathogenic E. coli by repeated acquisition of LEE pathogenicity island variants. Nat Microbiol. 1:15010. doi: 10.1038/nmicrobiol.2015.10

PubMed Abstract | CrossRef Full Text | Google Scholar

Iwu, C. J., Iweriebor, B. C., Obi, L. C., and Okoh, A. I. (2016). Occurrence of non-O157 Shiga toxin-producing Escherichia coli in two commercial swine farms in the Eastern Cape Province, South Africa. Comp. Immunol. Microbiol. Infect. Dis. 44, 48–53. doi: 10.1016/j.cimid.2015.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Jajarmi, M., Fooladi, A. A., Badouei, M. A., and Ahmadi, A. (2017). Virulence genes, Shiga toxin subtypes, major O-serogroups, and phylogenetic background of Shiga toxin-producing Escherichia coli strains isolated from cattle in Iran. Microb. Pathog. 109, 274–279. doi: 10.1016/j.micpath.2017.05.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaper, J. B., Nataro, J. P., and Mobley, H. L. (2004). Pathogenic Escherichia coli. Nature Rev. Microbiol. 2, 123–140. doi: 10.1038/nrmicro818

PubMed Abstract | CrossRef Full Text | Google Scholar

Karch, H., Meyer, T., Rüssmann, H., and Heesemann, J. (1992). Frequent loss of Shiga-like toxin genes in clinical isolates of Escherichia coli upon subcultivation. Infect. Immun. 60, 3464–3467.

PubMed Abstract | Google Scholar

Kase, J. A., Maounounen-Laasri, A., Son, I., Lin, A., and Hammack, T. S. (2015). Comparison of eight different agars for the recovery of clinically relevant non-O157 Shiga toxin-producing Escherichia coli from baby spinach, cilantro, alfalfa sprouts and raw milk. Food Microbiol. 46, 280–287. doi: 10.1016/j.fm.2014.08.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Malik, A., Nagy, B., Kugler, R., and Szmolka, A. (2017). Pathogenic potential and virulence genotypes of intestinal and faecal isolates of porcine post-weaning enteropathogenic Escherichia coli. Res. Vet. Sci. 115, 102–108. doi: 10.1016/j.rvsc.2017.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Martins, F. H., Guth, B. E., Piazza, R. M., Elias, W. P., Leão, S. C., Marzoa, J., et al. (2016). Lambs are an important source of atypical enteropathogenic Escherichia coli in southern Brazil. Vet. Microbiol. 196, 72–77. doi: 10.1016/j.vetmic.2016.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Momtaz, H., Dehkordi, F., Hosseini, M. J., Sarshar, M., and Heidari, M. (2013b). Serogroups, virulence genes and antibiotic resistance in Shiga toxin-producing Escherichia coli isolated from diarrheic and non-diarrheic pediatric patients in Iran. Gut pathogens 5:39. doi: 10.1186/1757-4749-5-39

PubMed Abstract | CrossRef Full Text | Google Scholar

Momtaz, H., Dehkordi, F. S., Rahimi, E., Ezadi, H., and Arab, R. (2013a). Incidence of Shiga toxin-producing Escherichia coli serogroups in ruminant's meat. Meat Sci. 95, 381–388. doi: 10.1016/j.meatsci.2013.04.051

PubMed Abstract | CrossRef Full Text | Google Scholar

Momtaz, H., Farzan, R., Rahimi, E., Safarpoor Dehkordi, F., and Souod, N. (2012). Molecular characterization of Shiga toxin-producing Escherichia coli isolated from ruminant and donkey raw milk samples and traditional dairy products in Iran. Sci. World J. 2012:231342. doi: 10.1100/2012/231342

PubMed Abstract | CrossRef Full Text | Google Scholar

Momtaz, H., Karimian, A., Madani, M., Dehkordi, F., Ranjbar, R., Sarshar, M., et al. (2013c). Uropathogenic Escherichia coli in Iran: serogroup distributions, virulence factors and antimicrobial resistance properties. Ann. Clin. Microbiol. Antimicrob. 12:8. doi: 10.1186/1476-0711-12-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Olowe, O. A., Idris, O. J., and Taiwo, S. S. (2013). Prevalence of TET genes mediating tetracycline resistance in Escherichia coli clinical isolates in Osun State, Nigeria. Eur. J. Microbiol. Immunol. 3, 135–140. doi: 10.1556/EuJMI.3.2013.2.7

PubMed Abstract | CrossRef Full Text | Google Scholar

Otero, V., Rodríguez-Calleja, J.-M., Otero, A., García-López, M.-L., and Santos, J. A. (2013). Genetic characterization of atypical enteropathogenic Escherichia coli isolates from ewes' milk, sheep farm environments, and humans by multilocus sequence typing and pulsed-field gel electrophoresis. Appl. Environ. Microbiol. 79, 5864–5869. doi: 10.1128/AEM.01809-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Paton, A. W., and Paton, J. C. (1998). Detection and characterization of shiga toxigenic Escherichia coli by using multiplex PCR assays for stx1, stx2, eaeA, enterohemorrhagic E. coli hlyA, rfbO111, and rfbO157. J. Clin. Microbiol. 36, 598–602.

PubMed Abstract | Google Scholar

Qiao, M., Ying, G. G., Singer, A. C., and Zhu, Y. G. (2017). Review of antibiotic resistance in China and its environment. Environ. Int. 110, 160–172. doi: 10.1016/j.envint.2017.10.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Ranjbar, R., Dehkordi, F. S., Shahreza, M. H., and Rahimi, E. (2018). Prevalence, identification of virulence factors, O-serogroups and antibiotic resistance properties of Shiga-toxin producing Escherichia coli strains isolated from raw milk and traditional dairy products. Antimicrob. Resist. Infect. Control 7:53. doi: 10.1186/s13756-018-0345-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ranjbar, R., Masoudimanesh, M., Dehkordi, F. S., Jonaidi-Jafari, N., and Rahimi, E. (2017). Shiga (Vero)-toxin producing Escherichia coli isolated from the hospital foods; virulence factors, o-serogroups and antimicrobial resistance properties. Antimicrob. Resist. Infect. Control 6:4. doi: 10.1186/s13756-016-0163-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Rantsiou, K., Alessandria, V., and Cocolin, L. (2012). Prevalence of Shiga toxin-producing Escherichia coli in food products of animal origin as determined by molecular methods. Int. J. Food Microbiol. 154, 37–43. doi: 10.1016/j.ijfoodmicro.2011.12.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Ryu, S. H., Lee, J. H., Park, S. H., Song, M. O., Park, S. H., Jung, H. W., et al. (2012). Antimicrobial resistance profiles among Escherichia coli strains isolated from commercial and cooked foods. Int. J. Food Microbiol. 159, 263–266. doi: 10.1016/j.ijfoodmicro.2012.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Samra, Z. Q., Naseem, M., Khan, S. J., Nadia, D., and Athar, M. A. (2009). PCR targeting of antibiotic resistant bacteria in public drinking water of Lahore metropolitan, Pakistan. BioMed. Environ. Sci. 22, 458–463. doi: 10.1016/S0895-3988(10)60002-5

PubMed Abstract | CrossRef Full Text | Google Scholar

SAS (2010). Users Guide: Statistics. Version 9.4. Cary, NC: SAS Institute.

Google Scholar

Shahrani, M., Dehkordi, F. S., and Momtaz, H. (2014). Characterization of Escherichia coli virulence genes, pathotypes and antibiotic resistance properties in diarrheic calves in Iran. Biol. Res. 47:28. doi: 10.1186/0717-6287-47-28

PubMed Abstract | CrossRef Full Text | Google Scholar

Srinivasan, V., Gillespie, B. E., Lewis, M. J., Nguyen, L. T., Headrick, S. I., Schukken, Y. H., et al. (2007). Phenotypic and genotypic antimicrobial resistance patterns of Escherichia coli isolated from dairy cows with mastitis. Vet. Microbiol. 124, 319–328. doi: 10.1016/j.vetmic.2007.04.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Tennant, S. M., Tauschek, M., Azzopardi, K., Bigham, A., Bennett-Wood, V., Hartland, E. L., et al. (2009). Characterisation of atypical enteropathogenic E. coli strains of clinical origin. BMC Microbiol. 9:117. doi: 10.1186/1471-2180-9-117

PubMed Abstract | CrossRef Full Text | Google Scholar

Toro, M., Rivera, D., Jiménez, M. F., Díaz, L., Navarrete, P., and Reyes-Jara, A. (2018). Isolation and characterization of non-O157 Shiga toxin-producing Escherichia coli (STEC) isolated from retail ground beef in Santiago, Chile. Food Microbiol. 75, 55–60. doi: 10.1016/j.fm.2017.10.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Trabulsi, L. R., Keller, R., and Gomes, T. A. (2002). Typical and atypical Enteropathogenic Escherichia coli. Emerg. Infect. Dis. 8, 508–513. doi: 10.3201/eid0805.010385

PubMed Abstract | CrossRef Full Text | Google Scholar

Ye, J., Coulouris, G., Zaretskaya, I., Cutcutache, I., Rozen, S., and Madden, T. L. (2012). Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformat. 13:134. doi: 10.1186/1471-2105-13-134

PubMed Abstract | CrossRef Full Text | Google Scholar