- Research article
- Open Access
Identification of Escherichia coli from broiler chickens in Jordan, their antimicrobial resistance, gene characterization and the associated risk factors
BMC Veterinary Researchvolume 15, Article number: 159 (2019)
Avian pathogenic Escherichia coli (APEC) is the principle cause of colibacillosis affecting poultry. The main challenge to the poultry industry is antimicrobial resistance and the emergence of multidrug resistant bacteria that threaten the safety of the food chain. Risk factors associated with emergence of antimicrobial resistance among avian pathogenic E. coli were correlated with the inappropriate use of antimicrobials along with inadequate hygienic practices, which encourages the selection pressure of antimicrobial resistant APEC. The aim of this study was to isolate, identify, serogroup and genotype APEC from broilers, assess their antibiotic resistance profile, expressed genes and the associated risk factors.
APEC was isolated from the visceral organs of sick chickens with a prevalence of 53.4%. The most prevalent serotypes were O1, O2, O25 and O78, in percentage of 14.8, 12.6, 4.4 and 23.7%, respectively. Virulence Associated Genes; SitA, iss, iucD, iucC, astA, tsh cvi and irp2 were detected in rate of 97.4, 93.3, 75, 74, 71, 46.5, 39 and 34%, respectively and 186 (69.2%) isolates possess > 5–10 genes. The highest resistance was found against sulphamethoxazole-trimethoprim, florfenicol, amoxicillin, doxycycline and spectinomycin in percentage; 95.5, 93.7, 93.3, 92.2 and 92.2%, respectively. Sixty-eight percent of APEC isolates were found to have at least 5 out of 8 antimicrobial resistant genes. The most predominant genes were Int1 97%, tetA 78.4%, bla TEM 72.9%, Sul1 72.4%, Sul2 70.2%. Two risk factors were found to be associated with the presence of multi-drug resistant APEC in broiler chickens, with a P value ≤0.05; the use of ground water as source of drinking water and farms located in proximity to other farms.
This study characterized the VAGs of avian pathogenic E. coli and establish their antimicrobial resistance patterns. The widespread of antimicrobial resistance of APEC isolates and detection of ARGs highlighted the need to monitor the spread of ARGs in poultry farms and the environment in Jordan. Use of ground water and closely located farms were significant risk factors associated with the presence of MDR APEC in broiler chickens in Jordan.
Avian pathogenic E. coli causes localized or systemic infection outside the avian gut, which indicates as Extraintestinal Pathogenic E. coli (ExPEC). The infection caused by ExPEC is termed colibacillosis which is an infectious disease characterized by acute fatal septicemia or sub- acute fibrinous pericarditis, airsacculitis, salpingitis, and peritonitis affect broiler chickens aged 4–6 weeks [1, 2]. Colibacillosis is a common bacterial disease of economic importance in poultry through decreasing the infected birds’ productivity, increase mortality, condemnation of infected carcasses at slaughter, and prophylaxis and treatment cost  and is reported worldwide.
APEC is considered a primary or secondary pathogen of poultry. Strains which carry virulence genes (adhesin, invasins, toxins, resistance to host serum, iron acquisition systems, temperature-sensitive hemagglutinin, and K1 capsule) have all been shown to contribute to APEC pathogenesis [3, 4] and could induce colibacillosis without previous immune suppression factors; stress or concurrent infections .
The control and prevention of bacterial diseases in food animals is achieved by the application of antimicrobials during the periods of high risk of infectious bacterial diseases, as prophylactic treatment, and as growth promoters .
Bacterial antimicrobial resistance develops naturally over time; the unprecedented increase of antimicrobial resistant organisms is linked to the massive use of antimicrobial agents for disease control and prevention in human and animal medicine . Several forces play a role in the spread of antimicrobial resistant bacteria include the presence of carrier animal moving between animal herds and through vector action .
The key points in controlling avian colibacillosis are management interventions, infections control and vaccination strategies . Wide range of antimicrobial agents is used in poultry colibacillosis treatment, which include: β-lactams (penicillins, cephalosporin), aminoglycosides, tetracycline, sulphonamides and fluoroquinolones . The frequent use of antimicrobial agents give rise to selective pressure that lead to antimicrobial resistance against APEC .
The development of resistance is a complex process associated with the presence of resistance encoding genes that are found inside plasmids or chromosomal genetic material. Integrons are the genetic material responsible for capturing resistance genes that spread via the genetic mobile elements; transposons and plasmid. The presence of integrons is detected by amplification of integrase genes (intI 1, intI2 and intI 3) . Resistance to tetracycline is mediated through efflux pump system which encoded by tetracycline resistance group of genes (tetA, tetB, tetC, tetD, tetE and tetG) . Phenicols resistance encoding genes are (cat1, cat2, cat3, cmlA and cmlB)  aminoglycosides resistance genes are (strA, strB, addA1, addA 2)  and genes responsible for sulphonamide resistance are (sul 1, sul 2 and sul 3) .
Antimicrobial resistant E. coli strains pose a serious problem for public health, since these strains could be passed to humans via the food chain or by direct contact with infected birds. In addition, resistant E. coli may act as transporters for antimicrobial resistant genes to other pathogens .
In many developed countries, administration of antimicrobial agents is not only restricted for treatment purpose. Antimicrobials can also be used to enhance animal productivity, feed conversion rate and growth rate in food producing animals . This type of farming practice allows antimicrobial drugs to eliminate sensitive bacterial strains and select strains with genetic traits that can resist antimicrobials, which provides favourable conditions for selected strain persistence and spread at the farm level .
The use of antimicrobial agents as feed additives, administered at low concentrations (sub-therapeutic dose) usually over long periods of time, may lead to development of resistance [19, 20]. Other risk factors include: the breed of the animal, dose, duration of treatment, capacity of the farm, and animal husbandry practices . Poor hygiene and lack of commitment with control measures and disease prevention have participated in the propagation and expansion of antimicrobial resistant strains .
Resistant bacteria could be shed in the faeces and passed into sewage systems, which are considered as suitable transporters for resistance genes and the spread of resistant bacteria into the wider environment. Antibiotic residues and by-products found in municipal sewage, waste water treatment plants, and soil, are flushed into rivers by surface water and reach ground water resources .
The use of disinfectants to limit infection transmission between animals subsequently increases animal health and productivity. Quaternary ammonium compounds (QACs) may have the potential to induce the emergence of antimicrobial resistance, which could be raised from cross-resistance between QACs and a range of antimicrobials [24, 25]. The use of chicken litter-based organic fertilizers in the presence of antimicrobial resistance pathogens are considered as a serious environmental hazard, as the spread of fertilizers on pasture could contaminate ground water sources and land that may facilitate the transmission of antimicrobial resistant pathogens to other animal species and humans. This highlights that proper waste management could be effective in controlling the spread of antimicrobial resistance pathogens [21, 26]. Antimicrobial resistance has also been reported in wildlife, indicating that the common habitat between wildlife, food animals, water sources and environmental contamination has resulted in the transmission of antimicrobial resistant bacterial pathogens into the food chain as well as their role in contaminating foods of plant origin .
Therefore, the objectives of the current study are to isolate and identify E. coli from live sick birds, establish their serotypes, their virulence associated genes, antibiotic resistance profiles and their associated genes and to identify risk factors and farming practice associated with the antimicrobial resistance E. coli.
E. coli isolation
A total of 504 broiler chicken samples (from 84 broilers farm) were cultured, 269 (53.4%) isolates were confirmed as E. coli by conventional and RapID™ ONE System and were used for further molecular and antimicrobial testing.
Molecular identification of E. coli by PCR
All isolates that were confirmed as E. coli by the RapID™ ONE system also underwent PCR to further confirm the isolates as E. coli. The universal primer for 16 s RNA with 585 bp band size was used. Escherichia coli ATCC 25922 was used as positive control as demonstrated by (Fig. 1).
All confirmed E. coli isolates were serotyped. One hundred eighty-nine (70.3%) were identified as eleven different serotypes using the available antisera; O1, O2, O9, O18, O25, O26, O78, O111, O114, O119, O127. Whereas, the remaining isolates; 54 (20%) were untypeable and 26 (9.66%) were rough strains that show autoagglutination, serotypes and their frequencies are shown in (Table 1).
Multiplex polymerase chain reaction method for detection of virulence associated genes (VAGs)
Sixteen virulence associated genes were investigated using multiplex PCR, for avian E. coli indicates that sitA is the most prevalent gene (262, 97.4%) followed by iss (251, 93.3%), iucC (199, 74%), iucD (203, 75%), astA (190, 71%), tsh (125, 46.5%), cvi (106, 39%), irp2 (91, 34%), KpsII (33, 12.3%), KPS (20, 7.4%), KpsIII (13, 4.8%) and vat (7, 2.6%). HlyD and ibeA were not detected and papC and sfa were detected in one isolate each among the 269 E. coli tested (Fig. 2a, b).
One hundred eighty-six (69.2%) of the 269 E. coli tested isolates possess > 5–10 VAGs. In detail; 3 isolates possessed 10 VAGs, 17 isolates revealed 9 genes, 38 isolates revealed 8 genes, 60 isolates revealed 7 genes, 40 isolates revealed 6 genes, 28 isolates revealed 5 genes, 25 isolates revealed 4 genes, 55 isolates revealed 3 genes, 2 isolates revealed 2 genes, 4 isolates revealed one gene and 4 isolates revealed no genes.
Antibiotic susceptibility test
Standard disc diffusion method
The highest levels of antimicrobial resistance were found against sulphamethoxazole-trimethoprim, florfenicol, amoxicillin, doxycycline and spectinomycin in percentage of; 95.5, 93.7, 93.3, 92.2 and 92.2%, respectively (Table 2).
Minimal inhibitory concentration (MIC)
MIC was performed on all APEC isolates using eight different antimicrobial agents based on their common use in poultry sector. Results illustrated in Table 3 were interpreted according to animal criteria by clinical and laboratory standard institute . Escherichia coli ATCC 25922 was used as control for each run of the test.
Detection of antimicrobial resistant genes by multiplex PCR
DNA’s templates from the extraction step were used to detect the prevalence of eight antimicrobial resistance genes (ARG) among APEC isolates by multiplex PCR (Table 4). The eight antimicrobial resistance genes were present in different combinations, ranging from two genes in some isolates to eight genes in others. All isolates had at least two ARGs, 183(68%) of E. coli isolates found to possess at least 5 out of 8 ARGs, while only 3(1.1%) were found to have all the eight tested genes (Figs. 3 and 4 ).
Risk factors analysis
After excluding variables with the same answers, chi-square (X2) and fisher exact tests were performed to screen association between outcome variable (resistance status of the farm) and risk factors in univariable analysis. Twenty-nine variables included in the univariable analysis screening, only variables with P ≤ 0.25 considered for further analysis (Table 5). Nine variables have P value ≤0.25 used to perform the final logistic regression model. Collinearity between variables was tested using chi-square and spearman rank correlation test in bivariate analysis. Results of (X2) are shown in (Table 5) and results of spearman rank correlation test (Additional file 2).
The resistance status of each individual farm was used as unit of comparison, out of 84 farms completed the questionnaire; 49 (58.3%) resistance farm (presence of one or more multidrug resistant APEC isolate) coded as (1) Multidrug resistance is defined as a single bacterial isolate resistant to 3 or more antimicrobial classes (43), and 35 (41.7%) susceptible farms (no MDR-APEC present in the farm) coded as (0).
Final multivariable logistic regression
Nine variables from univariate analysis step were used to perform multivariable logistic regression model for the outcome, risk factors were considered significant when P value ≤0.05, non-significant factors re-entered when a new variable become significant or removed. Two variables with P-value ≤0.05 (water source and distance in relation to other farms) and two variables approaching significance with P-value ≤0.10 (use of antimicrobials as growth promoters and get prescription before antimicrobial treatment) remain in the final model (Table 6). The final model was tested to fit Hosmer and Lemeshow-of-fit test .
E. coli isolation
Colibacillosis is caused by APEC, which considered as one of the major threat to poultry industry and public health. In present study, APEC was isolated from broiler chickens in northern Jordan, with a clinical manifestation of colibacillosis at a prevalence rate of 53.4%. In Jordan, two investigations of broiler chickens with colibacillosis have been previously carried out with prevalence rate of 88.2 and 77% [30, 31], respectively. In other countries, the prevalence rates of colibacillosis range from 52.26 to 86.7% [32,33,34,35].
The high prevalence of E. coli infections in broiler chickens could be associated with the accumulation of E. coli aerosols in the atmosphere of chicken barns that are inhaled by chickens into the respiratory tract. Samples that gave negative bacterial culture may be collected from farms that used early antibiotic treatment policy. E. coli isolation was from chicken visceral organs which are the last stage of the disease colonization . RapID one system conformation and molecular identification were performed to reduce the false positive results.
E. coli serotypes
In the current study, serotypes O78, O1, and O2 were identified at a prevalence of 23.79, 14.86, and 12.63%, respectively. In Jordan, a study by Al-Tarazi  demonstrated that prevalent serotypes were O78 (8%), O1a and O1b (5.2%), O8 (4.8%), O127aO127b (4.8%), and O45 (4.5%) which was isolated from cases of broiler colibacillosis and egg peritonitis. However, similar results to our findings were presented in Egypt and Iran [37, 38]. In China and Northern Ireland, O78 was found as a predominant serotype in cases of broiler colibacillosis [39, 40]. It is clear that the results from this study and other previous evidence that O78, O2 and O1 were the most prevalent APEC serotypes in broiler chickens. Table 1 indicates that all serotypes are present in the four governorates included in this study which highlight that is no control measures to prevent spreading of the APEC.
Serotype O18 was identified in 1.5% of E. coli isolates that share common phenotypic and genotypic characteristics with human ExPEC and NMEC strains. This may explain the zoonotic potential of those strains . Other serotypes were isolated in less frequency, and they are of less important for poultry industry.
Virulence associated genes (VAGs)
Screening multiplex PCR for sixteen VAGs was performed for all isolates; the most prevalent genes were SitA (97.4%), iss (93.3%), iucC & D (75%), astA (71%), tsh (46.5%) and cvi (39%) genes. Presence of three out of four of iss, iucC, tsh and cvi genes indicate that the isolate is avian pathogenic E. coli  Also, Timothy  reported that presence of these genes are associated with avian colibacillosis and indicates presence of APEC. Sixty-nine percent (186 E. coli isolates) of the current study considered as pAPEC according to  report that chicken E. coli isolates carrying > 5 VAGs were classified as pAPEC. Sit A and iuc genes both contributes to iron acquisition. Sit A is usually detected in APEC more than other commensal E. coli . In this study sit A gene was detected with a high prevalence (97.3%) which is higher than the prevalence previously reported in Brazil, . High prevalence of increased serum survival protein coded by iss gene (93.3%), was higher than what was detected in USA and Germany where 80.5 and 82.7% of APEC isolated from birds with colibacillosis possess such gene [46, 47]. Tsh genes were found in 46.4% of isolates, similar to the findings of Ewers et al.  and Dozois et al.  where Tsh genes were detected at a prevalence rate of 53.3 and 49.8%, respectively. Toxin-producing genes astA was detected in 71% of the isolates which is higher than the study of  were astA detected in 21% of the tested E. coli.
In general, VAGs are integrated within the plasmid, the pathogenicity islands (chromosomally or extra chromosomally) or the bacteriophages, the acquisition of VAGs is usually through horizontal gene transfer [50, 51] which may explain the absence or the low prevalence of the remaining VAGs.
This study found lower resistance rates against beta lactams, tetracycline and fosfomycin than a previously reported . However, a higher percentage of resistance was identified in isolates against enrofloxacin, spectinomycin, gentamicin and florfenicol .
In the present study, 93.3% of the APEC isolates were resistant to amoxicillin, which is lower than the resistance rate of 100% reported in Jordan by Abu-Basha et al.,  and higher than the 83.3% resistance rate reported by Qabajah and Ashhab . In this study, 5.1% of the isolates were resistant to aztreonam, which is significantly lower than the resistance rates (41.1%) previously reported by Ahmed et al.  in Eygpt. This lower rate of resistance is likely to be due to the fact azetronem is not used in poultry in Jordan. In this study, APEC isolates were found to be resistant to doxycycline (92.2%) and oxytetracycline (55%) which is lower than the 100% resistance rate reported by Abu-Basha et al., . APEC isolates (57.2%) were found to be resistant to gentamycin, which is higher than previously reported [34, 52]. APEC isolates were found to be highly-resistant to spectinomycin (92.2%) compared to resistance rates (47%) previously reported by . APEC isolate resistance to the cephalosporin’s; ceftazidime, ceftriaxone and cefepime showed the lowest resistance levels among the tested panel of antimicrobials this result is expected for these types of cephalosporins as they are not used in poultry industry.
Attention should be paid to those antimicrobials used in broilers feed, drinking water, and as growth promoter in suboptimal doses; chlortetracycline, erythromycin, enrofloxacin, oxytetracycline and sulfonamides. The high resistance levels observed for these antibiotic classes reflect the widespread use of them in poultry. In Jordan, high frequencies of antimicrobial resistance were found in chicken isolates that can be attributed to the large-scale use of antimicrobials for disease treatment and prevention without veterinary consultation.
Antimicrobial Resistance genes
The current study targeted eight ARGs, commonly associated with antimicrobial resistance among APEC. For tetracycline resistance genes, TetA and TetB, 90.7% of the isolates expressed at least one of the tetracycline resistance genes, with TetA was the most prevalent gene. This is similar to a study carried out in Egypt, where 91.8% of APEC isolates from broilers, possessed tetracycline resistance genes, with the most prevalent type being TetB . The high prevalence of tet genes are associated with high resistant against tetracycline class (resistance range from 55 to 92.2%). A high prevalence of the class 1 integron (int1) gene was expressed by 97% of the APEC isolates, which was higher than previously reported . This finding highlighted the ability of the APEC isolates to capture ARG from other pathogenic bacteria and the environment. Sulphonamide resistance genes sul1 and sul2 were both prevalent in 70% of the APEC isolates, higher than a previous Portuguese study which found that APEC sul1 gene prevalence was 47% and sul2 was not tested . Also, the relatively high prevalence of sul1 and 2 (70%) were associated with high resistant against Sulphamethoxazole (95.5%).
Genes encoding beta-lactamases; bla-SHV and bla- TEM was identified in the APEC isolates at a prevalence of 1.8 and 72.9%, respectively. This differs from the findings of Huijbers et al.,  in the Netherlands who assessed the prevalence of ESBL producing E. coli in broiler and people living or working with broiler farms; Huijbers et al.,  study reported much higher prevalence of bla-SHV (17%) but lower bla- TEM (9.1%). The prevalence of Cat1 gene was 61.7% which is not significantly (P > 0.5) associated with high resistant to florfenicol (93.7%), this is may be due to presence of other Cat genes which are not tested.
Plasmids are considered as the main vector for horizontal gene transfer of ARGs. Increased levels of ARGs sulI, intI, aphA and traF in the aquatic environment facilitate the spread of AMR through plasmids. The high prevalence of integrons among APEC isolates (97%) which is reported to be responsible for the horizontal gene transfer and highly responsive to antimicrobial stress in the environment could explain the abundance of ARGs among the isolated APEC .
This study correlates risk factors that were hypothesized to be associated with the presence of MDR E. coli in broiler farms in Jordan. The main risk factors associated with the presence of MDR E. coli were; farms using water from artesian wells, as poultry drinking water increases the incidence of having MDR E. coli compare to farms supplied by the municipalities’ drinking water. Jordan has 12 ground water basins that serve 282 million m3 of water. This water is used for both industrial and irrigation purposes  Water environments are considered as reservoirs and amplifying sources of antimicrobial resistant genes of clinical importance .
Previous studies, performed in Canada, tested the antimicrobial resistance of Enterococcus spp. Identify that 86, 58 and 100% of the isolates were resistant to more than one type of antibiotic in poultry litter, surface water and ground water isolates, respectively . This finding suggests that there is a high presence of antibiotic resistant genes in surface water, wastewater, and poultry litter.
Furthermore, this study found that farms located in close proximity to other poultry farms were at high risk of contamination with MDR E. coli which is similar to finding of Hartung & Schukz , emphasized that serious pathogens are transmitted by air, which is positively correlated to farm density, considering farmers have no control over farm location. Therefore, farmers should pay attention toward wind directions in their area. Personal movement, vehicles and instruments can also be considered as vectors for transmission of pathogens.
Other potential risk factors related to antimicrobial usage were the use of antimicrobial agents as growth promoters and the administration of antibiotic without veterinary consultation. Many studies support that the improper use of antibiotics for increasing productivity, enhances the selection pressure for antimicrobial resistant pathogens [7, 61].
Public health concerns regarding antimicrobial residues and antimicrobial resistance pathogens in food and the environment reinforce the need for more research on safer alternatives to antibiotics as feed additives . Netherlands was ranked as the highest antimicrobial consuming country in 2007, with an estimated 600 tons of therapeutic antimicrobials used in the veterinary sector. Therefore, the Netherlands set up a monitoring action plan to reduce the antimicrobial use in animals. The first step taken was to establish a veterinary medicine authority, whose main purpose was to record antimicrobial usage and prescription from farmers and Veterinarians, and to set species-specific annual targets for antimicrobial use. This action plan resulted in a 56% reduction in antimicrobial usage in the period between 2007 and 2012 .
This study characterised the VAGs of avian pathogenic E. coli and establish their antimicrobial resistance patterns. The widespread of antimicrobial resistance of APEC isolates and detection of ARGs highlighted the need to monitor the spread of ARGs in poultry farms and the environment in Jordan. Use of ground water and closely located farms were significant risk factors associated with the presence of MDR APEC in broiler chickens in Jordan.
Chicken samples were collected from farms located in northern Jordan; Irbid, Jerash, Ajlune, and Mafraq governorates, which contain 896 broiler farms with annual capacity 12, 064,600 bird .
Sample size determination
According to the sample size formula from an infinite population:
Where; p = estimated prevalence of disease in the population, q = (1-p), d = accepted margin of error and Z the value for specific confidence level.
The confidence level is 95%, Z value = (1.96), Estimate prevalence = 88.2% , d = (0.05) thus, n = 159.8 farms.
Eighty-four farms were visited and asked to fill in the questionnaire before samples collection. Five hundred and four sick birds’ samples were collected during the period from April to December 2016.
A questionnaire was designed with 42 questions divided into four sections, which covered the factors believed to be associated with antimicrobial resistance. The questionnaire was translated to Arabic and answered by the owners or the veterinarian of each farm during personal interviews while collecting the samples. The questionnaire was field pre-validated. (Additional file 1).
Isolation and conventional identification of APEC
Aseptic swabs from liver, heart, spleen and lungs of birds symptomatic of colibacillosis were cultured on 5% sheep blood agar and on MacConkey agar media (Oxoid), and subcultured on selective differential media eosin methylene blue agar (EMB) (Oxoid) . The isolated bacteria were identified as E. coli by observing their cultural characteristics, morphology by Gram’s stain, oxidase test, biochemical reactions using indole, methyl-red, Voges-Preuskuar and citrate tests (IMViC), Kligler Iron Agar (KIA) and motility test as described by Tonu et al. . The suspected isolates were maintained in cryostat tubes containing 20% glycerol with LB Luria Bertani broth at − 70 °C .
Confirmation of APEC using RapID™ ONE system
E. coli isolates were tested using RapID ONE system Kit (Remel, USA) as indicated in the kit catalogue, and results then were interpreted using ERIC (Remel RapID database).
Serotyping was conducted using E. coli polyvalent O antisera and mono-specific antisera prevalent in poultry; O1, O2, O78, O8, O9, O18, O26, O25, O45, O55, O86, O111, O114, O119, O127, and O128 [30, 31, 67]. All the E. coli isolates were subjected to serotyping according to the instructions of the manufacturer (SSI Diagnostica) using a micro titre plate agglutination test.
Molecular identification of APEC
DNA extraction and detection of 16 s rRNA gene of E. coli by PCR
Extraction of DNA from the Escherichia coli was carried out by boiling procedure and rapid cooling method. In brief, a single colony of E coli was resuspended in 100 μl of nuclease free water and boiled for 10 min and immediately cooled on icebox followed by centrifugation at 10,000 rpm for 10 min. The supernatant was collected, stored at − 20 °C and used as DNA template .
E. coli isolates were confirmed by detection of 16S rRNA gene using conventional PCR. As described by Hossain et al., . Oligonucleotide primers sequences used for the amplification of 16S rRNA gene of E. coli was 16 s-F: GAC CTC GGT TTA GTT CAC AGA and 16 s-R: CAC ACG CTG ACG CTG ACC A, location within the gene 4,267,278–4,267,845 and amplicon size 485 bp. PCR reaction mixture consisted of 12.5 μl of 2 × PCR master mixtures (Promega), 10 pmol primer of each and 2 μl of genomic DNA in a final volume of 25 μl adjusted by nuclease free water. The cycling conditions consisted of initial denaturation at 95 °C for 5 min., followed by 30 cycles of 94 °C for 1 min., 55 °C for 45 s min. and 72 °C for 1 min., with final extension at 72 °C for 7 min. The amplified products were electrophoresed into 1.8% agarose gel at 100 V visualized under Gel doc/UV trans-illuminator.
Multiplex polymerase chain reaction method for detection of virulence associated genes (VAGs)
Each DNA extract was screened for 16 VAGs associated with avian pathogenic E. coli; sfa, iss, tsh, kps, kpsII, kpsIII, iucC, iucD, hlyD, ibeA, sitA, astA, cvi, papC, irp2 and vat, using a multiplex PCR . Primers were obtained from GENEWIZ Company (USA) and Intron, South Korea supplied all PCR constituents used in this study. All sixteen primers sequences were given in . Briefly, each 50 μl PCR reaction contained: 12 μl of 25 mM MgCl2, 21.3 μl nuclease free water, 5 μl 10x PCR buffer, 4 μl of 20 mM dNTPs, 0.3 μl of each 100 pmol forward and reverse primer, 0.3 μl, 5 U/ μl Taq polymerase and 5 μl template DNA. Thermocycler conditions were: initial denaturation 95 °C for 5 min; nine cycles of 95 °C for 60 s, 55 °C for 30 s, 72 °C for 60 s; twenty-eight cycles of 94 °C for 30 s, 55 °C for 30 s, 72 °C for 30 s with a final extension 72 °C for 7 min. The mixture was held at 4 °C. PCR products were subject to electrophoresis on a 2% agarose gel in tris–acetate buffer (TAE) at 150 V for 60 min alongside a super Ladder-Low 100 bp ladder (Intron, South Korea).
Two separate m-PCR assays were performed; one multiplex PCR previously described by Ewers et al.  and one m-PCR assays for ibeA and sitA described by Timothy et al. . Briefly, for a 25 ml multiplex PCR, 4 μl of 25 mM MgCl2, 13.9 μl nuclease free water, 2.5 μl 10x PCR buffer, 0.5 μl 20 mM dNTPs, 0.1 μl of each 100 pmol forward and reverse primers, 0.5 μl 5 U/ μl Taq polymerase and 2 μl DNA templates were used. Multiplex PCR thermocycler conditions were as follows: initial denaturation 94 °C for 3 mints followed by 25 cycles of: 94 °C for 30 s, 58 °C for 30 s, 68 °C for 3 mints with a final extension 72 °C for 10 mints. The mixture was held at 4 °C. Each individual PCR contained 1 μl DNA template, 1 μl of each primer (100 pmol) and 22 μl of 1.1x Reddymix PCR master mix with 1.5 mM MgCl2. M-PCR thermocycler conditions for sitA and ibeA were; 95 °C for 12 min and 25 cycles of: 94 °C for 30 s, 63 °C for 30 s, 68 °C for 3 min; 72 °C for 10 min with a final hold 4 °C. PCR products were subject to electrophoresis as above. Isolates carrying > 5 VAGs were classified as APEC.
Standard disc diffusion method
The agar disk diffusion test was carried out according to . All E. coli isolates were tested for 19 antibiotics: amoxicillin (25 μg), doxycycline (30 μg), ciprofloxacin (5 μg), ceftriaxone (30 μg), gentamicin (10 μg), florfenicol (30 μg), cefepime (30 μg), aztreonam (30 μg), imipenem (10 μg), cephalexin (30 μg), ceftazidime (30 μg), sulphamethoxazole-trimethoprim (23.75/1.25 μg), Amoxicillin-clavulanate (20/10 μg), apramycin (15 μg), spectinomycin (25 μg), Enrofloxacin (5 μg), Oxytetracycline (30 μg), Chlortetracycline (10 μg), and Fosfomycin (50 μg). Escherichia coli ATCC 25922 was used as control strain.
Minimal inhibitory concentration (MIC)
Susceptibility to 8 antimicrobials was evaluated by broth microdilution  Cationic-adjusted Muller-Hinton broth (Cationic-adjusted Muller-Hinton, Fluka, Switzerland) was used to prepare the bacterial inoculum and dilute the antimicrobial agents (Table 7). According to the MIC breakpoints, E. coli isolates that were resistant to 3 or more antimicrobial classes were considered multidrug-resistant isolates . The reference Escherichia coli ATCC 25922 strain was used as a control strain.
Molecular detection of antimicrobial resistant genes by multiplex PCR
PCR was conducted for the E. coli isolates that were found resistant to one or more of the previously mentioned antimicrobials, as described by . The DNA templates from the DNA extraction step were used to detect resistance genes(Table 8).
Eighty-four broiler farms completed the questionnaire and were included in the analysis using SPSS 21.0 software. Questions with the same answers were excluded from the analysis (application of “all in all out” strategy, disinfection of farm building before introduction of new flocks, application of vaccination program, previous history of respiratory diseases, monitoring of mortality rate and use of antimicrobials for disease treatment).
Chi-square (X2) and Fisher exact tests were performed to screen association between outcome variable (resistance status of the farm) and risk factors in univariable analysis. Only variables with P ≤ 0.25 considered for further analysis, which were used to perform the final logistic regression model. Collinearity between variables was tested using chi-square and Spearman rank correlation test in bivariate analysis.
The resistance status of a farm was used as unit of comparison, farms were categorized into resistance according to presence of one or more multidrug resistant APEC isolate coded as (1) and susceptible isolates coded as (0) depending on the multidrug resistance definition. According to WHO  five antimicrobial agents (OT, CN, CIP, AML and FOS) were selected in order to categorize the isolates into multidrug resistant patterns (resistant to three or more antimicrobials) and sensitive isolates .
Final multivariable logistic regression
Variables from univariate analysis step were used to perform multivariable logistic regression model for the outcome, risk factors were considered significant when P value ≤0.05, non-significant factors re-entered when a new variable become significant or removed. The final model was tested to fit hosmer and lemeshow-of-fit test.
Avian Pathogenic E. coli
Antimicrobial Resistance Genes
Chronic Respiratory Disease
Eosin Methylene Blue Agar
Extraintestinal Pathogenic E. coli
Kligler Iron Agar
Luria Bertani Broth
Minimum Inhibitory Concentration
Multiplex Polymerase Chain Reaction
Neonatal Meningitis E. coli
Polymerase Chain Reaction
Quaternary ammonium compounds
Rounds per Minuet
Ribosomal Ribonucleic Acid
United States of America
Virulence associated genes
Alexander DJ, Senne DA. Newcastle disease, other avian paramyxoviruses, and Pneumovirus infections. Disseases of poultry; 2008.
Lutful Kabir SM. Avian colibacillosis and salmonellosis: a closer look at epidemiology, pathogenesis, diagnosis, control and public health concerns. Int J Environ Res Public Health. 2010;7(1):89–114. https://doi.org/10.3390/ijerph7010089.
Dziva F, Stevens MP. Colibacillosis in poultry: unravelling the molecular basis of virulence of avian pathogenic Escherichia coli in their natural hosts. Avian Pathology : Journal of the WVPA. 2008;37(4):355–66.
Mellata M, Dho-Moulin M, Dozois CM, Curtiss R, 3rd Brown PK, Arné P, … Fairbrother JM. Role of Virulence Factors in Resistance of Avian Pathogenic Escherichia coli to Serum and in Pathogenicity. Infect Immun. 2003;71(1):536–40. https://doi.org/10.1128/iai.71.1.536-540.2003.
Collingwood C, et al. Is the Concept of Avian Pathogenic Escherichia coli as a Single Pathotype Fundamentally Flawed? Front Vet Sci. 2014;1:1–4 Available at: https://www.frontiersin.org/articles/10.3389/fvets.2014.00005/full.
Schwarz S, Shen J, Kadlec K, Wang Y, Michael GB, Feßler AT, Vester B. Lincosamides, Streptogramins, Phenicols, and Pleuromutilins: Mode of Action and Mechanisms of Resistance. Cold Spring Harb Perspect Med. 2016;6(11):1–30.
Levy S. Reduced antibiotic use in livestock: how Denmark tackled resistance. Environ Health Perspect. 2014;122(6):160–5. https://doi.org/10.1289/ehp.122-A160.
Dargatz DA, Fedorka-Cray PJ, Ladely SR, Ferris KE. Survey of Salmonella serotypes shed in feces of beef cows and their antimicrobial susceptibility patterns. J Food Prot. 2000;63(12):1648–53.
Landoni MF, Albarellos G. The use of antimicrobial agents in broiler chickens. Vet J. 2015a;205(1):21–7. https://doi.org/10.1016/j.tvjl.2015.04.016.
Zakeri A, Kashefi P. Antimicrobial susceptibilities of avian Escherichia coli isolates in Tabriz, Iran. Afr J Biotechnol. 2012;11(19):4467–70. https://doi.org/10.5897/AJB11.3168.
Awad A, Arafat N, Elhadidy M. Genetic elements associated with antimicrobial resistance among avian pathogenic Escherichia coli. Ann Clin Microbiol Antimicrob. 2016;15(1):59. https://doi.org/10.1186/s12941-016-0174-9.
Shin SW, Shin MK, Jung M, Belaynehe KM, Yoo HS. Prevalence of antimicrobial resistance and transfer of tetracycline resistance genes in Escherichia coli isolates from beef cattle. Appl Environ Microbiol. 2015;81(16):5560–6. https://doi.org/10.1128/AEM.01511-15.
Adefisoye MA, Okoh AI. Identification and antimicrobial resistance prevalence of pathogenic Escherichia coli strains from treated wastewater effluents in eastern cape, South Africa. MicrobiologyOpen. 2016;5(1):143–51. https://doi.org/10.1002/mbo3.319.
Tamang MD, Oh JY, Seol SY, Kang HY, Lee JC, Lee YC, et al. Emergence of multidrug-resistant Salmonella enterica serovar Typhi associated with a class 1 integron carrying the dfrA7 gene cassette in Nepal. Int J Antimicrob Agents. 2007;30(4):330–5. https://doi.org/10.1016/j.ijantimicag.2007.05.009.
Sharma VK, Johnson N, Cizmas L, McDonald TJ, Kim H. A review of the influence of treatment strategies on antibiotic resistant bacteria and antibiotic resistance genes. Chemosphere. 2016;150:702–14. https://doi.org/10.1016/j.chemosphere.2015.12.084.
Akond M. Antibiotic Resistance of Escherichia Coli isolated from poultry and poultry environment of Bangladesh. Am J Environ Sci. 2009;5(1):47–52. https://doi.org/10.3844/ajessp.2009.47.52.
Snary EL, Kelly LA, Davison HC, Teale CJ, Wooldridge M. Antimicrobial resistance : a microbial risk assessment perspective. Journal of Antimicrobial Chemotherapy, 2004;53(6):906–17. https://doi.org/10.1093/jac/dkh182.
Castanon JIR. History of the use of antibiotic as growth promoters in European poultry feeds. Poult Sci. 2007;86(11):2466–71. https://doi.org/10.3382/ps.2007-00249.
Diarra MS, Malouin F. Antibiotics in Canadian poultry productions and anticipated alternatives. Front Microbiol. 2014;5:1–15. https://doi.org/10.3389/fmicb.2014.00282.
Apata DF. Antibiotic resistance in poultry. Inter J Poul Sci. 2009;8(4):404–8.
McEwen SA, Fedorka-Cray PJ. Antimicrobial use and resistance in animals. Clin Infect Dis. 2002;34(Suppl 3):S93–s106. https://doi.org/10.1086/340246.
Fletcher S, Fletcher S. Understanding the contribution of environmental factors in the spread of antimicrobial resistance. Environ Health Prev Med. 2015. https://doi.org/10.1007/s12199-015-0468-0.
Dancer SJ. How antibiotics can make us sick: the less obvious adverse effects of antimicrobial chemotherapy. Lancet Infect Dis. 2004;4(10):611–9. https://doi.org/10.1016/S1473-3099(04)01145-4.
Nhung N, Thuy C, Trung N, Campbell J, Baker S, Thwaites G, Carrique-Mas J. Induction of Antimicrobial Resistance in Escherichia coli and Non-Typhoidal Salmonella Strains after Adaptation to Disinfectant Commonly Used on Farms in Vietnam. Antibiotics. 2015;4(4):480–94. https://doi.org/10.3390/antibiotics4040480.
Buffet-Bataillon S, Tattevin P, Bonnaure-Mallet M, Jolivet-Gougeon A. Emergence of resistance to antibacterial agents: the role of quaternary ammonium compounds - a critical review. Int J Antimicrob Agents. 2012;39(5):381–9. https://doi.org/10.1016/j.ijantimicag.2012.01.011.
Chen Z, Jiang X. Microbiological safety of chicken litter or chicken litter-based organic fertilizers: a review. Agriculture. 2014;4(1):1–29. https://doi.org/10.3390/agriculture4010001.
Greig J, Rajić A, Young I, Mascarenhas M, Waddell L, Lejeune J. A scoping review of the role of wildlife in the transmission of bacterial pathogens and antimicrobial resistance to the food chain. Zoonoses Public Health. 2015;62(4):269–84. https://doi.org/10.1111/zph.12147.
Clinical and Laboratory Standards Institute, 2012. Supplement M100-S22, Vol.32 No. 3.
Tsai M-H, Chu S-M, Hsu J-F, Lien R, Huang H-R, Chiang M-C, et al. Risk factors and outcomes for multidrug-resistant gram-negative bacteremia in the NICU. Pediatrics. 2014;133(2):e322–9. https://doi.org/10.1542/peds.2013-1248.
El-Sukhon SN, Musa A, Al-Attar M. Studies on the bacterial etiology of airsaculitis of broilers in northern and middle Jordan with special reference to Escherichia coli, Ornithobacterium rhinotracheale, and Bordetella avium. Avian Dis. 2002;46(3):605–12.
Al-Tarazi YH, Taqi-Eddin AR. Master thesis, entitled: Biochemical, Serological and Pathogenicity Studies on Escherichia coli isolated from chicken in Jordan. Amman, Jordan: Department of Biological Sciences, Faculty of Science, Jordan University; 1983.
Halfaoui Z, Menoueri NM, Bendali LM. Serogrouping and antibiotic resistance of Escherichia coli isolated from broiler chicken with colibacillosis in center of Algeria. Vet World. 2017;10(7):830–5. https://doi.org/10.14202/vetworld.2017.830-835.
Sen A, Veterinary C, Sciences A, Kanti S, Chittagong N, Sciences A, Sen A. Isolation of Escherichia coli from the liver and yolk sac of day old chicks with their Antibiogram isolation of Escherichia coli from the liver and yolk sac of day old chicks with their Antibiogram, (January); 2017.
Ahmed AM, Shimamoto T, Shimamoto T. Molecular characterization of multidrug-resistant avian pathogenic Escherichia coli isolated from septicemic broilers. Int J Med Microbiol. 2013;303(8):475–83. https://doi.org/10.1016/j.ijmm.2013.06.009.
Ahmed MS, Sarker A, Rahman MM. Prevalence of infectious diseases of broiler chickens in Gazipur District, vol. 7; 2009. p. 326–31.
Filho HCK, Brito KCT, Cavalli LS, Brito BG. Avian pathogenic Escherichia coli ( APEC ) - an update on the control. In: The Battle against Microbial Pathogens: Basic Science, Technological Advances and Educational Programs; 2015. p. 598–618.
Younis G, Awad A, Mohamed N. Phenotypic and genotypic characterization of antimicrobial susceptibility of avian pathogenic Escherichia coli isolated from broiler chickens. Vet World. 2017;10(10):1167–72. https://doi.org/10.14202/vetworld.2017.1167-1172.
SEIFI S, KHOSHBAKHT R, ATABAK AR. Antibiotic susceptibility, serotyping and pathogenicity EVALUATION of avian ESCHERICHIA COLI isolated from broilers in northern Iran. Bulg J Vet Med. 2015;18(1):74–82. https://doi.org/10.15547/bjvm.8.
Dou X, Gong J, Han X, Xu M, Shen H, Zhang D, Zou J. Characterization of avian pathogenic Escherichia coli isolated in eastern China. Gene. 2016;576(1):244–8. https://doi.org/10.1016/j.gene.2015.10.012.
McPeake SJW, Smyth JA, Ball HJ. Characterisation of avian pathogenic Escherichia coli (APEC) associated with colisepticaemia compared to faecal isolates from healthy birds. Vet Microbiol. 2005;110(3–4):245–53. https://doi.org/10.1016/j.vetmic.2005.08.001.
Krishnan S, Chang AC, Hodges J, Couraud PO, Romero IA, Weksler B, Nicholson BA, Nolan LK, Prasadarao NV. Serotype O18 avian pathogenic and neonatal meningitis Escherichia coli strains employ similar pathogenic strategies for the onset of meningitis. Virulence. 2015;6(8):777–86.
Schouler C, Schaeffer B, Brée A, Mora A, Dahbi G, Biet F, Moulin-Schouleur M. Diagnostic strategy for identifying avian pathogenic Escherichia coli based on four patterns of virulence genes. J Clin Microbiol. 2012;50(5):1673–8. https://doi.org/10.1128/JCM.05057-11.
Timothy S, Shafi K, Leatherbarrow AH, Jordan FTW, Wigley P. Molecular epidemiology of a reproductive tract associated colibacillosis outbreak in a layer breeder flock associated with atypical avian pathogenic Escherichia coli. Avian Pathol. 2008;37(4):375–8. https://doi.org/10.1080/03079450802216579.
Kemmett K, Williams NJ, Chaloner G, Humphrey S, Wigley P, Humphrey T. The contribution of systemic Escherichia coli infection to the early mortalities of commercial broiler chickens. Avian Pathol. 2014;43:37–42.
Andrew JM, BSAC working party on susceptibility testing. BSAC standardized disc susceptibility testing method (version 7). J Antimicrob Chemother. 2008;62:256–78. https://doi.org/10.1093/jac/dkn194.
Dissanayake DR, Octavia S, Lan R. Population structure and virulence content of avian pathogenic Escherichia coli isolated from outbreaks in Sri Lanka. Vet Microbiol. 2014;31(168):403–12.
Ewers C, Janßen T, Kießling S, Philipp H, Wieler LH. Molecular epidemiology of avian pathogenic Escherichia coli (APEC) isolated from colisepticemia in poultry. Vet Microbiol. 2004;104(1–2):91–101.
Dozois CM, Dho-Moulin M, Bree A, Fairbrother JM, Desautels C, Curtiss R 3rd. Relationship between the Tsh autotransporter and pathogenicity of avian Escherichia coli and localization and analysis of the Tsh genetic region. Infect Immun. 2000;68:4145–54.
Skyberg JA, Horne SM, Giddings CW, Wooley RE, Gibbs PS, Nolan LK. 2003. Characterizing avian Escherichia coli isolates with multiplex polymerase chain reaction. Avian Dis. 2003;47(4):1441–7.
Hacker J, Kaper JB. Pathogenicity island and the evolution of microbes. Annu Rev Microbiol. 2000;54:641–79.
Ochman H, Lawrence JG, Groisman EA. Lateral gene transfer and the nature of bacterial innovation. Nature. 2000;405:299–304.
Abu-Basha EA, Gharaibeh SM, Thabet AM. In vitro susceptibility of resistant Escherichia coli field isolates to antimicrobial combinations; 2012.
Qabajah M, Ashhab Y. Avian pathogenic Escherichia coli (APEC) in Palestine. Vet Res. 2012;30(2–3):299–316 Available at: https://biotech.ppu.edu/sites/default/files/publications/Avian%20Pathogenic%20Escherichia%20coli.pdf.
Mendonça N, Mendonça N, Figueiredo R, Mendes C, Card RM, Anjum MF, da Silva GJ. Microarray Evaluation of Antimicrobial Resistance and Virulence of Escherichia coli Isolates from Portuguese Poultry. Antibiotics. 2016;5(1):4 Available at: http://www.mdpi.com/2079-6382/5/1/4/htm.
Huijbers PMC, Graat EA, Haenen AP, van Santen MG, van Essen-Zandbergen A, Mevius DJ, van Duijkeren E, van Hoek AH. Extended-spectrum and AmpC β-lactamase-producing Escherichia coli in broilers and people living and/or working on broiler farms: prevalence, risk factors and molecular characteristics. J Antimicrob Chemother. 2014;69(10):2669–75.
Wang Q, Lu Q, Mao D, Cui Y, Luo Y. The horizontal transfer of antibiotic resistance genes is enhanced by ionic liquid with different structure of varying alkyl chain length. Front Microbiol. 2015;6:1–8.
Ministry of agriculture. Animal production department annual report, 2014.
Pruden A, Joakim Larsson DG, Amézquita A, Collignon P, Brandt KK, Graham DW, Lazorchak JM, Suzuki S, Silley P, Snape JR, Topp E, Zhang T, Zhu Y-G. Management options for reducing the release of antibiotics and antibiotic resistance genes to the environment. Environ Health Perspect. 2013;121(8):878–85.
Furtula V, Jackson CR, Farrell EG, Barrett JB, Hiott LM, Chambers PA. Antimicrobial resistance in Enterococcus spp. isolated from environmental samples in an area of intensive poultry production. Int J Environ Res Public Health. 2013;10(3):1020–36.
Hartung J, Schulz J. Risks caused by bio-aerosols in poultry houses. International Conference: Poultry in the 21st century, avian influenza and beyond. Bangkok; 2007.
Van Boeckel TP, Brower C, Gilbert M, Grenfell BT, Levin SA, Robinson TP, Teillant A, Laxminarayan R. Global trends in antimicrobial use in food animals. Proc Natl Acad Sci. 2015;112(18):5649–54. Available at: https://www.pnas.org/content/pnas/112/18/5649.full.pdf.
Speksnijder DC, Mevius DJ, Bruschke CJ, Wagenaar JA. Reduction of veterinary antimicrobial use in the Netherlands. The dutch success model. Zoonoses Public Health. 2015;62(s1):79–87.
Ministry of Water & Irrigation. 2013. Jordan Water Sector Facts and Figures 2013 Opening Statement. Available at: http://www.mwi.gov.jo/sites/en-us/Documents/W.inFig.EFINALE.pdf.
Hossain MK, Rahman M, Nahar A, Khair A, Alam MM. Isolation and identification of diarrheagenic Escherichia coli causing colibacillosis in calf in selective areas of Bangladesh. Bangl J Vet Med. 2013;11(2):145–9.
Giovanardi D, Campagnari E, Ruffoni LS, Pesente P, Ortali G, Furlattini V. Avian pathogenic Escherichia coli transmission from broiler breeders to their progeny in an integrated poultry production chain. Avian Pathol. 2005;34(March 2016):313–8.
Tonu MS, Sufian S, Sarker, Kamal MM, Rahman MH, Hossain MM. Pathological study on Colibacillosis in chickens and detection of Escherichia coli by Pcr. Bangladesh Journal of Veterinarian Medicine. 2011;9(1):17–25.
Knöbl T, Moreno AM, Paixão R, Gomes TAT, Vieira MAM, da Silva Leite D, BlancoJ E, Ferreira AJP. Prevalence of avian pathogenic Escherichia coli (APEC) clone harboring sfa gene in Brazil. Sci World J. 2012:1–7. Available at: http://www.hindawi.com/journals/tswj/2012/437342/.
Dolejska M, Cizek A, Literak I. High prevalence of antimicrobial-resistant genes and integrons in Escherichia coli isolates from black-headed gulls in the Czech Republic. J Appl Microbiol. 2007;103(1):11–9.
W H O, Surveillance, I., & Resistance, A. (2016). Critically important antimicrobials for human medicine.
Borges CA, Maluta RP, Beraldo LG, Cardozo MV, Guastalli EAL, Kariyawasam S, et al. Captive and free-living urban pigeons (Columba livia) from Brazil as carriers of multidrug-resistant pathogenic Escherichia coli. Vet J. 2017;219:65–7. https://doi.org/10.1016/j.tvjl.2016.12.015.
The authors would like to thank Deanship of research at JUST and the OIE for funding the project, Dr. Ana Mateus, Royal Veterinary College, London, for her epidemiological expertise and statistical advice.
The project was approved and funded in 31. Aug 2016 by Jordan University of science and technology, research council under reference number 267/2016 and the World Organization for Animal Health (OIE), where funders have no impact on the study design or results.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
This study was approved by the animal care and use committee at Jordan University of science and technology, with reference number 16/3/3/797 and approval date 27. Nov.2017. Animals involved in this study were euthanized by cervical dislocation as approved by the American veterinary medical association the process of euthanasia were performed as high ethical and welfare standards as possible.
Consent for publication
Author Tillie L. Cryer previously worked for BMC veterinary research as an assistant editor. Tillie L. Cryer is now currently the editor of the BMC veterinary research. At the time of submission and peer review of this manuscript Tillie L. Cryer was not employed by BMC. Author Tillie L. Cryer has not played any editorial role in the handling of this manuscript and this manuscript has been independently peer reviewed.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.