Virulence Determinants and Plasmid-Mediated Colistin Resistance mcr Genes in Gram-Negative Bacteria Isolated From Bovine Milk

A major increase of bacterial resistance to colistin, a last-resort treatment for severe infections, was observed globally. Using colistin in livestock rearing is believed to be the ground of mobilized colistin resistance (mcr) gene circulation and is of crucial concern to public health. This study aimed to determine the frequency and virulence characteristics of colistin-resistant Gram-negative bacteria from the milk of mastitic cows and raw unpasteurized milk in Egypt. One hundred and seventeen strains belonging to Enterobacteriaceae (n = 90), Pseudomonas aeruginosa (n = 10), and Aeromonas hydrophila (n = 17) were screened for colistin resistance by antimicrobial susceptibility testing. The genetic characteristics of colistin-resistant strains were investigated for mcr-1–9 genes, phylogenetic groups, and virulence genes. Moreover, we evaluated four commonly used biocides in dairy farms for teat disinfection toward colistin-resistant strains. Multidrug-resistant (MDR) and extensive drug-resistant (XDR) phenotypes were detected in 82.91% (97/117) and 3.42% (4/117) of the isolates, respectively. Of the 117 tested isolates, 61 (52.14%) were colistin resistant (MIC >2 mg/L), distributed as 24/70 (34.29%) from clinical mastitis, 10/11 (90.91%) from subclinical mastitis, and 27/36 (75%) from raw milk. Of these 61 colistin-resistant isolates, 47 (19 from clinical mastitis, 8 from subclinical mastitis, and 20 from raw milk) harbored plasmid-borne mcr genes. The mcr-1 gene was identified in 31.91%, mcr-2 in 29.79%, mcr-3 in 34.04%, and each of mcr-4 and mcr-7 in 2.13% of the colistin-resistant isolates. Among these isolates, 42.55% (20/47) were E. coli, 21.28% (10/47) A. hydrophila, 19.12% (9/47) K. pneumoniae, and 17.02% (8/47) P. aeruginosa. This is the first report of mcr-3 and mcr-7 in P. aeruginosa. Conjugation experiments using the broth-mating technique showed successful transfer of colistin resistance to E. coli J53-recipient strain. Different combinations of virulence genes were observed among colistin-resistant isolates with almost all isolates harboring genes. Hydrogen peroxide has the best efficiency against all bacterial isolates even at a low concentration (10%). In conclusion, the dissemination of mobile colistin resistance mcr gene and its variants between MDR- and XDR-virulent Gram-negative isolates from dairy cattle confirms the spread of mcr genes at all levels; animals, humans, and environmental, and heralds the penetration of the last-resort antimicrobial against MDR bacteria. Consequently, a decision to ban colistin in food animals is urgently required to fight XDR and MDR bacteria.


INTRODUCTION
Escherichia coli, Klebsiella pneumoniae, and Pseudomonas a e r u g i n o s a a r e th e m o s t c o m m o n G r a m -n e g a t i v e environmental pathogens in dairy farms. They could infect the mammary glands of dairy cows through the teat end and colonize the mammary tissue, causing mastitis (Scatamburlo et al., 2015;Gao et al., 2017). Although Aeromonas species are autochthonous to aquatic environments, Aeromonas hydrophila has been reported as a common contaminant in a diverse variety of foods such as milk and milk products (Sharma and Kumar, 2011). Contamination of milk with pathogenic Gram-negative bacteria is often derived directly from the udder excretion of infected animals or the farm environment, thus putting the consumers who consume unpasteurized milk at a high risk for developing foodborne diarrheal diseases (Garcia et al., 2019).
In Egypt, there are no regulations that control the use of antibiotics in animal husbandry for growth promotion or the prevention and treatment of bacterial diseases (WHO, 2013). The misuse of antibiotics in dairy production is a matter of concern for the effective antibiotic therapy of infected cows, food safety, and occupational exposure, and a potential threat to public health due to the selection of multidrug-resistant (MDR) bacteria and entry of antibiotic residues into bulk tank milk (Oliveira and Ruegg, 2014;Locatelli et al., 2019). The emergence of MDR and extensive drug-resistant (XDR) Gramnegative bacteria in livestock production increased the interest in the use of colistin (polymyxin E) as a last-resort antibiotic for the treatment of these pathogens (Kempf et al., 2016). Although colistin was prohibited in developed countries, it is still used in animal husbandry in Egyptian dairy farms (Lima Barbieri et al., 2017;Dandachi et al., 2019). Hence, the extensive use of colistin increases colistin-resistant Gram-negative bacteria. Also, it imposes a major public health concern (Lima Barbieri et al., 2017) due to the dissemination of colistin resistance genes (mcr) through mobilized plasmids among animal strains and subsequently transmitted to humans through the food chain or direct contact (McEachran et al., 2015;Garcıá et al., 2018;Hmede and Kassem, 2019). Nine plasmid-borne mcr-family genes (mcr-1-mcr-9) have been reported in over 40 countries from different continents across the globe; most of them were detected in several enterobacteria, including E. coli and K. pneumoniae (Luo et al., 2020). The mcr-1 gene has been detected in most continents, while other genes have only been found in a few countries (Ling et al., 2020). The mcr-1 gene has been reported previously in colistin-resistant E. coli isolated from raw milk (Coton et al., 2012;Hassen et al., 2019), mastitis bovine milk (Filioussis et al., 2019;Liu G. et al., 2020), and cheese (Coton et al., 2012;Hammad et al., 2019). A higher frequency of mcr-1 compared with other mcr genes has been reported in colistinresistant E. coli (98.9%), K. pneumoniae (100%), and P. aeruginosa (100%) recovered from animal and human sources (Javed et al., 2020). The implementation of pre-and post-milking teat disinfection is a critical effective means of reducing the incidence of clinical and subclinical mastitis as well as new intramammary infections caused by Gram-negative environmental pathogens in dairy farms such as E. coli, K. pneumoniae, and P. aeruginosa (Tiwari, 2013;Kamal and Bayoumi, 2015;Böhm et al., 2017;Rowe et al., 2018). Choosing the ideal disinfectant with an accurate concentration is vital to control resistant bacteria in dairy farms. The resistant bacteria can survive at concentrations many folds below the minimum inhibitory concentration (MIC) of the biocidal agents (Hughes and Andersson, 2012). Notably, Gramnegative bacteria have a risk for promoting antibiotic resistance after exposure to sublethal concentrations of some disinfectants such as chlorhexidine (Kampf, 2018).

Gram-Negative Bacterial Strains
A total of 117 Gram-negative bacterial strains belonging to Enterobacteriaceae (n = 90), P. aeruginosa (n = 10), and A. hydrophila (n = 17) were included in this study. The bacterial isolates were recovered previously from milk samples of dairy cows showing mastitis as well as bulk tank raw milk during 2018-2020. Gram-negative isolates were selected based on the provided data from the Microbiology Department, Faculty of Veterinary Medicine, Zagazig University, Egypt. Enterobacteriaceae strains included non-O 157 E. coli (n = 30), O 157 E. coli (n = 12), K. pneumoniae (n = 18), Enterobacter species (n = 10), and Citrobacter species (n = 20). The isolated bacteria were sourced from different farms. All the isolates were independent and were not clonally duplicated from the same source. The presumptive bacterial isolates were confirmed adopting the standard microbiological procedures (Quinn et al., 1994). In brief, MacConkey's, eosin methylene blue (Oxoid, Cambridge, UK) and HiCrome klebsiella-selective agar media (Himedia, Mumbai, India) were used for the cultivation of Enterobacteriaceae strains. Aeromonas species were grown on Rimler-Shotts medium (HiMedia, India) with a novobiocin (Oxoid, UK) supplement, while P. aeruginosa strains were cultivated on Pseudomonas cetrimide agar (Oxoid, UK). Serotyping of E. coli was applied using diagnostic polyvalent and monovalent O and H antisera (Denka Seiken Co., Tokyo, Japan). Genomic DNA was extracted from fresh bacterial cultures using the QIAamp DNA Mini kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. Polymerase chain reaction (PCR)-based confirmation of Gram-negative bacteria was applied using genus-and speciesspecific primer sets depicted in Supplementary Table S1.
The multiple antibiotic resistance (MAR) indices were evaluated as documented elsewhere (Tambekar et al., 2006), while the MDR and XDR phenotypes were reported according to Magiorakos et al. (2012). E. coli ATCC 25922, P. aeruginosa ATCC 27853, and A. hydrophila ATCC 7966 were used as quality control strains.

Detection of Plasmid-Mediated Colistin Resistance Genes
Plasmid DNAs were extracted from colistin-resistant Gramnegative bacterial strains (E. coli, K. pneumoniae, P. aeruginosa, and A. hydrophila) using Plasmid DNA Miniprep Kits (Thermo Fisher Scientific, Waltham, MA, USA) following the manufacturer's instructions. PCR amplifications of colistin resistance genes (mcr-1 to mcr-9) from phenotypic colistinresistant Gram-negative bacteria isolates were performed using oligonucleotide primer sequences and their annealing temperatures listed in Supplementary Table S1. Conventional PCRs were performed using a programmable 2720 thermal cycler (Applied Biosystem, Waltham, MA, USA) in a total reaction volume of 25 ml containing 12.5 mL of EmeraldAmp Max PCR Master Mix (Takara, Japan), 1 ml of each primer (20 pmol) (Sigma-Aldrich, Co., St. Louis, MO, USA), 5 ml of template DNA, and 5.5 ml of nuclease-free water. The amplified products were visualized after 30 min of electrophoresis on a 2% agarose gel containing ethidium bromide.

DNA Sequencing and Phylogenetic
Analysis of mcr Genes DNA sequencing of the PCR products was applied to validate their genetic identity and to verify any mutations in the mcr genes. The amplified DNA fragments were purified with the QIAquick PCR purification kit (Qiagen, Courtaboeuf, France). Sanger sequencing was performed in both directions using Bigdye Terminator V3.1 cycle sequencing kit (Perkin-Elmer, Inc. Waltham, MA, USA) in an applied Biosystems 3130 genetic analyzer (Hitachi, Tokyo, Japan). The nucleotide sequences were compared with those available in the National Center for Biotechnology Information (NCBI; www.ncbi.nlm. nih.gov). Nucleotide and deduced protein sequences were analyzed using the MEGA7 program . To investigate the genetic relatedness and the evolutionary distance among Gram-negative strains harboring the mcr genes, a phylogenetic tree was constructed using the neighbor-joining (maximum composite likelihood) method. DNA sequences generated in the study were deposited into the GenBank database with accession numbers MW811398-MW811434 and MZ648218-MZ648227.

Conjugation Assay
The transmissibility of colistin-resistant mcr genes between donors (mcr-positive isolates) and the recipient bacteria (E. coli J53; Na azide-resistant) was evaluated by the conjugation assays using the broth mating technique (Wang et al., 2003). An equal ratio (1:1) of donor and recipient cells (0.5 ml of each) in a logarithmic phase were added to 4 ml of sterile Luria-Bertani (LB) broth (Sigma-Aldrich, USA) then incubated overnight at 37°C without shaking. The mating mixture was serially diluted then the transconjugants were cultured on LB agar (Sigma-Aldrich, USA) supplemented with colistin (2 mg/L) or Na azide (100 mg/L). MICs for the donors, recipients, and transconjugants were determined by the broth microdilution method following the CLSI guidelines (CLSI, 2019). Potential transconjugants were examined for the existence of mcr genes using PCR assays. The conjugation efficiency was calculated as the number of transconjugants per donor as previously documented (Wang et al., 2003).

Disinfectant Susceptibility Testing
Gram-negative bacterial strains (three of each species exhibiting colistin resistance) were evaluated against four commonly used biocides in dairy farms for teat disinfection, namely chlorhexidine gluconate (0.5% in alcohol 70%; Pfizer, Manhattan, NY, USA), hydrogen peroxide (H 2 O 2 ; 6%, w/v), iodine (0.5%, w/v), and alcohol (ethanol; 70% v/v) (El Nasr Pharmaceutical Chemicals Co., Cairo, Egypt). The disinfectants were diluted (10-90%) using sterile distilled water, in which the original concentrations were considered stock solutions. The agar well diffusion method was applied to evaluate their effectiveness against the suspected strains of equivalent turbidity to that of 0.5 McFarland standards. Diameters of inhibition zones from the well edge to the inhibition front equal to or larger than 15 mm were regarded as susceptible (Wanja et al., 2020). The MICs of teat disinfectants against bacterial strains were determined as described before (Baron and Feingold, 1990).

Bioinformatics and Data Analysis
Statistical analyses were performed using the software SAS (SAS Institute Inc., 2012). Results of conjugation efficiency of Gramnegative bacteria were graphed by a boxplot through GraphPad Software (version 8.0.1, GraphPad Software Inc., La Jolla, CA, USA). Resistance phenotypes and frequencies of virulence genes among isolates were compared using the Chi-squared (c 2 ) test. P < 0.05 was considered significant. A heatmap with hierarchical clustering was generated to visualize the overall distribution of mcr variants and antimicrobial resistance phenotypes in colistinresistant isolates using the "pheatmap" package in R software (version 3.4.2) (Kolde, 2019). To determine the shared mcr variants among colistin-resistant isolates from different sources, a Venn diagram was generated using the Venny 1.0 tool, https://bioinfogp.cnb.csic.es/tools/venny/index.html. Spearman correlations were used to provide an estimate for the association among various variables (resistance phenotypes, colistin resistance, and virulence genes). The correlation coefficients and their p-values were visualized using a correlation plot. The correlation analyses and visualization were done using R packages rcorr (https://hbiostat.org/R/ Hmisc/) and corrplot (https://github.com/taiyun/corrplot). Genes or phenotypes that were present or absent in all analyzed subjects were not considered in these analyses.

Trends in Bacterial Resistance and MAR Indices
The overall resistance rates of Gram-negative isolates from clinical mastitis, subclinical mastitis, and raw milk to each tested antimicrobial were analyzed using c 2 test. As revealed in Table 1, a significant (p < 0.05) difference in resistance rates of isolates from clinical, subclinical mastitis, and raw milk was observed against ampicillin, cefazolin, ceftriaxone, nalidixic acid, tobramycin, fosfomycin, colistin, and tetracycline. All isolates from different sources were highly resistant to the first and second generations of cephalosporins especially cephalothin (100%). Furthermore, the MAR index of each antimicrobial was calculated to understand the level of antibiotic use. High MAR index (0.04) of ampicillin, amoxicillin-clavulanic acid, and the first and second generations of cephalosporins was observed. This implies the excessive use of these antimicrobials in the veterinary sector and human medicine. Meanwhile, imipenem achieved a lower MAR index.
A Venn diagram (Supplementary Figure S1) was generated to determine the shared mcr variants between colistin-resistant isolates from clinical mastitis, subclinical mastitis, and raw milk. E. coli mcr-4 and both mcr-7 and mcr-2 of P. aeruginosa were unique for clinical mastitis; whereas K. pneumonia mcr-3 and A. hydrophila mcr-1 to mcr-3 variants were exclusive for raw milk. Only E. coli mcr-3 and K. pneumonia mcr-2 were common among all sources.

Phylogeny of Colistin Resistance mcr Variants
As presented in Figure 2, the phylogenetic tree was constructed based on nucleotide sequences for all identified colistin resistance mcr variants (mcr-1 to mcr-4 and mcr-7) in 47 Gram-negative bacteria isolates of various species. The phylogenetic analysis revealed five potential clusters indicating sequence heterogeneity among the mcr variants. Despite the overall diversity of bacterial isolates, which were originating from different on-farm locations/ sources, the phylogenetic tree demonstrated the presence of the same clonal lineages of each mcr variant separately. Nonetheless, two mcr variants (mcr-4, accession number MW811433, and mcr-7, accession number MW811434) were clonally related; whereas the latter shared the same lineage with another mcr-2 variant (accession number MW811412).

Correlations Between the Antimicrobial Resistance and Virulence Traits
The correlation plots among pairs of antimicrobial resistance phenotypes, colistin resistance genes (mcr), and virulence genes in E. coli, K. pneumoniae, P. aeruginosa, and A. hydrophila are displayed in Supplementary Figure S2 and Supplementary Tables S2A-H. Significant positive correlations were pronounced between fluoroquinolone resistance and the existence of mcr-1 gene in E. coli (r = 0.58-0.67; p < 0.001).

Efficacy of Disinfectants on Colistin-Resistant Gram-Negative Bacteria
According to the agar well-diffusion results, certain disinfectants had a wide efficacy against the test bacteria. However, some disinfectants were ineffective even at high concentrations (Supplementary Table S3). E. coli and K. pneumonia isolates were resistant to chlorohexdine gluconate at concentrations lower than 50% (MIC = 40%), whereas P. aeruginosa and A. hydrophila were sensitive to such disinfectant even at low concentrations (MIC = 10%). E. coli and K. pneumoniae were sensitive to iodine at high concentrations (MIC = 80%), while P. aeruginosa and A. hydrophila exhibited resistance at all tested concentrations. On the other hand, hydrogen peroxide had the best inhibitory effect against all tested bacterial isolates even at a low concentration (10%), except for E. coli, which was sensitive to H 2 O 2 at a concentration of 50%. All tested bacteria were resistant to ethanol at all concentrations.  TET, tetracycline; MAR, multiple antibiotic resistance; NA, not applicable. a MAR index for each antimicrobial = total number of resistance scored/total number of antimicrobials tested × total number of isolates (Tambekar et al., 2006). *Trends in bacterial resistance to the tested antimicrobials using c 2 test; p < 0.05 was considered significant.

DISCUSSION
The recent global dissemination of mcr genes has raised the alarm of colistin resistance as a serious and urgent threat to public health (WHO, 2019). Consequently, large numbers of studies were carried out to explore the distribution of mcr genes in different reservoirs. Herein, we studied colistin-resistant Gram-negative bacteria isolated from bovine milk concerning the mcr resistance determinants and virulence characteristics.
Unfortunately, 52% of the tested isolates were colistin resistant. The high proportion of colistin-resistant bacteria and the diversity of mcr variants detected here indicates a lot of independent acquisitions of mcr genes by colistin-sensitive bacteria and confirm the hypothesis that animals and food chain may play a key role in mcr transmission (Luo et al., 2020). This is considering an alarming ratio since previous studies have detected colistin resistance in 37%, 40%, and 41% of bacterial isolates from pigs, poultry, and human clinical samples, respectively (Savin et al., 2020;Lay et al., 2021;Qadi et al., 2021). Of note, we identified these isolates in milk samples, which included various bacterial species. These findings could be attributed to the frequent usage of colistin in dairy farms that in turn could be responsible for the development of colistin resistance. This is consistent with a previous report of a Canadian dairy farm, which uses colistin to increase the milk yield although its usage is not officially licensed in veterinary medicine (Saini et al., 2012;Webb et al., 2017). Colistin-resistant E. coli isolates exhibited resistance to cephalothin, cefazolin, ceftazidime, ampicillin, amoxicillinclavulanic acid, cefoxitin, cefotaxime, piperacillin-tazobactam, and fosfomycin. Concurrently, nearly similar patterns were obtained recently for E. coli against several antibiotics in India, Croatia, and Egypt (Zdolec et al., 2016;Singh et al., 2018;Ramadan et al., 2021;Soliman et al., 2021). For K. pneumoniae isolates, the antibiotic resistance pattern was closely like the pattern of Brazilian MDR K. pneumoniae isolates (Ferreira et al., 2019). They showed resistance to cefazolin, cephalothin, cefoxitin, ceftazidime, cefotaxime, ceftriaxone, sulfamethoxazole-trimethoprim, and tigecycline. On the contrary, our colistin-resistant P. aeruginosa isolates showed resistance to most groups of antibiotics; however, previous colistin-resistant P. aeruginosa isolates from Iran showed only resistance to ticarcillin, ciprofloxacin, aztreonam, ceftazidime, gentamicin, imipenem, and piperacillin/tazobactam (Goli et al., 2016). Our A. hydrophila isolates have an antimicrobial resistance pattern much akin to the antibiotic resistance profiles of pathogenic aeromonads isolated from ornamental fish with resistance to gentamicin, cephalothin, cefoxitin, ceftazidime, cefotaxime, ceftriaxone, nalidixic acid, sulfamethoxazole-  AK, CN, TOB, CEF, FOX, CAZ, CFZ, SXT, NA, CIP, LVX, AM  FIGURE 2 | A phylogenetic tree shows the genetic relationship between mcr gene sequences of colistin-resistant Gram-negative bacteria recovered from clinical, subclinical mastitis, and raw milk. The tree was constructed using a neighbor-joining method with a bootstrap value for 1,000 replicates.    Table 2.
trimethoprim, aztreonam, and chloramphenicol (Saengsitthisak et al., 2020). Colistin-resistant isolates found in this study exhibited extreme emergence and spread of antimicrobial resistance and high virulence as well.
It is almost established that E. coli is deemed to be the topmost bacteria for the mcr genes by carrying mcr-1 up to mcr-5 with the chance for coexistence of more than one mcr gene (Yin et al., 2017;Hammerl et al., 2018). Similarly, colistin-resistant E. coli isolates in our study were found to carry mcr-1, mcr-2, mcr-3, or mcr-4 genes; 20% out of which were from raw milk and carried mcr-1, mcr-2, or mcr-3 genes. It is the first detection of mcrcarrying isolates in raw milk samples. A previous study has identified mcr-1 gene in one E. coli isolate from a cow suffering from subclinical mastitis in Egypt (Khalifa et al., 2016). Although the mcr-1 gene was detected frequently from different samples (Khalifa et al., 2016;Lima Barbieri et al., 2017;Ramadan et al., 2021), it is foremost found here in the raw milk samples. Marvelously, one E. coli isolate carrying mcr-1 gene obtained from hard cheese made from raw milk was reported in Egypt very lately (Ombarak et al., 2021), which in turn lightens the importance of raw milk quality control checks in Egypt, focusing on farm workers as well as the animal itself.
The prevalence of mcr-3 gene in our isolated E. coli was the dominant one, followed by the mcr-2 gene. Only one E. coli isolate was mcr-4 positive. However, in this study, the mcr-1, mcr-2, and mcr-3 genes were found in E. coli isolated from raw milk as well as clinical and subclinical mastitis samples; however, the only mcr-4-positive E. coli was isolated from a clinical mastitis sample. The variants mcr-2, mcr-3, and mcr-4 were identified formerly in E. coli isolates from Belgium, China, and Italy, respectively, and were isolated from calves and pigs (Xavier et al., 2016;Carattoli et al., 2017;Yin et al., 2017). All these data confirmed that E. coli is the chief bacteria carrying mcr genes from different samples and in different geographical areas.
Notably, all our colistin-resistant A. hydrophila were isolated from raw milk samples, and the predominant gene in them was mcr-1 followed by mcr-2 and mcr-3. In contrast, colistin-resistant A. hydrophila were found to harbor mcr-7 gene in cooked and raw meats and in environmental samples (Liu Y. et al., 2020). No previous reports until now showed the isolation of colistinresistant A. hydrophila from raw milk samples. Our findings warrant that the source of raw milk contamination is most probably the farm workers rather than the animal source, particularly, that mcr-3 gene was previously identified in isolates from human rectal swabs (Shen et al., 2018).
Of the identified colistin-resistant isolates, 19.12% were K. pneumoniae and were carrying mcr-1, mcr-2, or mcr-3 genes. The three variants were identified in the isolates from raw milk samples, while mcr-1 was found in subclinical mastitis isolates and mcr-2 in both isolates from clinical and subclinical mastitis. It is almost agreed that the mcr-1 gene is one of the few and clear examples of animal origin of a resistance trait that could hit the entire human health system (Poirel and Nordmann, 2016). Appealingly, three mcr-1-positive K. pneumoniae isolates were isolated from raw milk samples ensuring the role of humans in contamination and the transmission from animals to humans.
On the other hand, no colistin-resistant P. aeruginosa were found in the raw milk samples. They were isolated only from clinical and subclinical mastitis and were carrying mcr-1, mcr-2, mcr-3, or mcr-7 genes. Although, we detected 37.5% and 12.5% of the colistin-resistant P. aeruginosa isolates carrying mcr-1 and mcr-2 genes; respectively, it was found lately in Egypt that colistin-resistant P. aeruginosa were negative for mcr-2 gene and 50% were positive for mcr-1 (Abd El-Baky et al., 2020). We also found that the mcr-3 variant was prevalent with 37.5% among our colistin-resistant P. aeruginosa isolates. Even, to our knowledge, there is no report detected on mcr-3 in P. aeruginosa isolates before. Moreover, this is the first-time reporting mcr-7 from Egyptian isolates, and it was identified in P. aeuriginosa from clinical mastitic milk. The first plasmid-mediated colistin resistance gene mcr-7 was detected in K. pneumoniae from chickens in China . Recently, the mcr-7 gene was detected in a water sample from an environment of an alligator (dos Santos et al., 2020). The detection of mcr-7 in our isolate from clinical mastitis confirms the dissemination of mcr gene and its variants at all levels; humans, animals, and environmental. This was confirmed here while application of the conjugation assay as the plasmid harboring mcr genes (mcr-1 to mcr-4, and mcr-7) successfully transferred colistin resistance to the recipient strain through a broth-mating assay, representing their capability to mobilize the mcr-genes between isolates (Migura- Garcia et al., 2020).
The virulome of our colistin-resistant isolates was screened in each species. Many of the identified virulence genes in our study were required for bacterial survival and had a role in causing mastitis (Tartor and El-Naenaeey, 2016;Gao et al., 2019;Guerra et al., 2020). This in turn raises the alarm that most of the colistinresistant Gram-negative isolates are highly pathogenic. Here, the establishment of an interplay between antimicrobial resistance and virulence traits corresponded to each analyzed Gram-negative bacteria, revealing positive correlations between the existence of certain virulence genes and resistance to antimicrobial agents/ genes. A previous study suggested that the mobile genetic elements of proficient pathogens favor the coselection of both resistance and virulence-associated traits (Hennequin and Robin, 2016).
The virulome, as well as the resistome patterns of our colistinresistant isolates, confirmed the elevated pathogenicity and the widespread antibiotic resistance that correlated to these isolates. The phylogenetic analysis of the nucleotide sequences of our colistin resistance mcr variants showed three distinguishable clusters. Interestingly, each cluster has the same mcr variant, and this is in accordance with what was reported before about mcr-1, mcr-2, and mcr-3 . The presence of each variant of the mcr gene in the same clonal lineage despite the different sources of isolation confirms the high spread of colistin resistance from one source to another. Each of the mcr-4 and mcr-7 variants was shown in a unique cluster. This was unsurprising since our study showed the first report of mcr-7 in P. aeruginosa. Similarly, mcr-4 was previously found in a unique phylogenetic group . Our findings proved that each mcr gene has a distinct evolution and spread path producing a characteristic phenotypic colistin resistance.
All the above information has raised the alert for controlling the spread of colistin-resistant isolates. Thus, we decided to evaluate four commonly used biocides in the dairy farms against colistin-resistant isolates. In accordance with Lanjri et al. (2017), all isolates were sensitive to chlorhexidine with different concentrations. This is attributed to its cationic nature that helps in electrostatic interaction with the anionic group lipopolysaccharide in the cytoplasmic membrane of Gramnegative bacteria, resulting in modification of the membrane permeability and causing membrane disruption and cell death due to coagulation of the cytoplasm (Estrela et al., 2003). Another tested biocide was iodine; it showed an inhibitory effect on E. coli and K. pneumoniae at high concentrations and no inhibition effect on P. aeuriginosa and A. hydrophila at all tested concentrations. Iodine is proven to be considered a broadspectrum, rapidly acting bactericidal, and effective against all mastitis-causing bacteria (Gleeson et al., 2009). It acts through its penetration into the proteins, nucleotides, and fatty acids of microorganisms causing cell death (Gupta et al., 2018), which makes the development of resistance unlikely. However, in the field, the presence of some residues in the environment like the organic matter of the teat skin or in the challenge suspension may have been reacting with the disinfectant and thus less germicidal activity may have been produced (Enger et al., 2015). These conditions permit some bacteria to be exposed too much to lower concentrations of the active ingredients consequently and thus increasing their resistance to the used disinfectants (Wales and Davies, 2015). These results were in agreement with what was reported earlier that a dip containing the highest free iodine provided a greater significantly superior reduction of pathogens compared with a lower concentration of iodine dip (Foret et al., 2005). Furthermore, all isolates showed a sensitivity to H 2 O 2 , and this had been confirmed previously. H 2 O 2 produces oxidants as singlet oxygen, superoxide radicals, and hydroxyl radicals with high penetration power and not specific in the points of contact by attacking bacteria at any essential cell components, leading to a decrease in the probability of bacteria to create resistance with these compounds (Nilima, 2011;Benavides, 2014). This was in contrast with Enger et al. (2015) who reported that H 2 O 2 was the least efficacious tested teat dip, and these differences may be due to the evolution of resistance by the gradual use of the disinfectant (Wales and Davies, 2015). However, different concentrations of ethanol showed no antimicrobial activity against all tested bacteria, and this is in agreement with what was proved before about ethanol effectiveness (Al-Dabbbagh et al., 2015). This ethanol tolerance mainly appears to result in large part from adaptive and evolutionary changes in cell membrane composition (Ingram, 1990).

CONCLUSION
Ultimately, the spread of mobile mcr genes between MDR and XDR virulent Gram-negative isolates from dairy cattle confirms the dissemination of the mcr gene and its variants at all levels: human, animal, and environmental. Consequently, a decision to stop using colistin in food animals is a principal demand to defy XDR and MDR emergence. Additionally, establishing a perfect quality control check system in Egypt, focusing on both farm workers and the animal itself is urgent.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

AUTHOR CONTRIBUTIONS
YT and NA contributed equally in the conception and design of the study and participated with RG, HE, EK, and AA in the application of classical microbiological techniques. YT and NA carried out the PCR assays, sequencing approaches, and participated with HR in sequence analysis. HR performed the bioinformatics and statistical analyses of the data. RG, HE, SE, EK, AA, SA, and HR conceived the study and participated in the design. YT, RG, NA, HE, SE, and AA wrote the initial draft of the manuscript. All authors revised the manuscript critically for important intellectual content and gave the final approval of the version to be published.

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2021. 761417/full#supplementary-material Supplementary Figure S1 | | A Venn diagram showing the intersection between the mcr variants of colistin-resistant isolates from clinical, subclinical mastitis, and raw milk.
Supplementary Figure S2 | | Correlation plot between antimicrobial resistance phenotypes, colistin resistance genes (mcr) and virulence genes in E. coli, K. pneumoniae, P. aeruginosa, and A. hydrophila. In each plot, each cell represents the R-value (correlation coefficient) for each pair of features. The color and circle size indicate the value of correlation (red = negative and blue = positive correlation). The more intense the color, the higher the correlation value. Stars referring to significant correlation determined at a cutoff = 0.05. Insignificant correlations are left without stars.