Abstract
Bovine mastitis remains the main cause of economic losses for dairy farmers. Mammary pathogenic Escherichia coli (MPEC) is related to an acute mastitis and its treatment is still based on the use of antibiotics. In the era of antimicrobial resistance (AMR), bacterial viruses (bacteriophages) present as an efficient treatment or prophylactic option. However, this makes it essential that its genetic structure, stability and interaction with the host immune system be thoroughly characterized. The present study analyzed a novel, broad host-range anti-mastitis agent, the T4virus vB_EcoM-UFV13 in genomic terms, and its activity against a MPEC strain in an experimental E. coli-induced mastitis mouse model. 4,975 Single Nucleotide Polymorphisms (SNPs) were assigned between vB_EcoM-UFV13 and E. coli phage T4 genomes with high impact on coding sequences (CDS) (37.60%) for virion proteins. Phylogenetic trees and genome analysis supported a recent infection mix between vB_EcoM-UFV13 and Shigella phage Shfl2. After a viral stability evaluation (e.g pH and temperature), intramammary administration (MOI 10) resulted in a 10-fold reduction in bacterial load. Furthermore, pro-inflammatory cytokines, such as IL-6 and TNF-α, were observed after viral treatment. This work brings the whole characterization and immune response to vB_EcoM-UFV13, a biocontrol candidate for bovine mastitis.
Similar content being viewed by others
Introduction
Bovine mastitis remains the main cause of economic losses for dairy farmers, estimated at $ (US)533 billion worldwide, as well as public health concerns, since low quality milk can be considered a vehicle for pathogen transmission1,2,3.
Mastitis treatment is typically based on the use of short and long-acting antibiotics, respectively, during the lactation and dry period4. In terms of the lactation period and specifically regarding clinical mastitis, Escherichia coli, Streptococcus uberis, Streptococcus dysgalactiae and Staphylococcus aureus are the main etiological agents involved that have been routinely isolated5. Among these pathogens, mammary pathogenic Escherichia coli (MPEC) is responsible for an acute mastitis characterized by inflammation, increased somatic cell count (SCC) and impaired milk quality even after the infection has been cured6.
Using E. coli strains obtained from different types of mastitis (e.g per-acute and persistent) and the non-pathogenic strain K71, Blum et al. (2017) performed a mammary immune response comparison in experimentally infected cows and noticed differences regarding TNF-α, IL-6 and IL-17 secretion levels for each MPEC.
Within the current scenario of widespread antibiotic-resistant bacteria, bacterial viruses present as an efficient therapeutic or prophylactic tool in order to control different pathogens in dairy cows at different lactation stages. This is supported by current in vitro and in vivo assays7,8,9,10,11,12,13.
Considered a model system in molecular biology, coliphage T4 has been studied since the 1940s and possesses about 300 genes organized in a 168.9 kb linear dsDNA with an average GC content of 34.5%14. Lytic viruses related to T4 have awoken interest for their application in phage therapy due the absence of lysogenic modules, a broad-host-range (from Proteobacteria to Cyanobacteria phyla) and, recently, the identification of virion-associated peptidoglycan hydrolases (VAPGHs), which are considered potential enzybiotics15,16,17.
Indeed, T4 was used by Bruttin18 in the first safety test of phage therapy in humans and in immunological assays in order to elucidate the interactions between viruses, its host and the immune system. Bocian et al.19 investigated how purified T4 phage and T4-generated E. coli lysate impact immune cells differentiation, highlighting that lysis of Gram-negative bacteria by phages might not trigger excessive monocyte induction.
Currently, several studies have been exploited the use of mouse models in the attempt to evaluate novels anti-mastitis drugs mainly against E. coli and S. aureus20,21,22,23,24,25,26. The use of animal models is time and cost effective approach, along with a previous step for pre-clinical and clinical assays27. However, the absence of intrinsic cow factors can be considered a bottleneck when data is analyzed28.
The aims of the present study were to characterize the Escherichia phage UFV13, a T4virus, in genomic, protein and physiological terms and evaluate the immune response in an experimental E. coli-induced mastitis mouse model with the aim to use it in clinical trials to control mastitis in dairy cows.
Materials and Methods
vB_EcoM-UFV13 isolation and purification
vB_EcoM-UFV13 (UFV13) was obtained from the sewerage system of Viçosa, Minas Gerais state, Brazil and was propagated on Escherichia coli 30 following the well-established Sambrook & Russell29 protocol. This virus belongs to the bacterial virus collection of the Laboratório de Imunovirologia Molecular (LIMV) at Universidade Federal de Viçosa (UFV), Viçosa city, Brazil.
After virus propagation, viral particles were purified by a three-step protocol using ion exchange and desalting columns in a chromatography system (ÄKTAprime plus, GE Healthcare Life Sciences, Uppsala, Sweden). Briefly, an initial step using HiTrap Desalting prepacked column (GE Healthcare Life Sciences, Uppsala, Sweden) was conducted to remove any salts used at virus propagation stages and the first two peaks were collected and purified using an ion exchange chromatography column, with specific fractions (7 to 11) driven to the final step, a new desalting process. For the anion exchange column, start (20 mM Tris-HC, pH 8.0) and elution buffers (20 mM Tris-HC, 1 M NaCl, pH 8.0) were used, whereas in desalting steps a phosphate buffer (20 mM sodium phosphate, 0.15 M NaCl, pH 7.0) was prepared. Flow rate of 5 mL.min−1 was adopted for both columns. Finally, viral titer was measured at 37 °C by double-agar overlay method30 using E. coli 30 as the plating host. Phage stocks were stored at 4 °C for further analysis.
Bioinformatic analysis
Phage genome extraction, sequencing and annotation methodologies are described according to Duarte et al.31. Hypothetical proteins were manually checked for homologs using UNIPROT database. Protein isoelectric point and molecular weight were obtained using ExPASy32. Putative tRNAs and Rho-independent transcription terminators were, respectively, predicted using tRNAscan-SE33 and ARNnold web tool34. The Database of Gene Regulation in Bacteriophages (phiSITE)35 was used in order to check the three major classes of promoters (early, middle and late) as well as to confirm putative Rho-independent terminators. The CGView Server was used to generate a UFV13 genome graphical map36.
Nucleotide differences between UFV13 and Enterobacteria phage T4 (accession number AF158101.6) genomes were performed by checking reading alignments. High-quality Illumina reads were filtered and adaptor sequences were removed using Trimmomatic software (ver 0.33)37 (parameters: LEADING:10 TRAILING:10 SLIDINGWINDOW:4:15 MINLEN:65) and aligned to Enterobacteria phage T4 using Bowtie2 software (v2.2.4)38. SAMtools39 was used to convert the output SAM format to BAM format and, subsequently, to sort the BAM file. The sorted BAM file was processed with mpileup tool (SAMtools package) in order to extract the variants. The Binary Call Format (BCF) created was converted to VCF format using BCFtools39. VCF file and Enterobacteria phage T4 genes were used as input for SnpEff program40. Only variants with predicted “high” or “moderate” impact on the protein-coding gene were analyzed.
Whole-genomes that represent each genus from the subfamily Tevenvirinae (T4virus, Cc31virus, S16virus, Js98virus and Sp18virus) deposited on the International Committee on Taxonomy of Viruses (ICTV) were downloaded from NCBI (accession numbers are provided in Supplementary Table 4) and aligned with UFV13 genome using Progressive MAUVE41. Whole-genome single nucleotide polymorphisms (SNPs) were extracted as described in Treu et al.42 and a SNP-based phylogenetic tree was constructed using PHYLIP package43 and visualized by dendroscope44. Furthermore, a whole-genome phylogenetic tree was formulated aligning previously cited genomes from the Tevenvirinae subfamily. Online tools ClustalW245 and MAFTTT version 746 were used for multiple alignment. A neighbor-joining tree was drafted with MEGA747 and visualized with FigTree (http://tree.bio.ed.ac.uk/software/figtree/).
In silico bacterial hosts of UFV13 were predicted using HostPhinder48.
vB_EcoM-UFV13 structural protein analysis
With the aim of obtaining the UFV13 structural protein profile, proteomic analysis was conducted.
Viral propagation/purification (see item 2.2) were conducted and phage particles concentrated by adding NaCl 1 M and polyethylene glycol 8000 (PEG8000) 10% (w/v). The mixture was kept overnight at 4 °C. After a centrifugation (10,000 × g, 15 min), viral pellet was resuspended in 1 mL of SM buffer and one volume of chloroform followed by centrifugation at 4,000 g for 10 min. An equal volume of trichloroacetic acid (TCA) 10% (v/v) was added to the supernatant and incubated on ice for 30 min. The precipitated viral proteins were collected by centrifugation (11,000 × g, 20 min), washed three times with acetone (11,000 × g, 10 min) and resuspended in water. The BCA Protein Assay Kit (Boster Biological Technology, Wuhan, China) was used to estimate protein concentrations.
For protein profiling, the viral proteins were separated via sodium dodecyl sulfate-polyacrilamide electrophoresis (SDS-PAGE) on a15% gel49. Protein bands were highlighted by Coomassie Brilliant Blue R-250 dye and removed as described by Shevchenko et al.50. Mass spectrum (MS) was acquired by matrix-assisted laser desorption-ionization time of flight (MALDI/TOF-TOF) (Ultraflex III - BRUKER DALTONICS). Spectra intervals between 500 and 3,500 kDa were selected and forwarded for MS/MS analysis. Results were obtained by Mascot™ (Matriz Science) software using the NCBl nr protein database. Only proteins and peptides indicated as significant by Mascot were considered for further analysis.
Physiological features
Physiological features were assessed in triplicate following the protocols described by Jurczak-Kurek et al.51. In each of the following cases, following treatment the phage preparation was diluted and tittered as described above (2.2).
pH stability
Viral capability at different acidic and alkaline pH values was evaluated. One mL of purified viruses were transferred to 9 mL LB medium with pH 2, pH 4, pH 7 (control), pH 10 and pH 12 at a 1:9 ratio and incubated for 1 h at 37 °C, being proceeded by a 10-fold serial dilution and plating as described in item 2.2. After overnight incubation at 37 °C, virus stability was determined by the percentage of viruses able to produce lysis plate.
Thermal stability
Thermal stability studies on LB-diluted phage suspensions were made at −20 °C, 40 °C, 62 °C and 95 °C, respectively, for 12 h, 40 min, 40 min and 5 min. After these periods, a serial 10-fold dilution was conducted and a specific aliquot was plated and incubated overnight at 37 °C. Viruses that did not undergo any thermal treatment were used as a control.
Viral ability to propagate in different temperatures
Viral replication at different temperatures (4 °C, 22 °C, 30 °C and 37 °C) was assessed by spot-assay after 24 hours of incubation using a 10-fold serially diluted viral stock in LB medium.
Osmotic shock effect
In order to verify the effect of the osmotic shock on virus particles, a stock aliquot was transferred to TM buffer (10 mM Tris–HCl, 10 mM MgSO4; pH 7.2) with sodium chloride 4.5 M, incubated at room temperature for 15 min and quickly diluted in TM buffer without sodium chloride. Bacterial viruses incubated in TM buffer without sodium chloride were used as a control.
Antiviral resistance
Antiviral activity of sodium dodecyl sulfate (SDS), sodium lauroyl sarcosinate (Sarkosyl) and cetyltrimethylammonium bromide (CTAB) on virus particles was determined incubating a standardized viral suspension with 0.09% SDS (20 min at 45 °C), 0.1% CTB (1 min at 22 °C) and 0.1% Sarkosyl (10 min at 22 °C). Controls were done considering the same conditions for each antiviral compound but they were substituted for TM buffer.
Organic solvent effect
To study the effect of four different organic solvents, a viral suspension was added to 63% ethanol, 90% acetone, 90% chloroform and 50% dimethyl sulfoxide (DMSO). Mixtures were incubated for 1 h at 22 °C (ethanol and acetone), 1.5 h at 4 °C (chloroform) and 10 min at 4 °C (DMSO). In the next step, 10-fold dilutions in TM buffer (10 mM Tris–HCl, 10 mM MgSO4; pH 7.2) were prepared and used for plating. Phages incubated in TM buffer under conditions described above, were used as a control.
Animal model and immune response
Escherichia coli 30 strain
The mammary-pathogenic E. coli 30 was isolated from a dairy cow with acute mastitis and kindly provided by Brazilian Agricultural Research Corporation (Empresa Brasileira de Pesquisa Agropecuária – EMBRAPA) Dairy Cattle (Juiz de Fora, Minas Gerais, Brazil). E. coli 30 was evaluated for its ability to form biofilm as well as for motility capacity. Motility assays were conducted according to Deziel et al.52. Briefly, an overnight E. coli 30 aliquot was washed, suspended in distilled and sterilized water, and inoculated in King B medium (peptone 20 g/L, MgSO4.7H2O 1.5 g/L, K2HPO4 1.5 g/L), supplemented with 1.5, 0.5 and 0.3% of agar for twitching, swarming and swimming tests, respectively. Biofilm assay followed the common crystal violet (CV) staining method53. After incubation (37 °C for 48 h), wells were washed 3 times with PBS buffer to remove not adherent cells. CV was added at a minimum volume capable overcome bacterial suspension volume, making possible quantify all bacterial biomass, and incubated for 30 minutes at room temperature. To biomass quantify, CV was removed, wells washed with PBS, added ethanol to solubilize biomass internal CV crystals, and optical density measured at 560 nm. All assays were performed in triplicate.
Antimicrobial susceptibility test of E. coli 30 was assessed by disc diffusion assay using polisensidiscs for 25 different antibiotics (DME, Araçatuba, São Paulo, Brazil) following the manufacture’s recommendation. The results were interpreted according to the standards of the Clinical and Laboratory Standards Institute (CLSI)54.
E. coli-induced mastitis in mouse
With the aim to evaluate UFV13 effectiveness and the immune response against E. coli 30 in vivo, an experimental E. coli-induced mastitis mouse model was used. The trial followed the methodology described by Chandler55, with some modifications, and was approved by the Ethics Committee (Comissão de ética no uso de animais/UFV) according to the protocol 64/2016. Lactating Balb/c female mice (5 to 15 days) was intraperitoneally anesthetized with 10% Ketamine and 2% Xylazine, with subsequent surgical operation of the mammary gland by cutting the teat canal of the last two abdominal ceilings (R5 and L5 of each animal were used) (Supplementary Figure 4). E. coli 30 (100 UFC/ml), PBS and phage served as the control groups, while phage plus bacteria was considered the treatment group (MOI 10). Viral addition was done four hours after bacterial inoculation. Three animals were used to perform each experimental group.
In order to assess the lytic activity of the phage vB_EcoM-UFV13 in mammary glands, animals were euthanized 48 hours after the treatment by anesthetic overdose. Infected mammary glands were removed and transferred to 1.5 ml of (PBS), macerated and serially diluted (1:10) in PBS buffer. A microdrop assay56 was performed to estimate E. coli 30 colony-forming unit (CFU).
Cytokines IL-6, TNF-α, IL-2, IFN-γ, IL-4, IL-10 and IL-17A obtained from macerated mammary gland (L5 and R5 from each animal was pooled) were simultaneously quantified by the Cytometric Bead Array kit (CBA, BD Bioscience) in a BD FACSVerse Flow cytometry following manufacturer’s recommendations.
Histological analysis
For histology, mice were euthanized 48 hours after phage treatment. The mammary glands were removed and fixed in Karnovsky fixative (paraformaldehyde 4% and glutaraldehyde 4%, pH 7.3). Further, the tissues were embedded in paraffin and a 5 µm section was obtained by microtomy, stained with hematoxylin and eosin (H & E) and observed under light microscopy. The images were acquired under an Olympus DP73 microscope.
Statistical analysis
Statistical analysis was performed with GraphPad Instat 3 software (GraphPad, La Jolla, CA, USA) using the one-way analysis of variance (ANOVA) at 95% accuracy level to evaluate the differences between mammary gland cytokine production under different conditions. This assay was set up in triplicate and, for each animal, the two abdominal ceilings were pooled into one sample. Tukey’s test was used as post hoc test.
Results and Discussion
vB_EcoM-UFV13 genome analysis
Bacteriophage UFV13 was isolated from samples obtained in the sewage system of Viçosa, Minas Gerais, Brazil, a well-known source of novel viruses57.
From the genome sequence (165,772 bp, GC content 34.8%), 269 ORFs were predicted and annotated. The size is similar to that of other members of the T4virus genus58. In addition, a total of ten tRNA encoding genes (Gln, Leu, Gly, Pro, Ser, Thr, Met, Tyr, Asn and Arg) were identified and are organized in a gene cluster without introns or pseudogenes (Fig. 1). 13 ORFs showed an identity below 70% with the reference phage T4, while 24 ORFs were annotated as hypothetical protein coding sequences. ORFs 50 and 232 are, respectively, coding sequences for T2 and T6 bacteriophage proteins. 13 ORFs encode for Shigella phage proteins (pSs-1, Shfl2, Shf125875, SH7, SHBML-50-1). Genes for IpII and IpX were not identified in the UFV13 phage genome. Gene and protein information such as genomic coordinates, protein weight, pI and putative function for each UFV13 ORF are reported in Supplementary Table 1.
Functional categorization of UFV13 genes (Supplementary Table 2) revealed that nonessential and auxiliary genes mainly related to homing endonucleases of introns such as I-TevI, I-TevII and I-TevIII are absent, as with some of their related introns like mob and seg genes (Table 1). Homing endonucleases of introns are considered systems associated with gene conversion events, transference of mobile elements and gene exclusion in mixed infections14. In fact, genome analysis suggests a recent mixed infection among UFV13 and Shigella phage Shfl2, once segF was substituted by soc.1 and soc.2 from Shigella phage Shfl2. According to Belle et al.58, segF was absent in T2 phages, but the region is occupied by soc.1 and soc.2.
Similarity analysis using the deposited genomes from Yersinia phage PST and Shigella phage Shfl2 showed that the bacteriophage UFV13 has 97% identity with both phages. However, it has a larger percentage of aligned genome (96%) for the Shfl2, and 90% of genome aligned to phage PST (Supplementary Figure 1). Comeau et al.59, characterized Yersinia phage PST genome, containing dsDNA with 167,785 bp, 35.3% GC content and 9 tRNAs, whose values are near to those found for the phage UFV13, with 34.8% GC content. Jun et al., (2016) identified 10 tRNAs for both pSs-1 and Shfl2 viruses infecting Shigella spp. as host. Despite little existing knowledge about tRNAs functioning in phage life cycles, high tRNAs numbers found in some bacteriophages could be proportional to phage genomes sizes and are related to a short latent period and high burst size value60.
In the total, 35, 21 and 31, early, middle and late promoter regions respectively, were predicted, which corresponds to 73% of the promoters identified for the T4 genome (Supplementary Table 3). As expected, no promoter regions for mob and seg genes were found, along with I-TevII and internal protein II (ipII).
Overall, termination of transcription in T4 bacteriophages occurs by an intrinsic termination signal, a stem-loop arrangement accompanied by a U-rich region, or a Rho-dependent termination mechanism61. For UFV13, 82 Rho-independent transcription terminators were predicted and compared against phiSITE database. 30 sequences displayed high identity with the reference phage, while five were not predicted but were found using the sequences available on the phiSITE. Terminators for the genes 45, 5.4/6′, 34/35, 35, 37, nrdA, uvsY.−2 and 56/segF were not found or even identified using phiSITE (Table 2), while for the genes such as regA, wac, 24(b) and 30.9(b) were, respectively, detected on the ORFs 52, 163, 178 and 208. According to Miller et al. (2003), generally 34 terminators are found in the T4 genome. The absence of terminators for the genes nrdA and 56/segF can be explained by their absence, since these genes were not annotated. The significance of terminators within genes is unknown but have been described for T4 phages14.
A total of 5,071 filtered variants (4,975 SNPs, 86 insertions and 10 deletions) were identified (1 variant every 33 bases) between UFV13 and T4. 83,656 effects were assigned and categorized (395 high (0.472%), 1,670 low (1.996%), 2,722 moderate (3.25%) and 78,869 modifier (94.28%). The number of effects by functional class is: missense 2,879 (60.79%), nonsense 187 (3.95%) and silent 1,670 (35.26%). Considering only genes with high or moderate variants, the mainly affected protein category was that related to virion structural proteins, with gp7 being the most affected (Fig. 2). This protein is considered the second largest protein in the baseplate and is one of the seven components associated with wedge assembly and stability62. The large number of SNPs found in virion proteins reflects the diversification of phage UFV13 from the classical T4 and can be associated with the high capability to infect different bacterial genus such as Escherichia, Morganella and Shigella, an uncommon phage feature (data not shown).
In order to evaluate the phylogenetic relationship between UFV13 and genera belonging to the subfamily Tevenvirinae (T4virus, Cc31virus, S16virus, Js98virus and Sp18virus), a whole-genome tree was constructed. The phylogenetic tree divided phages into four clusters with vB_EcoM-UFV13 possessing the closest relationship with Shigella Shfl2 and Yersinia phage PST (Fig. 3), a result also supported by the SNP-based phylogenetic tree (Supplementary Figure 2).
With the aim to establish the host range of phage UFV13 host range, in silico Host Phinder test was conducted. Four bacterial genera were predicted to be infected by UFV13: Escherichia (E-value: 5.9e−1), Yersinia (E-value: 5.8e−1), Shigella (E-value: 6.3e−1) and Salmonella (E-value: 1.1e−2).
Viral protein analysis
SDS PAGE analysis revealed the presence of eight proteins (Supplementary Figure 3). Among these, four proteins were chosen and analyzed by MALDI/TOF-TOF, with the following being identified: UFV13_gp243 long tail fiber proximal subunit (139.95 kDa), E. coli chaperonin GroL (57.36 kDa), Escherichia phage vB_EcoM_112 major capsid protein (56.09 KDa) and E. coli outer membrane protein C (36.77 kDa). These proteins are deposited at UNIPROT under accession numbers A0A160CBJ3, A0A017I9Q1, A0A160CBB0 and A0A148HSV3, correspondingly. The presence of a common viral receptor (OmpC) and the chaperonin GroL highlights the need for viral purification enhancement. Viral isoelectric points are usually below six and some of them are able to bind to anion and cation exchange matrixes. However, host cell DNA is efficiently removed by cation exchanger columns, while host cell proteins are effectively eliminated by anion exchange63.
Physiological features
Viral ability to survive in a wide range of adverse conditions is a desired characteristic for therapeutic as well as a biological control agents64. Thus, phage stability was evaluated in different physical and chemical conditions.
UFV13 was relatively stable within a pH range of 7.0–12.0. An approximate 1 log-fold reduction on viral titer was observed after incubation at pH 4.0 (18% of survivability) (Fig. 4A). Interestingly, phage UFV13 was inactivated after 1 h of incubation at pH 2.0, but not at pH 12 (68% viral viability), which is indicative of considerable virion stability at basic pH values.
UFV13 was able to lyse E. coli at 30 °C and 22 °C with a plating efficiency of 69 and 56%, respectively (Fig. 4B). No lysis plates were observed after storage at 4 °C, which may be indicative of DNA injection without host lysis.
Viral thermal inactivation occurred in the extreme tested temperatures 95 and −20 °C (0 and 24% survival, respectively), while viral titer dropped to 50% when UFV13 was incubated for 40 min at 62 °C (Fig. 4C). No significant reduction of phage titer was observed after thermal treatment at 40 °C.
It was found that osmotic pressure change, detergent and organic solvents determined a significant titer drop of UFV13. A survivability of 12% was observed when UFV13 underwent a rapid dilution from high-concentration NaCl buffer to low-concentration ones (Fig. 4D). The anionic detergent Sarkosyl reduced the number of viable viral particles by 74%, while CTAB and SDS resulted in a complete inactivation of this virus (Fig. 5). Anti-phage activity was also evidenced after exposure to ethanol, chloroform, acetone and DMSO.
Physiological features evaluated in this study are in accordance with Jurczak-Kurek et al. (2016). In a broad physiological study using 83 bacteriophages isolated from urban sewage, the same source of UFV13, it was verified that the vast majority of phages were sensitive to a temperature of 62 °C (survival below 70%), were able to survive in basic pH (ranging from 10 to 12), were susceptible to detergents and organic solvents, except chloroform, and were also resistant to osmotic shock.
Viral stability assays have also been shown that UFV13 can survive in raw milk and it has potential capability to survive in mastitic milk. In general, mastitis can dampen the quality of raw milk composition, which includes increased levels of Na+/Cl− and pH65. Evaluating a bacteriophage cocktail composed of two T4 phages in raw milk against E. coli, Porter et al. (2016) observed a 3.3- to 5.6-log reduction of bacterial growth over a 12-h physiologic temperature. As discussed below (item 3.4), this study showed UFV13 activity against E. coli in lactating female mice.
E. coli-induced mastitis mouse model
E. coli 30, an isolate obtained from a dairy cow with acute mastitis, displays resistance to 14 of 25 evaluated antibiotics (Supplementary Table 5). This bacterium is resistant to at least three different antimicrobial drug classes, which allows its classification as a multi-drug resistant strain66 Moreover, of the four types of virulence factors analyzed (biofilm forming capability and motility types: swarming, swimming and twitching) the E. coli 30 was positive for biofilm formation, swarming, and swimming, being negative only for the twitching (Supplementary Table 6). Studies correlated this virulence factors to an increase in pathogenicity and immune response evasion67,68.
Nowadays, several studies have been used mouse mastitis models with the aim to evaluate potential anti-mastitis agents for use in dairy cows20,21,22,23,24,25,69,70,71 although the authors are aware that cow factors are relevant in the mastitis establishment28.
In order to evaluate UFV13′s in vivo activity against an MPEC strain and immune response to this treatment, an E. coli-induced mastitis model was employed. A 10-fold reduction of bacterial load was observed after viral inoculation using MOI 10 (Supplementary Figure 5).
Seven different cytokines (IL-6, TNF-α, IL-2, IFN-γ, IL-4, IL-10 and IL-17A) were measured, however only IL-10, TNF-α and IL-6 were identified as statistically significant (p < 0.05) (Fig. 6).
Studying whether T4 bacteriophage and T4-generated E. coli lysate influence cultures of peripheral blood mononuclear cells (PBMCs) activated or not by lipopolysaccharide (LPS), Bocian et al. (2016) observed that both preparations considerably increased the percentage of CD14+CD16−CD40+ and CD14+CD16−CD80 + monocytes in LPS-unactivated PBMCs cultures, as well as the concentration of IL-6 and IL-12. Notwithstanding, this result suggests that T4 bacteriophages may act both as a pro-inflammatory inducer, and also as CD40 activator, for this reason the increased expression of IL-6, IL-10, and TNF-α as a consequence of the presence of contaminating LPS left after the purification step rather than a property of the virion.
In an In vivo assessment of cytokine patterns followed by a single-dose T4 bacteriophage application in an E. coli induced mastitis mouse model, our work indicates that IL-10 levels decreased locally in phage therapeutic group when compared to the negative control (PBS buffer) (Fig. 6A). IL-10 has been extensively studied due its immunosuppressive features associated with the downregulation of pro-inflammatory cytokines, such as TNF-α and IFN-γ, and the resolution of the inflammatory process72,73. Evaluation of cytokine expression in the mammary gland in a mouse model of Streptococcus agalactiae mastitis performed by Trigo et al. (2009) revealed that the maximum concentration of IL-10 occurred after 72 hours and was correlated with a decreased level of TNF-α. This finding is in accordance with our results. Reduced IL-10 levels (Fig. 6A) and increased TNF- α (Fig. 6B) and IL-6 (Fig. 6D) is indicative of an ongoing inflammatory process. In fact, when only E. coli 30 was inoculated, an increased abundance of TNF-α was also observed when compared with the negative control. Although phage and E. coli 30 were individually able to induce an inflammatory response, phage treatment (103 PFU) did not provoke an additive effect on the production of any pro-inflammatory cytokine. Indeed, the IL-6 level diminished (Fig. 6D) after phage treatment. Different to what was found by Trigo et al. (2009), IFN-γ levels were detected in all groups, however no statistical difference was observed between the groups, which was also applicable to IL-2 (Fig. 6C–E, respectively).
The absence of detectable IL-4 and IL-17A cytokines, respectively present in Th2 and Th17 inflammatory responses, indicates that the use of T4 phage in mammary glands infected with E. coli induces Th1 T-cell responses. Th1 pattern is involved in a cellular immune response that protects against intracellular infections by viruses and microorganisms that grow in macrophages74.
Histological analysis of the mammary gland in the control group (Fig. 7A) shows intact tissue morphology described by well delimited cells and acini, as well as the presence of fatty cells and a milky secretion, even when E. coli 30 (~100 cells) were inoculated (Fig. 7B). Inoculation of the phage UFV13 (MOI 10) in the treated group led to an inflammatory reaction characterized by tissue damage, neutrophil infiltrates and mischaracterization of the acini border, although low levels of pro-inflammatory cytokines have been found as discussed previously (Fig. 7C).
Conclusion
Viral genomic analysis is a crucial step in bacteriophage screening for their use as an antibacterial agent. The UFV13 virus has no lysogenic modules or genes conferring antibiotic resistance. Viral stability analysis in the presence of detergents and organic solvents and incubation at different pHs and temperatures, have shown that the UFV13 virus presents high survivability at basic pHs and is relatively resistant under incubation at 62 °C for 40 min. vB_EcoM-UFV13 used in an animal model for mastitis reduced the total bacterial load by 90%, as well as inducing pro-inflammatory cytokines such as IL-6 and TNF-α, which makes it a potential biological agent capable of controlling acute infections caused by E. coli dairy cows.
References
Shaheen, M., Tantary, H. & Nabi, S. A Treatise on Bovine Mastitis: Disease and Disease Economics, Etiological Basis, Risk Factors, Impact on Human Health, Therapeutic Management, Prevention and Control Strategy. Adv. Dairy Res. 4, 1–10 (2016).
Thomas, V. et al. Antimicrobial susceptibility monitoring of mastitis pathogens isolated from acute cases of clinical mastitis in dairy cows across Europe: VetPath results. Int. J. Antimicrob. Agents 46, 13–20 (2015).
Zeinhom, M. M. A. & Abdel-Latef, G. K. Public health risk of some milk borne pathogens. Beni-Suef Univ. J. Basic Appl. Sci. 3, 209–215 (2014).
Crispie, F., Flynn, J., Ross, R. P., Hill, C. & Meaney, W. J. Dry cow therapy with a non-antibiotic intramammary teat seal - a review. Ir. Vet. J 57, 412 (2004).
Down, P. M., Green, M. J. & Hudson, C. D. Rate of transmission: A major determinant of the cost of clinical mastitis. J. Dairy Sci. 96, 6301–6314 (2013).
Blum, S. E., Heller, E. D., Jacoby, S., Krifucks, O. & Leitner, G. Comparison of the immune responses associated with experimental bovine mastitis caused by different strains of Escherichia coli, 190–197, https://doi.org/10.1017/S0022029917000206 (2017).
Schmelcher, M., Powell, A. M., Camp, M. J., Pohl, C. S. & Donovan, D. M. Synergistic streptococcal phage λSA2 and B30 endolysins kill streptococci in cow milk and in a mouse model of mastitis. Appl. Microbiol. Biotechnol. 99, 8475–8486 (2015).
Bryan, D., El-Shibiny, A., Hobbs, Z., Porter, J. & Kutter, E. M. Bacteriophage T4 infection of stationary phase E. coli: Life after log from a phage perspective. Front. Microbiol. 7 (2016).
Porter, J., Anderson, J., Carter, L., Donjacour, E. & Paros, M. In vitro evaluation of a novel bacteriophage cocktail as a preventative for bovine coliform mastitis. J. Dairy Sci. 99, 2053–2062 (2016).
Dias, R. S. et al. Use of phages against antibiotic-resistant Staphylococcus aureus isolated from bovine mastitis 1. J. Anim. Sci 91, 3930–3939 (2013).
Sulakvelidze, A. The challenges of bacteriophage therapy. Ind. Pharm 45, 14–18 (2011).
Gill, J. J. et al. Efficacy and Pharmacokinetics of Bacteriophage Therapy in Treatment of Subclinical Staphylococcus aureus Mastitis in Lactating Dairy Cattle. Antimicrob. Agents Chemother. 50, 2912–2918 (2006).
Merril, C. R., Scholl, D. & Adhya, S. L. The prospect for bacteriophage therapy in Western medicine. Nat. Rev. Drug Discov. 2, 489–97 (2003).
Miller, E. S. et al. Bacteriophage T4 Genome Bacteriophage T4 Genome †. Microbiol. Mol. Biol. Rev. 67, 86–156 (2003).
Chibani-Chennoufi, S., Dillmann, M. L., Marvin-Guy, L., Rami-Shojaei, S. & Brüssow, H. Lactobacillus plantarum bacteriophage LP65: A new member of the SPO1-like genus of the family myoviridae. J. Bacteriol. 186, 7069–7083 (2004).
Abouhmad, A., Mamo, G., Dishisha, T., Amin, M. A. & Hatti-Kaul, R. T4 lysozyme fused with cellulose-binding module for antimicrobial cellulosic wound dressing materials. J. Appl. Microbiol. 121, 115–125 (2016).
Rodríguez-Rubio, L., Martínez, B., Donovan, D. M., Rodríguez, A. & García, P. Bacteriophage virion-associated peptidoglycan hydrolases: potential new enzybiotics. Crit. Rev. Microbiol. 39, 427–434 (2013).
Bruttin, A., Brüssow, H. & Bru, H. Human Volunteers Receiving Escherichia coli Phage T4 Orally: a Safety Test of Phage Therapy. Antimicrob. Agents Chemother. 49, 2874–2878 (2005).
Bocian, K. et al. LPS-activated monocytes are unresponsive to T4 phage and T4-generated Escherichia coli lysate. Front. Microbiol. 7 (2016).
Olson, M. A., Siebach, T. W., Griffitts, J. S., Wilson, E. & Erickson, D. L. Genome-wide identification of fitness factors in mastitisassociated Escherichia coli. Appl. Environ. Microbiol. 84 (2018).
Wang, J. et al. Oligopeptide Targeting Sortase A as Potential Anti-infective Therapy for Staphylococcus aureus. Front. Microbiol. 9, 1–10 (2018).
Iwano, H. et al. Bacteriophage ΦSA012 Has a Broad Host Range against Staphylococcus aureus and Effective Lytic Capacity in a Mouse Mastitis Model. Biology (Basel) 7, 8 (2018).
Hu, G. et al. Cynatratoside-C from Cynanchum atratum displays anti-inflammatory effect via suppressing TLR4 mediated NF-κB and MAPK signaling pathways in LPS-induced mastitis in mice. Chem. Biol. Interact. 279, 187–195 (2018).
Roussel, P. et al. Escherichia coli mastitis strains: In vitro phenotypes and severity of infection in vivo. Plos One 12, 1–20 (2017).
Johnzon, C. F. et al. Mastitis pathogens with high virulence in a mouse model produce a distinct cytokine profile in vivo. Front. Immunol. 7, 1–11 (2016).
Yu, Y. et al. Efficacy of cefquinome against Escherichia coli environmental mastitis assessed by pharmacokinetic and pharmacodynamic integration in lactating mouse model. Front. Microbiol. 8, 1–9 (2017).
Ingman, W. V., Glynn, D. J. & Hutchinson, M. R. Mouse models of mastitis – how physiological are they? Int. Breastfeed. J. 10, 12 (2015).
Burvenich, C., Van Merris, V., Mehrzad, J., Diez-Fraile, A. & Duchateau, L. Severity of E. coli mastitis is mainly determined by cow factors. Veterinary Research 34, 521–564 (2003).
Sambrook, J. & Russell, D. W. Molecular Cloning - Sambrook & Russel - Vol. 1, 2, 3. Cold Springs Harb. Lab. Press 3th Editio (2001).
Adams, M. Bacteriophages. Bacteriophages 620 (1959).
Duarte, V. S. et al. Complete genome sequence of vB_EcoM-UFV13, a new bacteriophage able to disrupt Trueperella pyogenes biofilm. Genome Announc. 4 (2016).
Gasteiger, E. et al. ExPASy: The proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 31, 3784–3788 (2003).
Lowe, T. M. & Eddy, S. R. TRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 25, 955–964 (1996).
Naville, M., Ghuillot-Gaudeffroy, A., Marchais, A. & Gautheret, D. ARNold: a web tool for the prediction of Rho-independent transcription terminators. RNA Biology 8, 11–13 (2011).
Klucar, L., Stano, M. & Hajduk, M. phiSITE: database of gene regulation in bacteriophages. Nucleic Acids Res. 38, D366–D370 (2010).
Grant, J. R. & Stothard, P. The CGView Server: a comparative genomics tool for circular genomes. Nucleic Acids Res. 36 (2008).
Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).
Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359 (2012).
Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).
Cingolani, P. et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin). 6, 80–92 (2012).
Darling, A. C. E., Mau, B., Blattner, F. R. & Perna, N. T. Mauve: Multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 14, 1394–1403 (2004).
Treu, L. et al. The impact of genomic variability on gene expression in environmental Saccharomyces cerevisiae strains. Environ. Microbiol. 16, 1378–1397 (2014).
Tuimala, J. A primer to phylogenetic analysis using the PHYLIP package. Espoo Finland Center for Scientific Computing Ltd 6 (2006).
Huson, D. H. et al. Dendroscope: An interactive viewer for large phylogenetic trees. BMC Bioinformatics 8, 460 (2007).
Larkin, M. A. et al. Clustal W and Clustal X version 2.0. Bioinformatics 23, 2947–2948 (2007).
Katoh, K. & Standley, D. M. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol. Biol. Evol. 30, 772–780 (2013).
Kumar, S., Stecher, G. & Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol. Biol. Evol. msw054, https://doi.org/10.1093/molbev/msw054 (2016).
Villarroel, J. et al. HostPhinder: A phage host prediction tool. Viruses 8 (2016).
Green, M. R. & Sambrook, J. Molecular Cloning: A Laboratory Manual (Fourth Edition). Cold Spring Harbor Laboratory Press; 4th edition (June 15, 2012) 2028, Available at, http://www.amazon.com/Molecular-Cloning-Laboratory-Edition-Three/dp/1936113422 (2012).
Shevchenko, A., Tomas, H., Havlis, J., Olsen, J. V. & Mann, M. In-gel digestion for mass spectrometric characterization of proteins and proteomes. Nat. Protoc. 1, 2856–2860 (2007).
Jurczak-Kurek, A. et al. Biodiversity of bacteriophages: morphological and biological properties of a large group of phages isolated from urban sewage. Sci. Rep. 6, 34338 (2016).
Déziel, E., Comeau, Y. & Villemur, R. Initiation of biofilm formation by Pseudomonas aeruginosa 57RP correlates with emergence of hyperpiliated and highly adherent phenotypic variants deficient in swimming, swarming, and twitching motilities. J. Bacteriol. 183, 1195–1204 (2001).
Belgini, D. R. B. et al. Culturable bacterial diversity from a feed water of a reverse osmosis system, evaluation of biofilm formation and biocontrol using phages. World J. Microbiol. Biotechnol. 30, 2689–2700 (2014).
CLSI. Performance Standards for Antimicrobial Susceptibility Testing; Twenty-Seventh Edition. Clinical and Laboratory Standards Institute (2017).
Chandler, R. L. Studies on experimental mouse mastitis relative to the assessment of pharmaceutical substances. J. Comp. Pathol. 81, 507–514 (1971).
Naghili, H. et al. Validation of drop plate technique for bacterial enumeration by parametric and nonparametric tests. Vet. Res. forum an Int. Q. J. 4, 179–83 (2013).
Ackermann, H. W. Bacteriophage observations and evolution. Research in Microbiology 154, 245–251 (2003).
Belle, A., Landthaler, M. & Shub, D. A. Intronless homing: Site-specific endonuclease SegF of bacteriophage T4 mediates localized marker exclusion analogous to homing endonucleases of group I introns. Genes Dev. 16, 351–362 (2002).
Comeau, A. M., Arbiol, C. & Krisch, H. M. Composite conserved promoter-terminator motifs (PeSLs) that mediate modular shuffling in the diverse T4-like myoviruses. Genome Biol. Evol 6, 1611–1619 (2014).
Jun, J. W. et al. Bacteriophage application to control the contaminated water with Shigella. Sci. Rep. 6, 22636 (2016).
Hinton, D. M. Transcriptional control in the prereplicative phase of T4 development. Virol. J. 7, 289 (2010).
Leiman, P. G. et al. Morphogenesis of the T4 tail and tail fibers. Virol. J. 7, 355 (2010).
Kramberger, P., Urbas, L. & Štrancar, A. Downstream processing and chromatography based analytical methods for production of vaccines, gene therapy vectors, and bacteriophages. Human Vaccines and Immunotherapeutics 11, 1010–1021 (2015).
Thung, T. Y. et al. Use of a lytic bacteriophage to control Salmonella Enteritidis in retail food. LWT - Food Sci. Technol. 78, 222–225 (2017).
Ogola, H., Shitandi, A. & Nanua, J. Effect of mastitis on raw milk compositional quality. J. Vet. Sci. 8, 237–242 (2007).
Zhanel, G. G., Zhanel, M. A. & Karlowsky, J. A. Oral Fosfomycin for the Treatment of Acute and Chronic Bacterial Prostatitis Caused by Multidrug-Resistant Escherichia coli. Can. J. Infect. Dis. Med. Microbiol. 2018, 1–9 (2018).
Domenech, M., Ramos-Sevillano, E., García, E., Moscoso, M. & Yuste, J. Biofilm formation avoids complement immunity and phagocytosis of Streptococcus pneumoniae. Infect. Immun. 81, 2606–2615 (2013).
Kao, C. Y. et al. The complex interplay among bacterial motility and virulence factors in different Escherichia coli infections. Eur. J. Clin. Microbiol. Infect. Dis. 33, 2157–2162 (2014).
Brouillette, E., Grondin, G., Talbot, B. G. & Malouin, F. Inflammatory cell infiltration as an indicator of Staphylococcus aureus infection and therapeutic efficacy in experimental mouse mastitis. Vet. Immunol. Immunopathol. 104, 163–169 (2005).
Brouillette, E. & Malouin, F. The pathogenesis and control of Staphylococcus aureus-induced mastitis: Study models in the mouse. Microbes Infect. 7, 560–568 (2005).
Nazemi, S. et al. Expression of acute phase proteins and inflammatory cytokines in mouse mammary gland following Staphylococcus aureus challenge and in response to milk accumulation. J. Dairy Res. 81, 445–54 (2014).
Trigo, G. et al. Leukocyte populations and cytokine expression in the mammary gland in a mouse model of Streptococcus agalactiae mastitis. J. Med. Microbiol. 58, 951–958 (2009).
Masso-Welch, P. A., Merhige, P. M., Veeranki, O. L. M. & Kuo, S.-M. Loss of IL-10 Decreases Mouse Postpubertal Mammary Gland Development in the Absence of Inflammation. Immunol. Invest. 41, 521–537 (2012).
Kaiko, G. E., Horvat, J. C., Beagley, K. W. & Hansbro, P. M. Immunological decision-making: How does the immune system decide to mount a helper T-cell response? Immunology 123, 326–338 (2008).
Acknowledgements
We are grateful to the Núcleo de Análise de Biomoléculas and Núcleo de Microscopia e Microanálise of the Universidade Federal de Viçosa for providing the facilities to conduct the experiments. We also acknowledge the financial support of the following Brazilian agencies: Fundação de Amparo à Pesquisa do Estado de Minas Gerais (Fapemig), Coordenacão de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Financiadora de Estudos e Projetos (Finep), Sistema Nacional de Laboratórios em Nanotecnologias (SisNANO)/Ministério da ciência, tecnologia e Informação (MCTI). We are also grateful to Empresa Brasileira de Pesquisa Agropecuária (EMBRAPA), Gado de Leite Farm, Juiz de Fora, Minas Gerais, Brazil, that kindly provided us with Escherichia coli 30.
Author information
Authors and Affiliations
Contributions
V.S.D., R.S.D. and S.O.P. conceived and designed the experiments; V.S.D., A.M.K., S.C., L.T., A.X., C.G.F. and P.M.P.V. performed genomic analysis; C.S. contributed to E. coli-induced mastitis mouse model; M.S.V., I.S.P., and G.R.S. conducted phage physiological features assays. F.M. took part in the histological analysis; J.S.C. was involved in protein profile experiments; C.C.S. contributed materials. V.S.D. wrote the main manuscript text. All the authors have read and contributed to the final version of the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
da Silva Duarte, V., Dias, R.S., Kropinski, A.M. et al. Genomic analysis and immune response in a murine mastitis model of vB_EcoM-UFV13, a potential biocontrol agent for use in dairy cows. Sci Rep 8, 6845 (2018). https://doi.org/10.1038/s41598-018-24896-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-018-24896-w
This article is cited by
-
Genomic and proteomic characterization of vB_SauM-UFV_DC4, a novel Staphylococcus jumbo phage
Applied Microbiology and Biotechnology (2023)
-
Milk microbial composition of Brazilian dairy cows entering the dry period and genomic comparison between Staphylococcus aureus strains susceptible to the bacteriophage vB_SauM-UFV_DC4
Scientific Reports (2020)
-
Evaluation of phage therapy in the treatment of Staphylococcus aureus-induced mastitis in mice
Folia Microbiologica (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.